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SPECIAL  FEATURES  OF  THE  LIMITATIONS  TO  CHINA'S  PATENT  LAW 

Beijing  GUOJI  MAOYI  WENTI  [INTERNATIONAL  TRADE  JOURNAL]  in  Chinese  No  3, 
May-Jun  85  pp  1-7 

[Article  by  Zheng  Chengsi  [6774  2052  1835],  Law  Institute,  the  Chinese  Academy 
of  Social  Sciences:  "The  Limitations  of  Patents  With  Regard  to  International 
Trade" ] 

[Excerpt]  II.  Some  Special  Features  of  the  Limitations  to  Chinese  Patents 
A.  The  "State  Plan  Permit"  System 

In  order  to  assist  in  the  popularization  of  advanced  technologies  in  China  and 
at  the  same  time  not  to  hinder  the  importation  of  foreign  advanced 
technologies,  China's  patent  law  has  adopted  a  unique  "state  plan  permit" 
system.  This  system  works  as  follows:  Relevant  departments  in  charge  under 
the  State  Council  and  the  people's  governments  of  provinces,  autonomous 
regions  and  centrally-administered  municipalities  are  authorized  to  allow, 
according  to  the  state  plan,  designated  units  to  use  the  patent  for  an 
important  invention  or  creation  held  by  units  under  state  ownership  that  are 
within  their  own  system  or  under  their  jurisdiction;  the  designated  unit  that 
uses  the  patent  shall  pay  royalties  to  the  unit  which  holds  the  patent,  in 
accordance  with  state  stipulations.  If  a  Chinese  unit  under  the  collective 
ownership  system  or  an  individual  owns  a  patent  which  is  of  major  significance 
to  the  national  or  public  interests,  and  which  should  be  popularized,  the 
relevant  departments  in  charge  under  the  State  Council  may  submit  a  report  to 
the  State  Council  and  on  approval  by  the  State  Council,  this  type  of  patent 
may  also  be  dealt  with  according  to  the  above-stated  provisions. 

To  implement  the  "state  plan  permit"  system  and  other  special  restrictions  of 
the  Chinese  system,  the  Chinese  patent  law  divides  patentees  into  three 
categories:  1.  Patent-holders,  i.e.  patentees  who  are  Chinese  enterprises  and 
units  under  state  ownership;  2.  Chinese  patent-owners,  i.e.  patentees  who  are 
Chinese  enterprises  and  units  under  the  collective  ownership  system  and 
Chinese  natural  persons;  3.  Foreign  owners  of  patents,  i.e.  Chinese-foreign 
joint  venture  enterprises,  wholly-owned  foreign  enterprises  and  foreign 
natural  persons  who  have  obtained  Chinese  patents. 
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The  "state  plan  permit"  system  applies  only  to  patentees  under  the  first  two 
categories.  Foreign  owners  of  advanced  technologies  need  not  worry  about 
having  their  patents,  brought  in  as  shares  in  joint  ventures,  as  investments, 
etc.,  "popularized"  and  their  exclusive  use  rights  invalidated. 

In  relation  to  the  State  Council  and  its  relevant  departments  in  charge  and 
the  relevant  departments  of  provinces,  autonomous  regions  and  centrally- 
administered  municipalities,  Chinese  holders  and  owners  of  patents  are  by  law 
restricted  in  their  exclusive  use  rights.  These  particular  restrictions  are 
different  from  most  of  the  restrictions  discussed  in  Chapter  1.  Under  most  of 
the  conditions  discussed  in  Chapter  1,  the  exclusive  use  right  of  the  patentee 
loses  all  binding  effect  on  everyone,  but  under  the  "state  plan  permit" 
system,  the  exclusive  use  right  of  the  patentee  loses  validity  only  with 
respect  to  the  state  organs  concerned  and  the  users  specially  designated  by 
these  organs.  This  means  that  the  holder  or  owner  of  a  patent  that  is  to  be 
"popularized"  still  retains  his  exclusive  rights  in  relation  to  all 
enterprises  and  individuals  other  than  those  specially  designated  by  the  state 
for  use  of  the  patent.  It  means  that  the  enterprise  or  person  has  the  right  of 
granting  licenses  and  collecting  use  fees,  as  well  as  the  right  to  file  suits 
against  unlicensed  users  and  to  claim  compensation. 

B.  Restrictions  in  the  Assignment  of  Patents 

The  right  of  assignment  is  not  a  right  peculiar  to  owners  of  patents,  but  one 
enjoyed  (under  general  conditions)  by  every  owner  of  a  property.  In  most 
countries,  a  patentee  is  not  at  all  restricted  in  the  assignment  of  his  patent 
right,  just  as  the  owner  of  goods  when  selling  his  goods.  However,  because 
the  "patent-holders"  among  the  Chinese  patentees  also  have  ownership  rights  in 
units  under  the  system  of  state  ownership,  they  are  not  legally  designated 
"owners"  of  the  patent,  but  merely  "holders."  When  assigning  patent  rights, 
this  type  of  unit  must,  therefore,  obtain  approval  from  its  higher  competent 
authority.  In  international  trade,  foreign  enterprises  that  claim  to  have 
obtained  Chinese  patents  must  pay  special  attention  to  this  matter.  They  must 
first  of  all  clarify  whether  the  patent  which  is  being  assigned  to  them  is 
held  by  a  unit  under  state  ownership,  and  if  so,  whether  the  transfer  was 
approved  by  its  superior  organ.  Otherwise,  they  would  possibly  be  obtaining  a 
patent  through  an  illegal  transfer.  At  the  same  time,  any  assignment  of  a 
patent  to  a  foreigner  must  be  approved  by  the  responsible  State  Council 
department.  This  provision  applies  to  the  patent  rights  of  all  Chinese 
juristic  persons  and  individuals. 

The  assignments  dealt  with  here  refer  to  the  transfer  of  patent  rights  as 
property  and  not  to  the  transfer  of  use  rights  (i.e.,  not  to  the  granting  of 
patent  licenses).  The  Chinese  patent  law  imposes  no  restrictions  at  all  on 
the  right  of  Chinese  holders  or  owners  of  patents  to  issue  licenses. 

C.  No  "Import  Rights"  Being  Conferred 

China's  patent  law  does  not  confer  "import  rights"  to  the  patentee;  this  is  a 
further  provision  that  limits  the  patent.  The  question  whether  an  import 
right  is  conferred  is  usually  brought  up  in  the  patent  legislation  of 
developing  countries.  The  granting  of  an  import  right  has  its  positive 
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aspects  with  regard  to  the  importation  of  advanced  technologies.  Because  the 
owner  of  the  advanced  technology  will  in  general  apply  for  a  patent  in 
developed  countries  as  well  as  in  developing  countries,  an  enterprise  that  has 
obtained  a  patent  license  in  a  developed  country  has  frequently  a  much  larger 
production  capacity  at  its  disposal  than  an  equivalent  enterprise  in  a 
developing  country.  If  the  patentee  in  the  latter  kind  of  country  cannot 
enjoy  the  controlling  right  over  the  importation  of  the  patented  product,  the 
assignee  in  the  first  type  of  country,  through  marketing  the  products  in  the 
latter  type,  can  make  it  unprofitable  to  manufacture  in  the  latter  type  of 
country.  In  consideration  of  this  consequence,  the  owner  of  the  technology 
will  abstain  from  applying  for  a  patent  in  the  latter  type  of  country,  but  if 
import  rights  would  be  conferred,  it  would  induce  the  foreign  patentee  to 
obstruct  the  importation  of  the  same  article,  manufactured  in  high  quality  and 
at  low  price,  into  the  developing  country  and  thereby  cause  economic  losses  to 
these  countries. 

However,  we  must  come  to  understands  countries  which  confer  import  rights  to 
the  patentee  will  also  frequently  not  allow  this  import  right  to  be  exercised 
at  will  by  the  patentee.  First,  the  import  right  is  not  being  granted  for 
designs;  second,  as  to  the  import  right  for  invention  patents,  it  is 
frequently  restricted  by  other  even  more  stringent  measures.  For  instance, 
Article  1,  Paragraph  2  of  the  Unified  Patent  Law  of  the  African  Intellectual 
Property  Organization  mentions  the  conferring  of  import  rights  to  the 
patentees,  but  also  stipulates  in  its  Article  55,  Paragraph  1:  For  the 
importation  of  certain  patented  products  for  the  important  need  of  national 
defense,  public  health  and  national  economy,  the  industrial  property 
department  of  any  member  nation  of  the  African  Intellectual  Property 
Organization  is  authorized  to  issue  "requisitioning  licenses,"  to  assist 
importation.  The  requisitioning  by  the  state  is  one  method  that  inflicts  on 
the  owner  of  the  technology  a  more  serious  restriction  than  not  having  the 
import  right.  China  did  not  institute  a  state  requisitioning  system,  rights 
are  therefore  balanced,  and  it  is  only  that  no  import  right  is  conferred, 
which  will  not  at  all  hinder  Chinese  importation  of  technologies. 

In  developed  countries  where  the  patent  law  confers  import  rights  to 
patentees,  restrictive  provisions  are  frequently  added  to  prevent  the  patentee 
using  this  monopolistic  right  from  creating  a  negative  influence  on  the  market 
of  one’s  own  country.  For  instance,  Article  48  of  the  British  patent  act  of 
1977  (revised  in  1979)  stipulates!  If  a  certain  patented  invention  is  capable 
of  being  used  in  the  United  Kingdom  but  has  its  use  hindered  by  importation  of 
the  said  patented  product,  anyone  expecting  to  use  that  invention  may  apply 
for  a  compulsory  license;  in  another  respect,  if  the  demand  on  the  British 
market  for  a  certain  patented  product  is  being  met  to  a  substantial  extent  by 
importation,  the  Patent  Office  can  change  this  situation  by  issuing  compulsory 
patent  licenses,  so  that  one’s  own  country  may  be  able  to  manufacture  the  same 
type  of  products. 

This  shows  that  not  granting  import  rights  to  patentees,  although  a  limitation 
of  a  special  kind  in  China,  is  not  at  all  excessive  compared  with  other 
countries;  it  is  a  reasonable  limitation. 
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D.  "No  Punishment  for  Ignorance" 


Article  62,  Paragraph  2,  of  the  Chinese  Patent  Law  stipulates:  "The  use  or 
sale  of  a  patented  product  without  any  knowledge  that  such  product  was 
manufactured  and  sold  without  the  permission  of  the  patentee"  shall  not  be 
deemed  to  infringe  a  patent.  Because  use  was  unknowing  use,  it  can  generally 
be  called  "use  in  good  faith".  Using  the  product  well  knowing  that  the 
article  is  patented  by  someone  else  can  be  called  "use  in  bad  faith."  In  the 
patent  laws  of  most  countries  such  "use  in  bad  faith"  is  considered  an  act  of 
patent  infringement,  but  they  do  not  explicitly  say  whether  "use  in  good 
faith"  constitutes  a  patent  infringement,  or  make  "use  in  good  faith"  into  a 
limitation  of  the  patent  law.  This  article  of  China’s  patent  law,  therefore, 
appears  to  be  unique. 

The  said  provision  in  the  Chinese  law  can  be  similarly  expressed  by  a 
traditional  saying  in  the  ancient  Chinese  legal  system,  namely  "no  punishment 
in  case  of  ignorance."  The  main  reason  for  this  provision  is,  in  my  opinion, 
that  under  the  present  conditions  in  China,  when  protection  of  patents  is  only 
in  its  initial  stage,  and  the  majority  of  Chinese,  at  this  stage,  have  no 
direct  understanding  of  the  actual  significance  of  the  monopoly  of  the 
patentee,  it  is  hard  to  prevent  patented  products  from  being  unknowingly  used 
or  sold  by  careless  oversight.  It  is,  therefore,  reasonable  not  to  treat  such 
actions  as  patent  infringements  at  the  present  time.  However,  I  believe  that 
the  specific  application  of  this  provision  will  at  least  require  some 
supplementation:  First,  on  gaining  the  knowledge  which  he  lacked  before 
(regardless  of  whether  he  was  warned  by  the  patentee  or  informed  by  others), 
the  offender  must  immediately  halt  the  activity  in  question.  Second,  if  the 
offending  party  has  gained  an  unjustifiable  profit  from  the  patented  products 
during  the  time  he  unknowingly  used  or  sold  them,  he  must,  on  gaining 
knowledge  of  the  state  of  things,  turn  the  profits  over  to  the  state  treasury 
or  use  them  to  compensate  the  patentee  for  losses  (if  there  had  been  losses). 
Third,  the  offending  party  must  submit  positive  evidence  of  his  "ignorance." 

At  the  same  time,  attention  should  also  be  paid  to  the  following:  many 
foreign  countries  merely  recognize  responsibility  for  "acts  committed  with 
full  knowledge,"  but  do  not  explicitly  negate  responsibility  for  "acts 
committed  in  ignorance."  If  the  commercial  transactions  that  possibly  relate 
to  patents  took  place  outside  of  China,  the  special  provisions  of  Article  62, 
Paragraph  2,  of  the  Chinese  patent  law  must  not,  therefore,  be 
indiscriminately  applied,  but  every  possible  effort  should  be  made  beforehand 
to  avoid  such  unknowing  use  or  sale,  as  otherwise  the  offender  will  possibly 
be  sued  (and,  furthermore,  adjudged)  as  having  committed  an  infringement  of 
patent  rights. 
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FOETUS  MODEL,  DEMOGRAPHIC  STABILITY  DESCRIBED 

Beijing  SHUXUE  DE  SHIJIAN  YU  RENSHI  [MATHEMATICS  IN  PRACTICE  AND  THEORY] 
in  Chinese  No  4,  Oct  84  pp  31-35 

[Article  by  Shi  Jiahung  [2508  0857  7703]  of  Fu  Dan  University:  "Foetus  Model 
and  Its  Stability  in  Demography"] 

[Text]  I.  Introduction 

The  huge-numbered  population  problem  is  often  induced  by  solving  the  initial 
boundary  value  problem  of  a  partial  differential  equation  of  the  first  order 
[1,  2].  Without  considering  population  migration,  the  partial  differential 
equation  is 


«»  +  «.  —  —/*(*><)*>  (1) 

In  the  equation,  the  unknown  function  u(x,  t)  represents  the  age  density 
of  the  population  at  time  "t"  and  age  "x",  and  p(x,  t)  represents  the  ratio 
of  the  death  rate  of  the  population.  Variables  "x"  and  t  represent, 
respectively,  age  and  time  with  units  of  age  and  year,  respectively. 

The  age  density  of  population  at  the  initial  time  is  known,  that  is, 

».  .  •  .. 

u(x,0 )  —  (2) 

Marginal  conditions  are  in  a  feedback  form  [3]: 

«(0,/)  -  ‘ h(x,t)1t(x,/)u(x,0dx,  ( 3 ) 

j^t ,  r  ,  r9  are  the  upper  and  lower  limits  of  the  female  child-bearing  age, 
h(x,  t)  Is  the  birth  model,  k(x,  t)  is  the  ratio  of  the  female  population 
and  (t)  is  the  female  birth  rate  at  time  "t". 

Here  in  birth  rate  there  is  no  distinction  between  the  number  of  times  of 
birth  and  the  number  of  foetuses,  and  for  a  more  precise  and  reasonable 
calculation,  it  is  appropriate  to  consider  separately  the  birth  rate  of 
each  birth  because  the  birth  model  of  each  birth  is  obviously  different. 
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The  birth  model  curve  of  i+t  is  almost  the  birth  curve  of  i  moving 
horizontally  several  years  to  the  right.  On  the  other  hand,  from  the 
practical  point  of  view,  the  foetus  model  is  also  convenient.  The  reason 
is  that  in  the  formulation  of  a  population  policy,  the  relationship  between 
the  birth  rate  and  the  birth  model  can  be  better  considered. 

Suppose  that  the  birth  model  of  "i"  birth  is  h.(x,  t)  and  the  birth  rate 
is  MO  •  Then  the  foetus  model  is  to  change  Equation  (3)  to 

«(o,  0  —  2  MO  f  *  M*>04(*>0«(*>0<te>  o  <  p,  <  l.  (4) 

i-j  .  . 


What  value  positive  integer  "N"  takes  depends  on  the  problem.  Based  on  the 
situation  of  our  nation,  "N"  may  take  the  value  of  3.  However,  whatever  value 
"N"  may  take  will  not  have  any  substantial  effect  on  our  following  conclusion. 
Therefore,  we  still  use  the  general  "N"  to  express  the  marginal  condition. 

Articles  [2],  [3]  and  [5]  on  the  stability  problem  of  systems  (1),  (2)  and 
(3)  under  the  Lyapunov  stability  function  have  been  discussed  in  detail, 
and  it  has  been  proved  that  there  exists  a  marginal  birth  rate  fie  ,  in  which 
when  constant  fi  <  fit  »  the  solution  of  (1),  (2)  and  (3)  is  stable,  and  when 
fi>  fit  ,  it  is  not. 

Now  let  us  discuss  the  solution  of  (1) ,  (2)  and  (4)  and  the  problem  of 
stability. 

II.  Solution  of  Problem  of  Initial  Marginal  Value 

Theorem  1:  When  [0,  rM),  0,  the  solution  for  (1),  (2)  and  (4)  is 


<p(x  —  0*  ,  x  —  t  >  0 


Proof:  Let 


2  M*  —  —  X 

i=  t  Jr> 


x)%(r ,l  — *  x)u  (r ,t  —  x)dre  f 


x  —  t<  0 


(5) 


then  function  v(x,  t)  satisfies 


v,  +  vx  0,  (6) 

v(x,  0)  —  q>(.x)e^°MSUt ,  (7) 

v(0 ,/)  —  2  v( r,t)dr  (7') 
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The  general  solution  representation  of  (6)  is: 


o(x,t )  —  f(x  —  t)» 


To  get  the  form  of  function  f,  let  t-0,  and  we  get,: 


/v*)  —  ^  °* 


When  *  —  /  >  0 


K*>0  —  <p(*  —  t)ei0  , 

—  v{x,  t)e  J0  —  <p(x  —  *)* 


Substituting  x=0  for  f(x  -  t) ,  we  get 


/(-,)  -  £>,(,)  (\(r »<)«»•  Kr,/Vr,  /  >0. 

i  =  1  ^  r' 


When  *  —  <  <  0  , 


,(«,,)  -  K,  - ,)  -  ±  m<  -  *)jy  r,«  -  *>«'.«  -  *<'•«  -  *>*■ 

-  V  8,0  -  x)  P*  h,(r,t  -  *)*0,/  -  x)«(r,t  -  x)dr . 

ir' 

X  (r,A,(r,<  -  *)40>  -  x)«(r,/  -  C 

■''l 


The  combination  of  (8)  and  (9)  is  (5),  and  the  theorem  is  proved.  Notice 
that  n  > 0  ,  so  (5)  is  a  recurrence  representation  of  "u". 

III.  Stability  of  Solution 

Here  we  discuss  the  continuity  of  the  solution's  dependence  on  the  t^O 
initial  state.  To  describe  it  precisely,  that  is,  if  e>0  ,  there  is  always 


S  >  0  ,  then  when 


||<pOOII  <  S 
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there  will  be 


sup  ||k(*,/)||  <  s,  (11) 

/>0 

Among  them  ||*||  is  a  norm  of x€.  I0,ru]  with  a  certain  meaning.  Then  the  system 
is  considered  stable.  Otherwise,  it  is  not  stable,  especially  in  a  case 
no  matter  how  small  S  is,  there  is  always  a  <p(»)  to  satisfy  (10)  to  make 

lim  sup  \\u(x,t)\\  —  oo, 
r~«  ««t 


Then  the  system  (1) (2) (4)  is  divergent  unstable. 

Notice  that  u(x,t)  is  the  continuous  function  of  *€  [0,  ry],  t  ^  0. 

Here  we  take  the  maximum  value  of  the  norm  r€  [0,rM],that  is, 

H«(;r,/)||  —  max  u(x,t). 

x€[0tfy) 

N  (t  fg 

Theorem  2:  i)  If  fiiO)  1  k/(x,i)Ji(x,t)e  .^ ****** dx^  1,  (12) 

J'i  . 

then  the  system  is  stable. 

ii)  Suppose  (.there  exists  n>0  ,  that  makes 

V  pr,(r)  I  *  h,(x,t)^(x,t)e  i9**'**  dx  >  1  +  n»  O3) 

■  *  1  •''< 

Then  the  system  (1) (2) (4)  is  divergent  unstable. 

Proof:  Suppose  (12)  stands,  write 

max  ©(*>  —  M. 

#€  (• pfti 

Let  /G[0,rt],  From  (8)  we  can  easily  see  that  when  *€  {t,  fj] 

mu  .<*,/)  -  »«  »(»  -  ,),  Ir-’M  <  M  04)  • 

When  *€  from  (5) ,  (12)  and  (14)  we  get 

*'(*>  0  —  y.  £/(<  “  *)  f  * hi(r,i  —  *M(r,/  —  x~)v(r ,t  —  x)*"^*5*  dr 

#-i  J't 

“  *)  (  —  *)$(M  —  x)*“^‘<6‘<*  rfr  <  M.  (15) 

/-1 

so  against  /€  [0,rt]  stands 


sup  v(x,  /)  *^  A/.  (16) 

x«[0^,3 
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Repeat  the  steps  in  (15),  and  it  will  be  easy  to  get  (16)  when  . 

Then  to  the  assigned  g  >  0  ,  we  only  take 

>  -  . 


Then  from  (16)  we  know  that  (1)  can  be  deduced  from  (10).  Theorem  (i)  is 
proved. 


Now  suppose  (13)  stands,  if  we  let 


then  when  *  —  >  ^  0  > 

u(x,t')  *■  a«. 

Let  us  use  finite  induction  to  prove  that  when 

—  (m  4-  1>S  <  *  —  t  <  —  mr2 
v(_x ,/)  ^  (1  + 

In  fact ,  when  —  /  s^O 


(17) 


v(x,0 


■  ■ 

<h  —  *) 

i  m  1  J 


x)A(r,/  —  x)v(r,t  —  *)<?  ^  ^ 

.  *)*<>,,  -  ir  3*  (1+  l)«o. 


(18) 


If  the  second  equality  sign  in  (18)  is  changed  to  then  (18)  stands 

against  — r,  <  *.—  0  and  this  is  (17)  when  m=0. 

Repeat  steps  in  (18)  and  change  aQ  to  (1  +  v)m<h  ,  and  using  finite  induction 
(17)  can  be  proved. 

From  (17)  the  conclusion  of  Theorem  (ii)  can  be  immediately  derived. 

The  theorem  is  proved. 

Note  1:  If  (12)  or  (13)  does  not  stand  at  any  condition  but  only  when  /  > /#  , 
and  tQ  is  a  positive  constant,  then  obviously  the  conclusion  of  the  theorem 
will  still  stand. 

Note  2:  We  call  the  situation  in  which  the  equality  sign  in  (12)  stands 
a  birth  marginal  saturation.  If  n >  Pn  hi  and  k  have  no  relation  with  time 
"t",  then  the  marginal  saturation  becomes 

'  Sftf  "!•  '  (19) 

iml  J'l 
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We  call  the  birth  rate  PuPu' ’ '  >fa  that  satisfies  the  relational  formula  (19) 
a  marginal  birth  rate.  In  this  situation,  the  (ii)  of  Theorem  2  can  be 
changed  to  (ii) '  on  condition 

im  %  J'l  l 

Then  the  system  is  divergent  stable.  Therefore,  in  this  situation, 
when  the  birth  rate  is  greater  than  the  marginal  value,  it  becomes  unstable, 
and  when  the  birth  rate  is  less  or  equal  to  the  marginal  value,  it  becomes 
stable. 

Note  3:  If  we  let 

N 

A/Or,/)  -A(»),  »- 1,2,  £>,  -/?,  (20) 

then  (19)  changes  to 

jr,  •  ;  •  •  • 

This  is  just  the  representation  of  the  marginal  birth  rate  mentioned  in 
Articles  [3]  arid  [5].  The  first  formula  in  (20)  represents  the  birth  model 
without  foetus  differentiation,  and  the  in  Formula  2  is  the  birth  rate 
defined  in  Article  [3].  Therefore,  Theorem  2  is  the  extension  of  the 
stability  concluded  in  Articles  [3]  and  [5]. 
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APPLICATION  OF  FUZZY  MATHEMATICS  IN  CORE  DRILLING  DISCUSSED 

Beijing  SHUXUE  DE  SHIJIAN  YU  RENSHI  in  Chinese  No  1,  Jan  85  pp  14-18 

[Article  by  Wang  Shaohao  [3769  1421  6275]  and  Gao  Shen  [7559  2773],  Zhongnan 
Institute  of  Metallurgy] 

[Text]  Most  of  the  studies  in  petrographic  grading  by  the  geological  explora¬ 
tion  departments  in  China  and  foreign  countries  are  limited  to  studies  of  the 
drillability  of  rocks  which  offer  various  programs  in  rock  grading  according  to 
drillability.  This  kind  of  grading  gives  consideration  only  to  the  degree  of 
difficulty  of  the  mechanical  disaggregation  of  the  rocks  without  considering 
the  effects  of  the  stratigraphic  complexity,  wall  stability,  hole  deviation 
and  difficulty  of  core  sampling  on  the  drilling  technology  process  as  a  whole. 
Therefore,  this  kind  of  grading  can  reflect  the  drillability  of  the  rocks  but 
cannot  reflect  the  difficulty  of  the  core  drilling  of  a  certain  area  as  a 
whole. 

In  estimating  investment  needed  to  prospect  in  a  certain  area,  the  management 
department  must  make  a  comprehensive  evaluation  of  the  difficulties  of  drilling. 
We  believe  that  besides  drillability,  one  must  also  consider  factors  of  sta¬ 
bility,  difficulty  of  core  sampling  and  hole  deviation.  Only  with  the  four 
factors  considered  can  one  have  an  overall  view  of  the  difficulties  of  the 
drilling  operation. 

Lacking  a  more  scientific  method  in  our  country  for  the  calculation  of  invest¬ 
ment  needed  to  prospect  a  certain  area,  we  present  for  discussion  in  this 
article  a  comprehensive  evaluation  method  using  fuzzy  mathematics,  with  con¬ 
sideration  of  the  four  factors  of  drillability,  stability,  difficulty  of  core 
sampling  and  hole  deviation  in  a  certain  site  to  provide  a  comprehensive 
evaluation  and  a  classification  of  the  difficulty  of  an  exploration  drilling. 

We  also  offer  a  simple  model  of  calculation  for  the  evaluation  of  exploration 
drilling  investment  in  a  certain  area. 

Let  the  aggregate  of  the  target  factors  be 

U  =  (drillability,  stability,  difficulty  of 
core  sampling,  hole  inclination) 
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The  result  of  the  evaluation  will  be  the  difficulty  of  the  exploration  drilling, 
which  is  divided  into  five  grades  which  also  mean  the  five  levels  of  investment: 

V  =  (I,  II,  III,  IV,  V). 

The  criteria  for  single-factor  evaluation  are  shown  in  Table  1.  For  better 
understanding  and  simpler  calculation,  in  our  evaluation  considering  the  single 
factor  of  drillability,  we  use  only  the  index  of  mechanical  drilling  speed. 

For  the  same  reason,  in  stability,  we  use  the  index  of  the  longitudinal  wave 
velocity  Vp  measured  from  the  core  sample;  in  difficulty  of  core  sampling,  we 
use  the  percentage  of  the  core  sampling  rate;  and  in  hole  inclination,  we  use 
the  result  of  boring  electrical  logging,  the  index  of  axial  declination  y 
calculated  from  the  bearing  and  the  apical  angle. 

Table  1.  Classification  of  Drillability  in  Exploration 


B  *  (2)^>^ 

-  I 

ii 

III 

IV 

V 

4®k 

( 

(4) 

(5) 

& 

rn 

wrz  ” 

R 

W  “ 

M.  K 

10>  1-3 

4-5 

6-7  i 

8—9 

10—12 

Ll)«tti§ 

<*/**> 

>2.-45 

2.45-1.15 

1.15-0.57 

0.57-0.25 

<0.25 

Qiik 

.OMP. 

c  rmwd 

V,  <=f 

19)  >5 

4-5 

3—4 

2-3 

<2 

(20) 

(22) 

mm 

T7¥V~ 

W— - 

. 

m  * 

(W) . — 

?R«S: 

(27)  ~ 

>90% 

75-90% 

60-75% 

45—60% 

<45% 

(28) 

(29 ) 

(30) 

*  m 

~wr~- 

(34)  ■' 

n 

<K/g*> 

5)  <1.5 

2±0 . 5 

3±0.5 

4±0.5 

>4.5 

Key 

11. 

Mechanical  drilling  speed 

1. 

Class 

12. 

Stability 

2. 

Factor 

13. 

Extent 

3. 

Drillability 

14. 

Long-term  stable 

4. 

Description  of  extent 

15. 

Stable 

5. 

Soft 

16. 

Medium  stable 

6. 

Mid-hard  (1) 

17. 

Unstable 

7. 

Mid-hard  (2) 

18. 

Seriously  unstable  (serious  fall- in) 

8. 

Hard 

19. 

Longitudinal  wave  velocity  of  rock 

9. 

Very  hard 

sample  Vp<km/sec> 

10. 

Twelve  rock  grades 

20. 

Drill  sampling 
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Table  1  continued: 


21. 

Extent 

29. 

Extent 

22. 

Very  easy 

30. 

Very  unlikely  deviation 

23. 

Easy 

31. 

Unlikely 

24. 

Medium 

32. 

Medium 

25. 

Hard 

33. 

Likely  deviation 

26. 

Very  hard 

34. 

Very  likely  to  deviate 

27. 

28. 

Core  sampling  rate 

Hole  deviation 

35. 

Axial  inclination  <degrees  per 
100  meters> 

In  any  specific  drilling  area,  though  information  can  be  obtained  from  the 
production  record  of  one  or  several  holes  drilled  and  the  results  of  rock 
specimen  tests,  rock  formations  are  different  in  nature  in  the  exploration 
area,  and  there  is  no  way  to  classify  the  whole  exploration  area  according  to 
the  classifications  in  Table  1,  no  matter  which  single  factor  we  choose  to  use. 
For  example,  if  drillability  is  used,  in  that  area  some  rocks  belong  to  incom¬ 
petent  beds,  some  to  mid— competent  beds  and  some  to  competent  beds.  Therefore, 
drillability  in  a  certain  area  can  only  be  a  fuzzy  idea  and  should  be  expressed 
with  a  fuzzy  aggregate: 

A  ■■  «<a/I  +  #<o/n  +  ttjlll  +  uajTV  +  m„5/V » 

0 

Here  the  Uai  is  the  extent  of  subordination  of  that  element  in  the  area  of 
Class  i  to  the  subaggregate  of  "area  drillability."  The  value  of  Uai  is  de¬ 
rived  as  follows:  among  m  times  the  calculations  done  (or  of  rock  samples 
tested  or  the  segments  of  holes  measured  or  other  tested  values) ,  there  are  n 
times  whose  mechanical  drilling  speed  (or  some  other  tested  value)  belongs  to 
Class  i  drilling  area.  The  total  length  of  the  former  is  Lm  and  the  total 
length  of  the  latter  is  Ln.  The  rate  of  subordination  should  be  their  ratio: 

Mi !<i>/ I'm. 

5 

Obviously,  and  the  fuzzy  vector  is  inductive.  Thus  we  get  the  four 

subaggregates  of  the  single-factor  evaluation: 

A  “  t<a2>  wb3»  W*j)» 

=  (k* is  Uh ««)> 

C  “=  Cucl>  uel>  uc3> 

Q  =  («i!>  «<J2>  Uiit  Ui 4,  Uds)t 

and  from  here  we  get  the  fuzzy  relational  matrix.  From  U  to  V 
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!ut t>  •••»  Ug,\ 

/rn>  i"a>.  ***»  fi5\ 

Kil)  ««,  *.*•»  Ub5 

rn>  *  •  * ,  fijl 

«eli  tle 2,  •••,,  Ucl 

i- . . . . 

»*>  ‘  UdJ 

r^,  •  •  rj 

This  is  the  matrix  derived  from  the  single  factor  evaluation.  The  element  rij 
represents  the  feasibility  in  the  class  j  prospecting  area  with  the  factor  i 
in  consideration. 

Again,  based  on  statistics,  if  the  percentage  of  the  cost  of  drilling,  wall 
protection,  leak  repair,  core  collection,  inclination  prevention  and  correction, 
etc.  in  the  total  investment  can  be  determined  by  the  operational  conditions 
and  technological  requirements ,  then  we  get  the  order  of  subordination 

U  “  (mi*  «a>  «39  m±). 

Finally,  through  combined  calculation  of  the  two  fuzzy  matrices,  we  get  the 
result  of  the  comprehensive  evaluation: 

Z-M'S- 

Now,  we  use  the  statistics  of  the  exploration  drilling  in  a  certain  area  and 
the  data  gathered  on  the  spot  as  our  bases  and  apply  the  above  method  to  make 
an  investment  calculation.  The  data  from  a  certain  drilling  are  listed  as 
Table  2. 

From  the  data  in  that  table  we  can  get  the  four  fuzzy  subaggregates  of  the 
single-factor  evaluation. 

For  drillability,  we  calculated  a  total  of  338.50  meters.  Included  in  it, 
when  drilling  speed  is  greater  than  2.45  m/hr,  progress  is  0  meters;  at  speeds 
of  2.45-1.15  m/hr,  it  is  2.15  meters;  at  1.15-0.57  m/hr  it  is  278.32  meters; 
at  0.57-0.25  m/hr,  it  is  58.03  meters;  and  when  the  speed  is  less  0.25  m/hr, 
the  length  is  also  0  meters.  Then  we  get 

k*  —  *“  0, 

u*  -  2.15/338.50  —  0.01, 

«„  —  278.32/338.50  —  0.82,  ■ 

-  58.03/338.50  -  0.17, 

A  —  (0,  0.01,  0.82,  0.17,  0). 

Similarly  we  can  calculate  according  to  the  stability:  0  meters  belongs  to 
Class  I,  210.83  meters  to  Class  II,  125.52  meters  to  Class  III,  2.15  meters 
to  Class  IV  and  0  meters  to  Class  V,  that  is, 
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Table  2.  Data  From  XX  Prospecting  Area  Well  ZK  7-22 


2.  Hole  Depth  £  %  M 

<M>  Ofc> 

3 .  Segment  . - - 

Length  M  1  86.65 

4.  Name  of  Rock  - - - - 

5.  Siliceous  2  123.41 

limestone _ _ 

6.  Siliceous  3  127. 00 

mudstone  _ _ 

7 •  Siderite  ^  j^9 

limestone  _ ' 

8.  Siliceous  ,  „ 

mudstone 

9.  Siliceous 

limestone  6  161*40 

10.  Siltstone  ' 

11.  Medium  fine-  7  216.71 

grain  stone  - - 

bind  8  218.70 

12.  Siltstone  — - - 

13.  Pelitic  '9  249.50 

siltstone  — - - 

14.  Silstone  10  278.14 

15.  Medium  fine- _ 

grain  stone  u  295.46 

bind _ _ 

16.  Silstone  12  298.70 

17 .  Sandy _ _ 

mudstone  13  299.8I 

18.  Silstone  _ _ _ 

19'  F1"®  ,  H  3«.« 

sandstone 

20.  Siltstone  T  ’  " 

21.  Sandy  15  339*3' 

mudstone  " 

22.  Fine  16  359-5 

sandstone  - 

23.  Mechanical  17  373.8' 

Drilling - - — 

speed  18  376.6 

(m/hr)  — — - - - 

24.  Longitudinal 

Wave  Velocity  <m/ sec> 

25.  Stratified  Drilling  Rate 

26.  Depth 

27.  Sampling  Rate 

28.  Record  of  Inclination 


.  (4) 

grok 

<*> 

<*/*> 

W 

83.45 

0.95 

4618 

86.65 

123.41 

31% 

37% 

36.76 

1.14 

4691 

127.00 
.  144.21 

56% 

13% 

3.59 

n&ani* 

0.77 

4874 

149.15 

149.95 

22% 

0% 

22.15  ■ 

0.96 

3835 

150.70 

151.30 

0% 

0% 

2.15 

1.31 

2666 

161.40 

216.71 

65% 

22% 

.  10.10 

0.97 

3789 

218.70 

219.00 

5% 

100% 

55.31 

0.47 

3676 

222.53 

223.73 

40% 

100% 

0,84 

mUdbg  i.oi 


"  0.57 


^  0.88- 


277.34 

278.14 

278.64 

295.46 


Depth 

Sampling  Rate 
Record  of  Inclination 


29.  M 

30.  Optical 

31.  Bearing 


(28) 

M413K 


30)50  *(29) 
-  1°20* 

100  * 

9  =  3°20' 
a  =  282° 

130 

9  =  4°0' 
a  =  290° 

150  # 

6°0' 
a  =  277° 

200  % 

9  =  5°40' 
a  =  285° 

250  * 

9=*  6°  10' 
a  =  280° 
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299.80 

1.10 

0.82 

1 

3623 

298.70 

299.80 

98% 

100% 

300  * 

9  =  8°0' 

14 

305.80 

6.00, 

(iH) 

»5J>g 

1.00 

4061 

305.80 

306.30 

83% 

0 

a  =  295° 

'  15 

339.36 

12.51 

0.66 

4060 

317.10 

321.70 

93% 

98% 

350  # 
6  =  9°30' 

16 

359.50 

20.14 

UUj 

0.90 

4253 

325.65 
,  326.85 

100% 

75% 

a  =  292° 

17 

373.89 

14.39 

1.08 

.3580 

339.36 

359.50 

48% 

80% 

360  * 

9  =  9°30' 

18 

376.61 

2.72 

0.37 

4522 

373.89 

376.61 

83% 

100% 

a  =  290° 
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«il  “  #iS*“  0> 

uu  T  210.83/338.50  —  0.62, 

125.52/338.50  —  0.37,  , 

««  -  2.15/338.50  -  0.01* 

§  -  (0,  0.62,  0.37,  0.01,  0). 

For  core  recovery,  373.61  meters  are  calculated.  Among  them,  66.13  meters 
with  a  core  recovery  rate  greater  than  90  percent  belongs  to  Class  I,  50.28 
meters  with  a  rate  of  90-75  percent  belongs  to  Class  IV,  17.74  meters  with  a 
rate  of  75-60  percent  belongs  to  Class  III,  32.92  meters  with  a  recovery  rate 
of  60-45  percent  belongs  to  Class  IV  and  206.54  meters  with  a  core  rate  of 


less  than  40  percent  belongs  to  Class  V. 

Thus, 

-  66.13/373.61 

—  '0.18,' 

—  50.28/373.61 

-  0.13, 

wrt  -  17.74/373.61 

—  0.05, 

—  32.92/373.61 

—  0.09, 

ua  —  206.54/373.61  —  0.55, 

Q  —  (0.18,  0.13,  0.05,  0.09,  0.55). 

The  hole  deviation,  first  we  calculate  the  axial  inclination  according  to  the 
apical  angle  and  the  bearing  used  for  the  electrical  logging,  as  shown  in 
Table  3.  The  value  of  the  axial  inclination  is  calculated  with  the  cosine 
formula : 

cost  —  cos  Si  cos  &2  +  sm  9  j  sin  cos  Aa. 

From  Table  3  and  by  calculation,  we  know  that  when  a  hole  deviation  of  360 
meters  has  been  measured,  we  get 

«*  —  50/360  —  0.14, 
uti  —  50/360  —  0.14, 

«*  —  140/360  —  0.39, 

uh  —  0, 

#*  —  120/360  —  0.33, 

Q  —  (0.14,  0.14,  0.39,  0,  0.33). 

Now  we  get 

7  0  ,  0.01,  0.82,  0.17,  0  V 

Rmt  0.  ,  0.62,  0.37,  0.01,  0  | 

0.18,  0.13,  0.05,  0.09,  0.55  J 

*0.14,  0.14,  0.39,  0  ,  0.33  / 
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Table  3.  Results  of  Hole  Deviation  Calculation  of  Well  ZK7-22 


U  i)  LL L-, - , - , - , - i-IP-T-m-  r 


$«(*:> 

0, 

-  9t 

a. 

r'  a 

rS/S* 

»  % 

0-50 

50 

0. 

1°20' 

0 

137“ 

1.33° 

2.66° 

III 

50-100 

50 

1620' 

3°20' 

■  137o 

282° 

4.49° 

8.98° 

v 

100-130 

30 

3°20' 

4°0' 

282° 

290° 

0.84° 

2.80° 

ra 

130-150 

20 

4°0' 

6°0' 

290° 

277° 

2.29° 

11.45° 

v 

,  150—200 

50 

6°o; 

5°40' 

277° 

285° 

0.88° 

1.76° 

u 

200-250 

50’ 

5°40' 

.  6°  10' 

285° 

280° 

0.72° 

1.44° 

l 

250—300 

50 

6°10' 

8°0' 

280° 

295° 

2.59° 

5.18° 

v 

300-350 

50 

8°0' 

9°30' 

295° 

292° 

1.57° 

3.14° 

in 

350-360 

10 

9°30' 

9°30' 

292° 

290° 

0.33° 

3.30° 

in 

Key: 

1.  Depth  (M)  4.  Degree/ 100M 

2.  Length  (M)  5.  Classification 

3.  Degree 

According  to  the  prospecting  class  and  the  actual  situation  of  the  example,  we 
presume  the  subordination  of  each  factor  as  follows: 

A/  —  (0.50,  0.15,  0.10,  0.25). 

Finally,  we  get  the  result  of  the  comprehensive  evaluation: 


M°B 

(0.50, 0.15, 0.10, 0.25)o 

*  • 

(0.14,  0.15,  0.50,  0.17,  0.25). 


/  0  , 

a.oi. 

0.82, 

0.17, 

0 

\o  , 

0.62, 

0.37, 

0.01, 

0 

10.18, 

0.13, 

0.05, 

0.09, 

0.55 

'0.14, 

0.14, 

0.39, 

0  , 

0.33 

Factors  of  the  comprehensive  evaluation  must  also  be  summed  up  as 

0,14  +  0.15  +  0.50  +  0.17  +  0.25  -  1.21. 

so  that  the  final  result  is 

/  0.14  0.15  0.50  0.17  0.25  \ 

\i;21*  1.21*  1.21*  1.21*  1.21/ 

—  (0.12,0.12,0.41,0.14,0.21). 
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The  above  result  indicates  that  this  drilling  area  should  belong  to  Class  IV 
investment  for  prospecting  area. 
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[Article  by  Yang  Xiaozhen  [2799  1420  7201],  Hu  Shiru  [7579  0013  1172],  ^ 

Lu  Yafei  [0712  0068  7236],  Zhu  Shannong  [2612  0810  6593]  and  Li  Xiangkui 
[2621  4161  7608]  of  the  Chemistry  Institute  of  the  Chinese  Academy  of 
Sciences:  "A  Study  of  the  Structure  of  Nylon  1010";  paper  received  on 

22  April  1983] 

[Text]  Abstract:  Nylon  1010  specimens  prepared  from 
both  vitreous  and  melted  crystalline  states  were  studied 
by  x-ray  diffraction,  infrared  spectroscopy,  differential 
scanning  calorimetry,  polarized  optical  microscopy,  viso- 
elasticity  and  density.  It  was  found  that  its  structure 
is  similar  to  that  of  Nylon  66.  However,  the  lattice  is 
greatly  distorted  when  crystallized  from  the  melt.  Related 
problems  were  discussed  and  a  hydrogen  bond  plane  distor¬ 
tion  model  was  introduced. 

Nylon  1010  is  a  special  engineering  plastic  in  China.  Its  chain  structure  and 
polymeric  structure  have  been  studied. 1,2  In  order  to  understand  the  struc¬ 
ture  changes  of  Nylon  1010  crystallized  from  different  states,  we  investigated 
the  polymeric  structures  and  crystals  of  specimens  prepared  from  three  differ¬ 
ent  scanning  calorimetry,  polarized  optical  microscopy,  viscoelasticity  and 
density.  The  first  type  is  quenched  from  a  melt.  The  second  type  is  crystal¬ 
lized  isothermally  from  a  vitreous  state.  The  third  type  is  crystallized  from 
a  melt.  We  found  that,  in  addition  to  its  crystal  structure  being  similar  to 
that  of  Nylon  66,  the  lattice  of  the  third  type  crystal  is  severely  distorted, 
resulting  in  its  unique  x-ray  diffraction  pattern  as  well  as  changes  of  the 
infrared  band,  melting  peak  and  kinetic  behavior.  We  believe  that  it  is 
reasonable  to  explain  this  phenomenon  by  the  distortion  of  the  hydrogen  bond 
plane. 


^Current  address  is  Chemistry  Department,  Zhengzhou  University. 
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Experiments  and  Results 
Sample  Preparation 

Nylon  1010  resin  used  was  produced  by  Zhengzhou  University  (castor  oil  line) . 
It  was  dissolved  in  methanol  and  cast  into  a  film  on  a  glass  plate.  Then,  it 
was  placed  in  a  fixture  and  immersed  in  a  230°C  salt  bath  for  3-8  minutes. 

The  three  types  of  specimens  were  prepared  as  follows:  type  I  was  quenched  in 
dry  ice  and  heptane,  type  II  was  quenched  in  dry  ice  and  heptane,  crystallized 
isothermally  at  150-180°C  for  several  minutes  to  1  day  (under  nitrogen)  ,  and 
then  air  cooled  at  room  teperature,  and  type  III  was  crystallized  isothermally 
at  170-180°C  under  nitrogen  for  several  minutes  to  2  hours  and  then  air  cooled 
at  room  temperature. 

Thickness  Determination 

A  micrometer  was  used  to  measure  six  points  and  report  the  average.  Sample 
thickness  is  80-120  y. 

Density  Determination 

A  density  gradient  tube  with  carbon  tetrachloride  and  heptane  was  used.  The 
density  range  is  1.008-1.105  g/cm3  and  the  resolution  is  0.0003  (g/cm3)/cm. 
Specimens  were  measured  at  20 °C. 

Polarized  Optical  Microscopy 

A  Model  XPT-6  optical  microscope  manufactured  by  Jiangnan  Optical  Instrument 
Plant  with  a  magnifying  power  of  10  x  25  was  used.  Figures  1(A)  and  (B)  show 
type  I  and  II  specimens.  We  did  not  see  any  spherocrystal.  However,  in  type 
III  specimen  (C)  the  presence  of  spherocrystals  is  apparent.  Based  on  color, 
they  were  identified  as  positive  spherocrystals.  The  white  line  in  (A)  indi¬ 
cates  some  kind  of  molecular  order  which  suggests  that  Nylon  1010  can  easily 
crystallize.  Even  though  it  is  quenched  immediately  after  melting,  it  is  not 
possible  to  obtain  totally  amorphous  specimens. 


Figure  1.  The  Three  Specimens  Under  a  Polarized  Optical  Microscope 
(A)  Type  I;  (B)  Type  II;  (C)  Type  III 
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X-ray  Diffraction 

A  Scientific  Electrical  Equipment  Model  2028  x-ray  dif fratometer  with  a  Cu 
target,  a  Ni  filter  and  a  1°  incident  slit  was  used.  The  sample  film  was 
adhered  to  a  glass  plate  using  Canadian  glue.  Typical  x-ray  diffraction 
patterns  of  these  three  types  of  specimens  are  shown  in  Figure  2.  The  crystal 
structure  of  Nylon  1010  exhibits  two  distinct  diffraction  peaks  at  20°  and  24°, 
similar  to  those  of  Nylon  66. 3  The  analog  diffractionoplanes  are  100  and 
010  +  110,  respectively.  The  spacings  are  dioo  =  4.44A  and  doioq.no  ~  3.7lA. 
There  is  an  unnoticeable  peak  near  8.5°.  The  type  I  peak  is  diffused,  indi¬ 
cating  that  the  specimen  is  almost  amorphous.  The  type  II  peak  is  similar  to 
that  of  Nylon  66 . 3  The  100  plane  diffraction  peak  is  very  high  and  sharp  for 
type  III  and  its  010  +  110  peak  is  small  and  diffused. 


Figure  2.  X-ray  Diffraction  Patterns  for  Three  Types  of  Specimens 


Infrared  Spectroscopy 

A  Perkin-Elmer  Model  180  infrared  spectrophotometer  was  used  to  record  the 
spectra.  The  resolution  is  2  cm-1.  The  bands  at  940  and  1007  cm-1  are 
related  to  the  vibration  of  the  nylon  skeleton.1’2  Figure  3  was  prepared  by 
normalizing  absorbance  and  density  as  well  as  by  taking  sample  thickness  into 
account.  Type  I  and  II  specimens  seem  to  vary  linearly  with  density,  while 
the  density  of  type  III  specimen  seems  to  be  independent  of  absorptivity. 

Differential  Scanning  Calorimetry 

A  Perkin-Elmer  Model  DSC-II  calorimeter  was  used.  The  temperature  rising  rate 
was  20°C,  the  range  was  5  mCal/sec  and  the  nitrogen  flow  rate  was  45  ml/min. 
Results  obtained  in  the  400-500°K  range  are  as  follows: 

Figure  4  shows  a  high  melting  peak  at  near  471°K  for  type  I  and  II  samples, 
while  the  main  melting  peak  for  type  III  specimens  appears  at  463°K.  The 
difference  is  apparent.  This  is  because  in  the  temperature  rising  process  we 
have  to  experience  the  optimal  crystallization  temperature,  which  is  estimated 
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to  be  165°C,  despite  that  the  temperature  rising  rate  is  high  (20°C).  Thus, 
crystallization  occurred  in  DSC  measurements.  Of  course,  this  happened  to  all 
samples.  The  melting  peak  of  type  I  is  a  few  degrees  lower  and  the  peak  is 
wider  than  that  of  type  II,  which  is  more  crystalline.  This  effect  may  be 
caused  by  the  superposition  of  the  crystal  and  the  melt,  as  well  as  by  the 
superposition  of  the  relaxation  of  the  crystallization  process  and  the  melting 
process.  The  small  peak  before  the  melting  peak  of  type  II  specimens  is 
basically  related  to  the  temperature  of  the  thermal  treatment.  It  is  approxi¬ 
mately  10°C  higher  than  the  thermal  treatment  temperature. 

Viscoelasticity 

A  Model  DDV-II  Dongyang  instrument  was  used.  The  frequency  was  110  Hz,  the 
temperature  rising  rate  was  2°C/min  and  the  temperature  range  was  from  liquid 
nitrogen  temperature  to  100°C.  Figure  5  shows  typical  viscoelasticity  patterns 
of  these  three  types  of  samples.  Within  the  temperature,  there  are  three 
transition  peaks.  At  approximately  70°C,  it  is  the  vitreous  or  a  transition, 
-48°C  corresponds  to  a  6  transition  and  -118°C  is  the  y  transition.  The  two 
latter  transitions  correspond  to  the  movement  of  the  amid  group  and  the  carbon 
chain,  respectively.4  Among  all  specimens  tested,  the  a  transition  peaks  of 
type  III  samples  are  always  less  than  those  of  types  I  and  II. 


Figure  5,  Viscoelasticity  Patterns  of  Three  Types  of  Samples 
®  =  I;  o  =  II;  A  =  III 


Discussion 

From  the  above  results  one  can  see  that  there  are  obvious  differences  between 
types  II  and  III  in  terms  of  crystal  structure,  x-ray  diffraction  pattern, 
relation  between  skeletal  vibration  and  density,  enthalpy  variation  and  kine¬ 
tic  characteristics.  We  believe  that  these  differences  indicate  that  they 
are  different. 
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Type  II  is  similar  to  other  crystalline  polymers  crystallized  from  the  vitreous 
state.  The  resulting  spherocrystal  is  smaller.  Normally  a  polymer  forms  a 
larger  spherocrystal  from  the  melting  state.  The  lattice  is  orderly  over  a 
long  span  in  three  dimensions.  Although  type  III  specimens  crystallized  from 
the  melt  shows  large  spherocrystals ,  however,  the  height  and  half  width  of  the 
two  diffraction  peaks  are  quite  different.  This  indicates  that  it  only  grows 
orderly  in  one  dimension  and  there  is  no  apparent  change  in  other  directions. 

It  may  even  decrease  in  the  direction  of  the  010  +  110  plane,  which  is  the 
hydrogen  bond  plane.  It  is  a  mesh  created  by  linking  hydrogen  bonds  among  the 
stretched  molecular  chain.  The  direction  of  the  hydrogen  bond  is  the  direction 
of  the  100  plane.  Therefore,  the  degree  of  orderliness  in  large  type  III 
crystals  is  nonuniform. 

The  loss  peak  tg6max  obtained  from  the  viscoelasticity  due  to  Tg  transition 
is  compared  to  the  absorbance  at  1007  and  940  cm  1  on  the  infrared  spectra 
(see  Figure  6).  The  intensity  of  the  loss  peak  decreases  with  increasing 
intensities  of  these  two  absorption  bands  for  type  II  specimens.  However, 
type  III  samples  are  different.  Although  the  infrared  spectrum  peak  height 
varies  a  lot,  yet  the  loss  peak  remains  unchanged.  Their  loss  peaks  are  less 
than  those  of  type  II  samples,  indicating  that  there  is  less  noncrystalline 
area  in  type  III  samples.  At  the  two  upper  points  for  type  II  specimens,  the 
densities  are  higher  than  those  of  type  III  samples.  This  abnormal  effect  is 
attributed  to  the  distortion  of  the  crystal  structure  of  type  III  specimens. 

The  DSC  melting  temperature  of  type  III,  which  characterizes  its  crystal 
structure,  is  approximately  10°C  lower  than  that  of  type  II,  reflecting  a 
distorted  structure  as  well.  The  structure  distortion  in  the  geometric  space 
leads  to  the  lowering  of  internal  polymerization  energy.  The  density  decrease 
in  the  crystalline  region  is  also  understandable. 


Figure  6.  Infrared  Absorbance  vs.  tg6max  for  Two  Types  of 
Samples 

(Absorbance  x  0.1) 
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Hosemann  investigated  the  distortion  of  crystalline  polymers  in  depth.  It 
usually  belongs  to  type  II  defects.  The  lattice  defect  we  detected  is  more 
severe  than  that  he  described.  Moreover,  the  distortion  in  various  directions 
is  different.  The  cause  and  explanation  of  the  formation  of  this  structure  and 
a  description  of  the  picture  of  the  distorted  structure  are  as  follows. 

Because  the  activation  energy  of  molecular  chain  movement  is  lower  for  type  III 
crystals,  the  chain  segment  moves  faster  and  easier.  In  this  case,  the  free 
energy  for  nucleus  formation  is  relatively  high.  The  number  of  nuclei  formed 
is  less,  causing  the  molecular  chain  to  stretch  forming  a  larger  hydrogen 
bond  cross-linked  plane.  This  naturally  makes  the  diffraction  peak  of  the 
100  crystal  plane  high  and  narrow.  The  hydrogen  bond  mesh  plane  is  orderly 
aligned  under  the  influence  of  van  der  Waal  force  which  is  much  weaker  than 
hydrogen  bonding  force.  It  is  comprehensible  that  this  effect  causes  a  wide 
distribution  of  crystal  plane  spacing.  However,  experimental  results  showed 
much  more.  The  fact  that  the  absorbance  of  the  infrared  band  corresponding 
to  the  skeletal  vibration  of  type  III  specimens  is  independent  of  density 
indicates  that  Che  molecular  chains  in  the  crystalline  region  are  not  all  in 
a  planar  saw-tooth  structure.  This  makes  it  impossible  to  form  a  real  hydrogen 
bond  plane.  At  around  180°C,  it  already  reaches  the  melting  temperature  of 
hydrocarbon  chain  macromolecules.  In  the  hydrogen  bond  plane  of  Nylon  1010 
with  a  long  hydrocarbon  chain,  van  der  Waal  force  is  not  enough  to  control 
the  hydrocarbon  movement.  However,  the  hydrogen  bond  is  still  very  strong 
which  leads  to  the  wide  fluctuation  of  spacing  of  hydrogen  bond  plane  and  the 
significant  distortion  of  the  hydrogen  bond  plane.  Consequently,  the  effect 
that  the  diffraction  peak  of  the  010  +  110  crystal  plane  of  the  larger  type  III 
spherocrystal  is  not  increased  can  be  explained.  This  effect  is  not  apparent 
for  type  II  specimens  because  a  large  number  of  nuclei  is  formed  after  quench¬ 
ing.  When  it  is  isothermally  recrystallized,  even  at  180°C,  it  is  not  possi¬ 
ble  to  make  all  the  nuclei  disappear.  Consequently,  it  is  difficult  to  form 
a  larger  hydrogen  bond  mesh  surface.  Because  the  hydrogen  bond  surface  is 
smaller,  the  scale  of  orderliness  is  small  in  all  directions.  The  effect  of 
hydrogen  bond  surface  distortion  is  reduced  and  the  drastic  difference  of  dis¬ 
tortion  along  various  directions  is  also  reduced.  Of  course,  it  belongs  to  the 
second  kind  of  lattice  distortion  specified  by  Hosemann. 

Meaningful  information  can  also  be  obtained  from  the  8  and  y  transition  peaks 
measured  kinetically.  8  peak  corresponds  to  the  motion  of  the  nonhydrogen 
bonded  amid  group.  In  Figure  5,  the  8  peak  for  type  III  specimens  is  the 
smallest,  indicating  that  there  are  more  hydrogen  bonds  in  the  crystalline 
region  which  is  in  agreement  with  its  high  degree  of  crystallinity,  y  peak 
corresponds  to  the  carbon  chain  movement.  The  y  peak  for  type  III  specimen  is 
not  the  smallest,  indicating  that  the  number  Of  orderly  carbon  chains  in  the 
lowest  ground  state  is  not  proportional  to  the  large  crystalline  region.  In 
other  words,  there  are  irregular  carbon  chains  in  the  crystalline  region. 

This  kinetic  characteristic  also  supports  the  above  hydrogen  bond  plane  dis¬ 
tortion  model. 
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DIFFERENTIAL  ELASTIC  SCATTERING  CROSS  SECTIONS  OF  14.2  MeV  NEUTRONS  FROM 
NIOBIUM 

Beijing  YUANZIHE  WULI  [CHINESE  JOURNAL  OF  NUCLEAR  PHYSICS]  in  Chinese  Vol  7  No  2, 
May  85  pp  106-111 

[Article  by  Li  Jingde  [2621  2529  1795] ,  Xie  Daquan  [6200  1129  0356],  Li  Yexiang 
[2621  2814  4382],  Zhang  Di  [1728  2769],  Wang  Shiming  [3769  0013  2494]  and  Cao 
Jianhua  [2580  1696  5478]] 

[Text]  I.  Introduction 

Elastic  scattering  is  the  primary  physical  process  in  the  interaction  between 
fast  neutrons  and  matter.  Accurate  data  of  differential  elastic  scattering 
cross  sections  not  only  is  valuable  in  the  design  of  nuclear  reactors  and  in 
other  nuclear  projects,  but  also  provides  a  reliable  basis  for  adjusting 
optical  model  parameters,  verifying  theoretical  models,  and  predicting  measured 
data  in  energy  deficient  zones. 

The  superior  physical  properties  of  niobium  makes  it  an  attractive  construction 
material  for  reactors.  In  fusion  reactors,  there  is  growing  interest  in  using 
niobium  as  the  cover  layer  for  the  fusion  plasma;  consequently,  niobium  neutron 
data  is  becoming  increasingly  more  important.  In  the  literature,  there  have 
been  two  reported  measurements  of  differential  elastic  scattering  cross  sections 
of  niobium  fast  neutrons  in  the  neighborhood  of  14  MeV:  the  measurement  by 
Western,  et  al.,  at  20  different  angles  between  20°  and  140°, 1  and  the  measure¬ 
ment  by  Kammerdiener  at  18  different  angles  between  25°  and  145°. 2  The  two  sets 
of  data  differ  substantially  for  angles  less  than  45°,  and  there  is  no  experi¬ 
mental  data  below  20°.  There  are  also  significant  differences  between  the 
experimental  values ,  the  theoretical  values ,  and  recommended  values  of  the  total 
elastic  scattering  cross  sections.  The  integrated  elastic  scattering  cross 
section  given  by  Western,  et  al.,  is  1.77b,  which  obviously  does  not  agree  with 
the  value  2.12b  obtained  by  Angeli,  et  al.,3  using  mass  system  theory;  the  cal¬ 
culated  cross  sections  based  on  the  optical  model  given  by  Smith,  et  al. ,  fall 
between  2.3b  and  3b;  the  recommended  values  are  between  1.77b  and  2.3b;  thus, 
these  results  differ  by  25  percent  or  greater.4-7  The  large  discrepancies  indi¬ 
cate  that  there  is  clearly  a  need  to  improve  the  quality  of  the  experimental 
data.  In  the  evaluation  report,  Smith,  et  al. ,  pointed  out  that  by  improving 
the  measurement  of  14  MeV  elastic  scattering  cross  sections,  particularly  in 
the  range  0°-45°,  the  quality  of  evaluation  can  be  significantly  improved. 
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In  view  of  the  importance  of  the  data  in  the  14  MeV  energy  zone  for  applied  and 
theoretical  research,  and  in  order  to  clarify  the  differences  of  experimental 
results  and  to  fill  the  void  of  small  angle  data,  we  have  made  laboratory  mea¬ 
surements  of  the  differential  elastic  scattering  cross  sections  of  14.2  MeV 
neutrons  in  niobium  at  28  different  angles  between  6°  and  150°;  these  angular 
data  are  fitted  with  Legendre  polynomials  to  obtain  the  integrated  elastic 
scattering  cross  sections. 

II.  Experimental  Setup  and  Measurement  Technique 

The  experiment  was  conducted  on  the  400  KV  high-voltage  multiplier  at  the 
Sichuan  University  Nuclear  Science  Research  Institute  using  the  associated 
particle  time  of  flight  method;  the  experimental  setup  is  shown  in  Figure  1. 

The  200  KeV  deuterium  beam  generated  by  the  high-voltage  multiplier  bombards 
the  tritium-titanium  target,  and  produces  neutrons  through  T(d,n)  He  reaction. 
The  "associated"  a  particles  are  detected  using  a  50  ym  thin-film  plastic 
scintillator  which  is  coupled  to  the  56  AVP  type  photoelectric  multiplier. 

This  detector  is  positioned  at  90°  relative  to  the  deuterium  beam;  the  corre¬ 
sponding  neutron  beam  which  is  constrained  by  the  rectangular  diaphragm  is 
located  82?  relative  to  the  deuterium  beam;  the  average  neutron  energy  is  14. 
MeV. 


In  this  experiment,  two  different  neutron  detectors  were  used.  The  first  con¬ 
sists  of  a  ST451  liquid  scintillator  5  cm  in  diameter  and  5  cm  thick  which  is 
coupled  to  the  56  AVP  photoelectric  multiplier;  it  is  used  to  measure  the  posi¬ 
tion  and  shape  of  the  associated  neutron  beam.  The  second  consists  of  a  ST451 
liquid  scintillator  10.4  cm  in  diameter  and  5  cm  thick  which  is  coupled  to  the 
XP2041  photoelectric  multiplier;  it  is  used  to  measure  the  differential  elastic 
scattering  cross  section.  In  order  to  reduce  background  noise,  the  neutron 
detector  is  located  inside  the  shield  structure.  The  shield  structure  is  1.3  m 
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in  diameter,  2  m  long,  and  filled  with  paraffin  and  lithium  carbonate  in  equal 
proportion  plus  an  appropriate  amount  of  ethylene  trichloride;  the  inner  layer 
is  a  5  cm  thick  lead  ring.  To  minimize  the  interference  due  to  hole-wall 

scattering,  the  alignment  hole  is  made  into  a  shape  of  a  double  frustum  cone. 

The  shield  weighs  3.5  x  10 3  kg,  and  is  carried  on  an  air-cushioned  cart;  it 
can  be  rotated  freely  about  the  axis  of  the  test  sample. 

The  position  and  shape  of  the  associated  neutron  beam  are  carefully  measured 
in  both  the  horizontal  and  vertical  directions.  The  associated  beam  has  a 
horizontal  opening  of  ±1.6°  (see  Figure  2)  and  a  vertical  opening  of  ±2.2°. 

With  the  aid  of  a  bidirectional  laser  collimator,  the  scattering  test  sample 
can  be  placed  in  the  center  of  the  associated  beam  so  it  is  surrounded  by  the 

beam.  The  axis  of  symmetry  of  the  test  sample  is  perpendicular  to  the  beam 

axis  of  symmetry;  the  center  of  the  test  sample  is  71  cm  from  the  target  center, 
and  the  flight  distance  of  the  scattering  neutron  is  230  cm. 


e 

Figure  2.  Position  and  Shape  of  the  Associated  Beam  in  the  Horizontal 
Direction 

The  shadow  cone  absorber  is  made  of  stacked  iron  bricks;  it  is  longer  than 
40  cm,  and  its  shape  and  position  are  determined  by  the  geometry  of  the 
measuring  angles. 

The  electronic  circuit  diagram  is  shown  in  Figure  3.  A  fast  signal  with  very 
short  time  constant  is  generated  from  the  anodes  of  the  photoelectric  multi¬ 
pliers  of  the  a  detector  and  the  neutron  detector,  and  passed  through  a 
constant-ratio  timing  device  to  provide  timing  signals  for  the  a  particles 
and  neutrons.  The  neutron  timing  signal  with  lower  count  rate  is  used  as  the 
"start"  signal  for  the  timing  analyzer,  and  the  a  timing  signal  with  higher 
count  rate  is  used  as  the  "stop"  signal  after  an  appropriate  delay.  The 
distribution  of  the  neutron  flight  time  is  converted  by  the  timing  analyzer 
into  corresponding  pulse  amplitude  distribution,  and  then  processed  by  the 
1024-channel  pulse  amplitude  analyzer. 
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Semiconductor  detection 


From  the  13th  pole  of  the  a  detector  photoelectric  multiplier,  a  linear  signal 
is  extracted  and  passed  through  the  timing  channel  with  the  proper  threshold 
and  window  width  to  produce  the  required  ot  particles;  they  are  recorded  by  a 
calibrator  which  is  also  used  to  monitor  the  flux  of  the  neutron  beam.  From 
the  12th  pole  of  the  neutron  detector  photoelectric  multiplier,  a  linear  slow 
signal  is  extracted  and  filtered  through  the  single  timing  channel;  the 
threshold  value  of  the  spectrometer  is  selected  based  on  measurement  needs. 

A  pulse  shape  discriminator  is  used  to  separate  the  neutrons  from  y  particles 
and  to  eliminate  y  background.  The  output  from  the  ot  timing  channel,  the 
neutron  linear  timing  channel,  and  the  pulse  shape  discriminator  are  passed 
through  a  two-stage  adjustable-delay  double-coincidence  circuit  to  provide 
the  open  signal  for  the  multichannel  pulse  amplitude  analyzer  where  coinci¬ 
dence  spectral  measurements  are  made. 

The  neutron  energy  threshold  is  chosen  to  be  1.5  MeV,  and  the  corresponding 
time  resolution  is  1  ns.  The  a  count  rate  is  4  x  104  count/s.  The  performance 
and  operating  condition  of  the  spectrometer  are  tested  using  a  C  and 
1H  (polyethylene)  test  sample. 

The  scattering  sample  is  a  93Nb  metallic  cylinder  2  cm  in  diameter  and  4  cm 
in  height;  its  purity  is  better  than  99.5  percent.  Alternate  measurements 
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of  the  flight  time  spectrum  with  the  test  sample  and  the  background  spectrum 
without  the  test  sample  are  made  at  28  different  angles  between  6°  and  150°, 
where  the  scattering  angle  is  measured  from  the  axis  of  the  associated  neutron 
beam.  Two  or  three  repeated  measurements  are  made  at  most  angles  with  good 
results  of  repeatability.  All  measurements  are  normalized  with  respect  to  the 
a  monitor  count.  During  the  entire  measurement  process,  the  stability  of  the 
energy  threshold  of  the  spectrometer  is  being  constantly  monitored. 

The  neutron  count  at  0°  angle  without  the  test  sample  is  also  measured  at  the 
same  a  count  rate. 

III.  Data  Processing  and  Correction 


For  the  geometric  configuration  of  this  experiment,  the  differential  elastic 
scattering  cross  section  is  computed  according  to  the  following  formula: 


d o  M( 0)  s(0)  (  IL  \,  1 

dflw“M(0)  '  N(  0)  ‘  a(0)  \l+LJ  'NT 


(1) 


where  N(0)  is  the  net  elastic  scattering  neutron  count  at  angle  after  compen¬ 
sating  for  the  background;  N(0)  is  the  neutron  count  at  0°  without  the  test 
sample;  M(0)  and  M(0)  are,  respectively,  the  total  count  of  a  associated  par¬ 
ticles  corresponding  to  the  measurements  N(0)  and  N(0);  e(0)  is  the  detection 
efficiency  of  elastic  scattering  neutrons  at  angle  0;  e(0)  is  the  detection 
efficiency  of  primary  neutrons  at  0°;  i  is  the  distance  between  the  test 
sample  and  the  target  center;  L  is  the  distance  between  the  test  sample  and 
the  detector;  Nt  is  the  total  number  of  nuclei  of  the  test  sample.  Nt  can  be 
determined  by  weighing  the  test  sample  using  a  precision  scale.  Thus,  the 
problem  reduces  to  one  of  measuring  the  other  physical  quantities  mentioned 
above . 


In  this  experiment,  we  made  the  approximation  e(O)/e(0)  =  1  because  the 
energy  of  the  scattering  neutrons  is  very  close  to  that  of  the  incident 
neutrons;  the  resulting  error  in  the  differential  cross  section  is  less  than 
1  percent. 

The  large  amount  of  measurements  of  background  spectra  indicate  that  above 
20°,  the  background  is  randomly  distributed  with  respect  to  the  flight  time; 
there  is  no  evidence  of  time-correlated  structure.  The  accidental  coinci¬ 
dence  background  on  the  right  side  of  the  elastic  peak  with  test  sample  is 
generally  consistent  with  that  with  the  test  sample  removed.  Below  20°,  due 
to  geometric  constraints,  the  scattering  effect  between  the  target  lining 
and  air  cannot  be  totally  eliminated;  also,  the  angular  distribution  of  this 
type  of  scattering  is  front  loaded,  hence  the  corresponding  elastic 
peak  position  on  the  background  spectrum  has  a  weak  structure.  The  back¬ 
ground  attenuation  due  to  the  presence  of  test  sample  is  taken  into  consider¬ 
ation  in  data  processing. 

Correction  is  made  to  the  elastic  scattering  neutron  peak  based  on  the  shape 
of  the  primary  spectrum,  thus  the  contribution  of  the  nonelastic  scattering 
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neutrons  is  removed.  The  loss  in  count  caused  by  the  dead  time  in  the  elec¬ 
tronic  circuit  is  also  corrected.  In  addition,  Monte-Carlo  method  is  used 
to  correct  for  neutron  flux  attenuation,,  multiple  scattering,  and  limited 
angular  resolution  of  the  system. 

IV.  Results  and  Discussion 

The  calculated  and  corrected  differential  elastic  scattering  cross  section  is 
shown  in  Figure  4.  The  total  error  is  between  4.7  and  8.5  percent.  An  error 
analysis  is  presented  in  Table  1. 


Table  1.  Error  Analysis 


Error  source 

Error,  percent 

Count  statistics 

1.5-7. 2 

Nonelastic  correction 

<2 

Geometric  factor 

<1 

Correction  for  dead  time 

in  electronic  circuit 

<1 

Detector  efficiency 

<1 

Correction  for  flux  attenuation  and  multiple  scattering 

3.6 

Number  of  nuclei  in  test 
weight  measurement) 

sample  (including  error  in 

0.5 

Total 

4. 7-8. 5 

For  comparison.  Figure  4  also  shows  the  experimental  data  of  References  [1,2] 
and  the  theoretical  angular  distribution  calculated  using  the  optical  model 
given  in  Reference  [9]. 


Figure  4.  Differential  Elastic  Scattering  Cross  Section  of  14.2  MeV  Neutron 
From  9  3Nb 

x — Ref  1  14.7  MeV  A— Ref  2  14  MeV 

o — current  experiment  solid  curve  represents  the  theoretical  result 

based  on  optical  method 


The  measured  angle  distribution  data  from  this  experiment  are  supplemented  by 
computed  data  for  small  and  large  angles  using  the  optical  mode;  at  0°,  they 
satisfy  the  Wick  limit10: 


where  k  is  the  wave  number,  and  cr^  is  the  total  cross  section. 

By  approximating  the  above  data  using  optimum  Legendre  polynomials,  and  inte¬ 
grating  the  approximating  curve,  we  can  calculate  the  integral  elastic  scatter¬ 
ing  cross  section  to  be  2083  ±  76.6  mb.  The  expression  of  Legendre  polynomials 
used  to  represent  the  differential  elastic  cross  section  is 

■  P,W-  (3) 

where  0  is  the  arc  length  of  the  mass  system  (O<0<?r),  =  and  the 

integral  elastic  scattering  cross  section  is  a  (2?)  =  2a\ (£•  •  sinfldfl 

J  o  di3 

The  measured  integral  elastic  scattering  cross  section  values  from  this  experi¬ 
ment  are  in  good  agreement  with  the  value  of  2036.65  mb  computed  in  Reference 
[9]  using  optical  model  parameters  and  with  the  value  of  2120  mb  obtained  in 
Reference  [3]  using  mass  system  theory. 

In  order  to  obtain  reliable  data  on  integral  elastic  scattering  cross  sections, 
it  is  necessary  not  only  to  improve  the  measurement  of  differential  cross 
sections,  but  also  to  make  reasonable  assumptions  about  differential  cross 
sections  in  large  and  small  angular  regions  where  experimental  data  are 
lacking.  This  issue  requires  further  discussion. 

Due  to  constraints  in  experimental  conditions,  there  is  a  lack  of  experimental 
data  in  both  the  small-angle  region  (0°-20°)  and  the  large-angle  region  (160°- 
180°).  In  the  past,  manual  extrapolation  techniques,  extrapolation  using 
Legendre  polynomial  approximation,  or  analytical  methods  based  on  theoretical 
models  have  been  used  to  fill  these  gaps.  For  low  incident  neutron  energy, 
these  methods  yield  similar  results  because  variations  in  the  angular  distri¬ 
bution  are  generally  quite  small.  However,  as  the  energy  level  increases,  the 
angular  distribution  becomes  more  and  more  front  loaded  and  large 
differences  between  these  methods  become  apparent;  generally  speaking,  the 
higher  the  energy  level  and  the  heavier  the  nuclear  weight,  the  larger  the 
differences.  For  example,  our  calculations  show  that  in  the  integral  elas¬ 
tic  scattering  cross  section  of  Niobium  in  the  neighborhood  of  14  MeV,  more 
than  70  percent  of  the  contribution  comes  from  the  region  below  20°,  where 
no  experimental  data  existed  prior  to  this  measurement.  Western,  et  al., 
obtained  the  elastic  scattering  cross  section  of  1770  mb  by  extrapolating 
from  the  blackcore  model1;  Smith,  et  al. ,  applied  extrapolation  and  interpo¬ 
lation  processing  to  the  same  angular  distribution  data  to  obtain  the  value 
of  2050  mb4;  the  difference  between  the  two  is  280  mb. 
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We  believe  that  a  reasonable  approach  would  be  to  use  the  Wick  limit  at  0°,  and 
fill  the  data  gaps  at  small  and  large  angles  using  optical  model  calculations, 
then  apply  Legendre  polynomial  approximation  to  obtain  the  integral  elastic 
scattering  cross  sections. 

In  summary,  the  main  contribution  of  this  paper  is  to  provide  experimental 
data  in  the  small-angle  region,  and  to  improve  the  reliability  of  integral 
elastic  scattering  cross  sections  by  using  an  improved  method  of  calculation. 

Thanks  are  due  to  Comrade  Wan  Hongsheng  for  calculating  the  theoretical  opti¬ 
cal  model  curve,  and  to  members  of  the  high-voltage  multiplier  group  for  their 
support  in  carrying  out  this  experiment. 
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Fudan  University;  "Differential  Circuit  Model  of  Turbomolecular  Pump  Blade  and 
Calculation  of  Combination  Turbine  Compression  Ratio"] 

[Text]  Based  on  a  differential  circuit  model  for  turbomolecular  pump  blades, 
the  simple  relation  of  its  compression  ratio  versus  the  leak  gap  and  impeller 
spacing  was  obtained.  In  addition,  a  comprehensive  correction  method  for 
compression  ratio  was  introduced.  Calculated  results  agreed  with  experimental 
data  well.  Furthermore,  the  relation  between  compression  ratio  and  impeller 
impeller  spacing  is  not  as  close  as  predicted  by  the  molecular  drag  theory. 

The  compression  ratio  of  a  single  stage  turbomolecular  pump  has  already  been 
derived  from  the  Monte  Carlo  method[l]and  statistical  theory[2,3].  The  com¬ 
pression  ratio  of  a  combination  turbomolecular  pump  is  generally  believed  to 
be  simply  the  product  of  the  compression  ratios  of  various  stages.  However, 
the  calculation  usually  deviates  from  the  experimental  data.  On  the  one  hand, 
it  is  because  the  effect  of  leak  gap  and  separation  is  not  included  in  the 
calculation.  Another  reason  is  due  to  error  made  in  measuring  the  compression 
ratio. 

Becker used  the  molecular  drag  theory  to  analyze  the  leakage  problem.  How¬ 
ever,  this  is  not  solved  because  he  introduced  too  many  correction  factors  and 
the  molecular  drag  theory  itself  is  not  rigorous[2,3] .  in  this  work,  based  on 
equilibrium  equations,  the  relation  of  compression  ratio  versus  leak  gap  and 
impeller  separation  was  obtained  by  introducing  a  differential  circuit  model 
for  the  blades.  Then,  the  "dead  area"  correction  for  the  compression  ratio 
was  discussed.  Corrected  results  showed  that  the  effect  of  impeller  separa¬ 
tion  on  compression  is  not  as  high  as  predicted  by  Becker's  molecular  drag 
theory.  The  theoretical  values  also  agreed  very  well  with  the  experimental 
data. 

I.  Effect  of  Leak  Gap  on  Compression  Ratio  of  Single  Stage  Turbine 


In  order  to  ensure  the  safe  operation  of  a  turbomolecular  pump,  there  is 
usually  a  0. 5-1.0  mm  gap  between  the  rotating  part  and  the  stationary  part  of 
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the  pump.  This  gap  becomes  a  leakage  path  which  lowers  the  compression  ratio 
of  the  turbine.  Because  draw  rate  is  defined  as  the  volume  of  flow  per  unit 
time  when  there  is  no  pressure  difference  across  the  turbine,  therefore,  the 
leak  in  principle  will  not  affect  the  pumping  rate. 

Figure  1  shows  a  schematic  diagram  of  a  cross  section  of  a  single  stage  tur¬ 
bine.  From  the  figure  one  can  see  that  the  leak  gap  A1  is  equivalent  to  a 
short  circuit  in  the  pumping  line. 

In  equilibrium,  the  pumping  rate  is  zero  and  the  following  equation  is  valid: 


N1AZn+NlAA^i=N1AAZ,  (1) 

where  and  N2  are  the  number  of  molecules  incident  upon  a  unit  area  of  the 
blade  per  unit  time  on  the  inlet  and  exhaust  side,  respectively.  E12  and  l2 1 
are  the  forward  and  reverse  transport  probabilities,  respectively.  t*ie 

transport  probability  of  the  leak  gap.  A  and  AA  represent  the  suction  area 
and  leakage  area,  respectively. 

From  equation  (1)  we  can  find  that  the  compression  ratio  of  the  turbine  is 


-Ni_=  j  4-  K- 

1  N,  1  +  ,  ,  AA1,  ( 2 

1  +  AS 


where  K=I12 
221 

Therefore,  as  long  as  the  transport  probability  of  the  leak  gap  is  known,  the 
relation  between  the  compression  ratio  of  a  single  stage  turbine  and  the  leak 
gap  can  be  quantitatively  determined.  As  a  matter  of  fact,  with  the  exception 
of  increased  leak  gap  Ap  the  situation  is  identical  to  that  of  finite  blades 
without  leakage.  In  approximation,  we  consider  that  this  is  a  problem  to  add 
a  thin  wall  flow  conductance  with  an  area  AA  to  the  finite  blade  leak-free 
situation,  i.e.  Eq=l.  The  compression  ratio  in  equation  should  be  understood 
as  that  of  the  finite  blade [2 s3] 


K,  =  Kl+!1 


where  ^  =  Yl’  ^  t'ie  ^^-stance  between  two  neighboring  blades,  and  1  is 

the  length  of  the  blade. 
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Figure  1.  Cross-section  of  A  Single  Stage  Turbomolecular  Pump  Impeller 

From  equation  (2)  we  know  that  the  leak  will  cause  the  compression  ratio  to 
go  down.  With  the  same  leak  gap,  its  effect  on  a  higher  compression  ratio 
turbine  is  larger. 

Let  1=12  mm,  1=0.5  mm,  angle  of  inclination  of  the  blade  6=20°,  specific  speed 
3=0.5,  aspect  ratio  t/b-1.0,  and  blade  width  b=10  mm.  The  corresponding  trans¬ 
port  probabilities  are [5]  From  equations  (2)  and  (3) 

212  =  0.522,  2j,  =  0.262 

we  can  determine  that  K=1.99,  Kf=2.2  and  K^=2.03,  i.e.  the  leak  lowers  the 
compression  ratio  by  8%. 

II.  Effect  of  Leaks  in  Multi-stage  Turbomolecular  Pump  on  Its  Compression 
Ratio  -  Differential  Circuit  Model  of  the  Blade 

Figure  2  shows  the  leak  in  a  combination  turbomolecular  pump.  The  leak  of  the 
moving  turbine  is  on  the  outer  side  and  that  of  the  stationary  turbine  is  on 
the  inner  side. 


pumping  direction 


Figure  2.  Cross-section  of  a  Combination  Turbomolecular  Pump 

Because  neighboring  impellers  are  very  close,  the  pressure  drop  caused  by  the 
radial  flow  of  the  gas  molecules  along  the  gap  can  no  longer  be  neglected. 

For  this  reason,  the  pressures  on  both  sides  of  an  impeller  are  no  longer  con¬ 
stants.  Instead,  the  compression  ratio  of  the  turbine  is  a  function  of  the 
blade  length.  In  the  following  we  will  introduce  a  differential  circuit 
model  to  discuss  this  problem. 


Let  us  assume  that 


(1)  the  geometric  and  operating  parameters  along  the  blade  are  identical,  i.e. 
pumping  capability  and  forward  and  reverse  transport  probabilities  are  the 
same;  and 


(2)  the  radial  flow  conductance  of  the  leak  gap  is  the  same. 


As  we  know  that  a  vacuum  system  and  an  electronic  circuit  have  the  following 
corresponding  relations  as  shown  in  the  following  table. 


vacuum  system 


electronic  circuit 


flow  conductance  F 
pressure  P 
flow  rate  U 
vacuum  pump 


electrical  conductance  Y 
voltage  V 
current  I 
power  source 


Let  the  pressure  intensities  on  the  inlet  and  exhaust  side  of  turbine  A  in 
Figure  2  be  Pj_  and  P2(X),  respectively.  The  spacing  between  the  pumping 
elements  (i.e.  the  space  between  two  neighboring  blades)  is  t  =  1.  Based  on 
the  theory  of  transport  probability,  for  any  differential  element  dx  of  each 
pumping  element  along  the  length  of  the  blade  we  have 


NJV •dx  =  N1Zi2dx—N2Sildx  (  4  ) 

where  W  is  the  coefficient  of  the  differential  pumping  element  dx. 

When  the  compression  ratio  is  calculated,  the  effective  drawing  rate  is  zero. 
Hence,  the  net  flow  rate  corresponding  to  the  differential  pumping  element 
should  be  equal  to  the  increment  of  the  radial  flow  rate  at  the  gap,  i.e. 


dUi(x )  =?!*—  va,Wdx=Pi»-^  Vo^n'dx—  P2(x)'  ^j-v,'2udx  (  5  ) 

4  4  4 

where  U^(x)  is  the  radial  flow  rate  assigned  to  each  pumping  element  by  the 
gap  and  va  is  the  average  thermal  velocity  of  the  gas  molecule. 

For  convenience,  we  define 


Fn  — j 

Fzi  —  —vJZi  i 

where  F-^  and  F2^  have  the  dimension  of  flow  rate. 
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Yt 


Figure  3.  Differential  Circuit  Model  of  an  Edge  Turbine  Blade 
Substituting  equations  (6)  and  (7)  into  equation  (5),  we  get 


<*2l=P)Fn4-P2<*>-FM.| 
d(x)  l  1 

where  Un—P\(.Fn,~~Fii) 

which  is  the  maximum  pumping  rate  of  the  element.  When  the  pressure  on  the 
inlet  side  remains  unchanged,  U0  is  a  constant. 

According  to  equation  (8) ,  we  can  make  an  equivalent  differential  electronic 
as  the  one  shown  in  Figure  3  for  the  pumping  element.  The  parameters  of  the 
circuit  in  Figure  3  and  those  of  the  vacuum  system  are  related  as  follows: 


(8) 

(9) 


dIt-!f-Uot  dYz  —  —~»Fui 

dY,  =  ~L’Fii  Y,  =  F,i 

It(x)=Ut(x)>  7,=P„ 

Vz(x)  =  P2(*)e 

From  the  equivalent  circuit  shown  in  Figure  3,  we  can  list  the  following 
equations 

dU(x)  =  dI<>-CVz(x)-Vl')'dY{)  (i0) 

dV2(x)  (11) 

and  condition 

Ji(0)=  0  (12) 
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By  re-writing  equation  (10)-(12)  into  vacuum  parameters,  we  get 


*&,(*}_  =VjL-c?tlx)  -Pt3  ^ 
d(x )  /  * 

dPz(x)  _  Ui  (x) 
dx  l*F  i 

Ui(  0  )  =  0 

By  differentiating  equation  (13)  and  substituting  it  into 


(13) 

(14) 

(15) 

(14) ,  we  get 


dWiix)] 

dxi 


(16) 


Then,  from  equations  (15)  and  (16),  U^(x)  is  determined.  Afterward,  U^(x)  is 
plugged  back  into  equation  (14)  to  determine  P2(x).  Then,  the  compression 
ratio  can  be  determined  by 


.P,U>  _  E  . 

P‘  ■  (■yW^F.+FrcthjiuyhJff 


(17) 


By  averaging  equation  (17)  with  respect  to  the  blade  length,  we  get 


Ri  =  K  — 


K -  1 


F zi  j_ 

/?2i  •  ct  h 

/Ftl 

Fi  ^ 

1  Fi  \ 

F, 

Because  any  practical  pumping  element  satisfies 


cth^ 

Equation  (18)  can  be  simplified  as 

R,=K—  K~  1 


1  , 


F^i  + 

f  Fn 

Fi  "V 

F, 

(IS) 


Fn  2  ■  ,  we  have 

Ft 


(19) 
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In  the  following,  we  will  discuss  a  few  special  cases: 


1.  The  impeller  spacing  h=°°,  i.e.  a  single  blade  turbine, 
radial  flow  conductance  of  the  gap  is  infinite,  i.e.  Fj=°°. 
second  term  of  the  denominator  in  equation  (18)  becomes  , 


In  this  case,  the 
Therefore,  the 


The  compression  ratio  of  a  pumping  element,  which  is  the  compression  ratio  of 
the  impeller,  is 


K,=K- 


F,(K- 1) 
Fn  +  Fi 


=  1  + 


K—  1 


1  + 


AAZ, 

A‘2  2i 


which  is  identical  to  equation  (2) .  Notice  that  equation  (19)  does  not  apply 
to  this  special  case,  therefore,  it  cannot  be  used  to  find  the  compression 
ratio. 


2.  The  impeller  spacing  h=0,  i.e.  the  gap  of  a  combination  turbomolecular 
pump  is  the  minimum.  In  this  case,  F^=0.  By  substituting  it  into  equation 
(19)  we  get 


R,=K 

which  shows  that  the  leak  has  no  effect  on  the  compression  ratio.  Our  under¬ 
standing  of  this  conclusion  is  that  when  the  separation  is  zero,  the  effect  of 
the  leak  is  concentrated  at  the  gap.  It  is  not  spread  outward  and  the  averag¬ 
ing  effect  can  be  neglected. 

The  middle  turbine  shown  in  Figure  2  can  be  expressed  by  the  equivalent  cir¬ 
cuit  shown  in  Figure  4.  According  to  electronic  circuit  theory.  Figure  4  is 
equivalent  to  Figure  3.  By  changing  the  electrical  conductance  dY-p  to 
1/2 (dYi)  in  Figure  3,  we  find  that  the  radial  flow  rate  Ui(x)  and  the  pressure 
on  either  side  of  the  impeller  are  identical  to  those  in  Figure  4.  However, 
the  radial  pressure  drop  P2(0)-P(x)  on  the  exhaust  side  is  increased  by  one 
fold.  In  order  to  distinguish  them,  the  voltage  in  Figure  3,  i.e.  the  pres¬ 
sure,  carries  an  *  when  we  use  the  equivalent  circuit  in  Figure  3  to  solve 
the  problem.  The  real  pressure  intensities  on  both  sides  of  the  impeller 
should  be 


P2(x)  =  P2*(*)+!-CP2*(  0  )-P**U)3 

(20) 

P,<*>  =  P»*+-5-CP**(0)-P.*(*>3 

(21) 
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(22) 


Hence,  the  compression  ratio  is 


K,(x)=i‘4-5-=-  2  - 


Pi(x)  p* 0 
2 


for  ease  of  use,  equation  (22)  can  be  simplified  as 

P2*(x) 


K,(x) 


Pi* 


(23) 


Obviously,  the  condition  for  equation  (23)  to  be  valid  is 


P2*(x )  »-^-CP2*(  0  )-Pt*(x)3 

Zi 

Pi*»:-J-W(0)-Pt*(*)D 

Z 

Most  practical  turbines  can  meet  the  above  condition.  Therefore,  we  can 
determine  the  compression  ratio  by  directly  using  equation  (17). 


(. K -  1  )  F,  .ch 


K,(x)=K- 


U2-ti) 


By  averaging  equation  (24)  with  respect  to  the  blade  length,  we  get 


(24) 


Ki=K- 


_ (K-l)Fi 

Fzl+F,j2lucthJZl*i 
V  F,  V  F( 


(25) 


dYi  dYi  Vs(X)  d Yi  dYi 


Figure  4.  Differential  Circuit  Model  for  the  Middle  Turbine  Blade 
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Similarly, 

therefore, 


because  a  practical  pumping  element  satisfies  the  condition  /  pu 
equation  (25)  can  also  be  simplified  as  v  Ft 


>2 


R,  =  K- 


K-  1 


F_u 
F , 


u,  l2F 

\  V  F( 


21 


When  the  spacing  h  =  00 ,  from  equation  (25)  we  get 

K-l 


K,=  1  +- 


1  + 


AA2  x 
Al 


At*  21 

When  the  spacing  h  =  0,  from  equation  (26)  we  get 


(26) 


Rt  =  K 

which  is  identical  to  turbine  A  on  the  edge. 

We  will  discuss  the  choice  of  the  value  of  the  flow  conductance  From 

earlier  analysis  we  understand  that  the  radial  flow  in  the  gap  due  to  leakage 
is  an  internal  motion  of  the  gas  molecules.  It  is  not  limited  by  the  inlet 
aperture.  Therefore,  its  flow  conductance  should  be  calculated  by  using  the 
formula  for  a  long  and  narrow  slit [6] 


F,= 


8  /  KT \V 

3  \2xm  )  l 


(27) 


where  tQ  is  the  sum  of  the  blade  separation  t  and  thickness  on  the  cross- 
section  of  the  blade;  i.e. 


t0  =  t+- 


sin  9 


Let  the  angle  of  inclination  of  the  blade  be  6-20°,  the  blade  length  1-20  mm, 
the  leak  gap  Al  =  0.5  mm,  impeller  separation  h  =  2  mm,  the  blade  width  b  = 

10  mm,  and  the  thickness  of  the  blade  be  neglected.  The  distribution  of  com¬ 
pression  ratio  along  the  length  of  the  blade  obtained  based  on  equations  (17) 
and  (24)  is  shown  in  Figure  5.  Curves  1  and  2  represent  the  middle  and  edge 
turbines,  respectively.  The  dotted  lines  on  the  top  and  bottom  represent  the 
two  special  cases  that  h  =  0  and  h  =  °°.  From  Figure  5  we  can  see  that  the 
compression  ratio  of  the  middle  turbine  is  higher  than  that  on  the  edge  when 
the  impeller  separation  is  the  same.  Furthermore,  the  leak  affected  area  is 
also  smaller. 
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In  the  above  discussion  we  did  not  discuss  the  mutual  interaction  between  two 
neighboring  leak  gaps.  In  reality,  the  spacing  between  turbomolecular  pump 
impellers  is  very  small,  usually  less  than  2  mm.  Therefore,  leak  affected 
areas  are  mostly  concentrated  near  the  gap  (see  Figure  5).  Leak  affected 
areas  of  neighboring  turbines  almost  do  not  overlap.  Thus,  their  mutual 
interaction  can  be  neglected. 


Figure  5. 
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Figure  6.  Axial  Cross-section  of  a  Practical  Turbomolecular  Pump  Impeller 
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III.  Correction  for  the  "Dead  Space"  in  the  Gap 

In  addition  to  the  leakage  effect,  the  impeller  gap  also  affects  the  relative 
velocity  of  the  incident  molecule  relative  to  the  impeller. 

Figure  6  shows  the  axial  cross-section  of  a  practical  impeller,  h^  and  h2 
represent  the  inner  and  outer  spacings  of  the  impeller,  respectively.  For  the 
moving  impeller,  incident  molecules  are  composed  of  molecules  scattered  from 
surfaces  AB,  AC  and  BD.  Among  them,  molecules  scattered  from  AB  and  BD  have  a 
relative  velocity.  This  means  that  the  relative  velocity  of  molecules  scat¬ 
tered  from  AC  with  respect  to  the  blade  with  a  linear  velocity  v^  is  zero. 
Furthermore,  only  one— half  of  the  molecules  scattered  from  AC  and  BD  can 
incident  upon  the  moving  impeller.  Therefore,  the  average  relative  velocity 
of  incident  molecules  should  be 


Vd  — 


2l  +  hl 
2l~\~hi  +  h2 


•Vb  —  P'Vb 


(28) 


2l  +  hx  . 

where  P — 2l  +  hk+lh^  is  the  correction  factor. 

From  equation  (28)  we  know  that  the  presence  of  the  gap  lowers  the  relative 
velocity  of  the  incident  molecule.  For  the  fixed  impeller,  all  we  have  to  do 
is  to  exchange  hi  and  h2  in  equation  (28) . 

This  correction  term  is  required  because  of  the  loss  of  relative  velocity  of 
the  molecule  scattered  from  AC.  Therefore,  we  call  it  the  "dead  space"  cor¬ 
rection.  This  "dead  space"  on  the  inside  of  the  moving  impeller  is  particular¬ 
ly  large  when  the  blade  is  made  by  a  twisting  technique,  which  makes  this 
correction  term  very  important. 

IV.  Overall  Correction  for  the  Compression  Ratio  of  a  Turbomolecular  Pump 

According  to  the  above  discussion,  it  is  concluded  that  the  correction  of  the 
compression  ratio  of  a  combination  turbomolecular  pump  should  also  include 
leak  correction  and  "dead  space"  correction,  in  addition  to  the  finite  blade 
length  correction.  Every  correction  term  has  been  derived.  The  question  re¬ 
mained  is  the  sequence  of  correction.  Based  on  the  physical  significance,  it 
is  more  reasonable  to  perform  corrections  in  the  order  of  dead  space  ,  finite 
blade  length  and  leak  correction.  However,  the  computation  work  load  is  larg¬ 
er  according  to  this  sequence.  In  order  to  save  time,  the  order  we  are  ^ 
taking  is  leak  correction,  finite  blade  length  correction  and  dead  space 
correction.  The  final  result  was  verified  experimentally. 

V.  Experimental  Apparatus  and  Results 

Figure  7  shows  a  diagram  of  the  experimental  apparatus.  The  experimental  pump 
comprised  of  a  5-stage  moving  impeller  and  a  4-stage  stationary  impeller  made 
of  an  aluminum  alloy.  The  dimensions  of  the  moving  impeller  are  shown  in 
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Figure  8.  The  blade  parameters  are:  0=30°,  —  =  0.9,  b=10mm,  l=25mm,  Al=0.5mm 
and  the  impeller  outer  diameter  D=169mm.  b  The  impeller  separation  is  ad¬ 
justable.  The  blade  parameters  of  the  fixed  impeller  are  the  same  as  those  of 
the  moving  impeller. 

In  the  experimental  apparatus,  a  combination  of  a  2x-8  mechanical  pump  and  a 
F-203  diffusion  pump  in  series  is  used  as  the  fore-pump.  The  bearing  of  the 
molecular  pump  was  lubricated  by  a  high  temperature  instrument  oil  which  has 
lasted  for  approximately  100  hours  in  use. 


Figure  7.  Experimental  Apparatus 

Key:  1  -  molecular  pump;  2  -  medium  frequency  power  supply;  3  -  Model  F-203 
diffusion  pump;  4  -  Model  2x-8  mechanical  pump;  5  -  photo  electric 
counter;  6  -  ionization  gauge;  7  -  needle  valve;  8  -  dessicator; 

9  -  cold  trap;  10  -  air  release  valve;  11  -  cooling  water  pipe. 


Figure  8.  Cross-section  of  the  Experimental  Impeller 
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The  vacuum  system  is  sealed  by  vacuum  rubber  o  rings.  The  pressure  on  either 
the  inlet  side  and  the  exhaust  side  of  the  impeller  was  measured  by  ionization 
gauges.  The  relative  error  of  the  two  vacuum  gauges,  including  calibrating 
the  amplifiers,  is  less  than  5  percent. 


In  the  experiment,  the  compression  ratio  was  measured  according  to  the  "rising 
pressure  method"  specified  by  the  International  Organization  of  Standardiza¬ 
tion.  The  gas  to  be  measured  is  dry  air.  In  order  to  ascertain  that  the 
impeller  is  operational  under  molecular  conditions,  the  primary  pressure  was 
always  less  than  1x10“^  torr.  It  was  discovered  that  the  use  of  a  liquid 
nitrogen  trap  for  the  primary  vacuum  gauge  had  no  effect  on  the  experimental 
results.  For  this  reason,  all  experimental  data  reported  was  obtained  without 
the  cold  trap.  The  impeller  separations  are  1,  2  and  3  mm.  The  impeller 
rotating  speed  is  2,400  rev/min.  The  relative  specific  velocity  is  3  =  0.446. 


Figure  9.  Compression  Ratio  vs.  Impeller  Separation  h. 

The  experimental  results  are  shown  in  Figure  9.  The  curve  is  the  theoretical 
value  after  a  comprehensive  correction,  which  agrees  well  with  the  experi¬ 
mental  data.  The  relative  error  is  less  than  10  percent  which  within  experi¬ 
mental  error.  Figure  9  also  shows  that  the  relation  between  compression  ratio 
and  impeller  separation  is  not  as  close  as  predicted  by  Becker. 

VI.  Conclusions 

The  leak  caused  by  the  gap  in  a  turbomolecular  pump  lowers  its  compression 
ratio.  The  effect  of  this  leak  on  the  middle  impeller  is  less  than  that  on  a 
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side  or  single  stage  turbine.  Both  corrected  and  measured  compression  ratio 
values  show  that  the  relation  between  compression  ratio  and  impeller  separa¬ 
tion  is  not  very  close.  The  current  turbomolecular  pump  has  a  1  mm  operating 
gap,  which  is  reasonable.  To  further  reduce  the  operating  gap  will  increase 
the  technical  difficulty.  The  compression  ratio,  however,  cannot  be  signi¬ 
ficantly  raised. 

This  work  was  completed  under  the  direction  of  Professor  Hua  Zhongyi  [5478 
0022  0001].  The  machine  shop  of  Zhejiang  University  provided  all  the  experi¬ 
mental  apparatus  and  testing  equipment.  The  experimental  work  was  jointly 
completed  by  the  comrades  of  molecular  pump  development  group  of  the  machine 
shop  of  Zhejiang  University.  Among  them,  comrades  Xu  Shouyi  [1776  1343  5030], 
Wang  Fujiu  [3769  4395  3773]  and  Zheng  Jianqing  [6774  1696  7230]  participated 
in  all  the  experimental  work.  Comrades  Wang  Demiao  [3769  1795  5379]  and 
Ren  Gaochao  [0088  7559  3390]  of  Zhejiang  offered  many  valuable  suggestions 
and  facilitated  the  experimental  work.  Comrade  Jiang  Xiangqi  [5637  4382  4388] 
of  Fudan  University  also  provided  beneficial  guidance  to  the  experimental 
work. 
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[Text]  Abstract:  A  theoretical  model  for  pulsed  active- 
passive  mode-locked  solid  state  lasers  is  presented.  The 
characteristics  of  these  lasers  can  be  solved  by  computer 
simulation  with  this  model. 

I.  Introduction 

Passive  mode-locked  pulsed  solid  state  lasers  are  widely  used  as  a  source  for 
picosecond  laser  pulses,  but  these  lasers  are  intrinsically  unstable.  A 
simple  method  for  improving  their  stability  is  to  add  an  active  modulator  in 
the  cavity.  Such  active-passive  mode- locked  lasers  are  first  reported  in 
1974  and  subsequently  studied  by  a  number  of  researchers.2  5  According  to 
the  reported  experimental  results ,  the  stability  and  repeatability  of  these 
lasers  are  considerably  better  than  the  pure  passive  mode- locked  lasers. 

Although  pulsed  active-passive  mode-locked  solid  state  lasers  have  been  in 
existence  for  10  years,  there  has  not  been  a  theoretical  model,  to  our 
knowledge,  that  explains  the  characteristics  of  this  type  of  laser.  Good 
theoretical  models6-8  have  been  developed  for  pure  active  mode-locked  lasers. 
Such  model  assumes  the  mode- locking  take  place  in  a  uniform  single  mode  field 
in  the  cavity  and,  after  repeated  passages  through  the  modulator,  a  pulse  is 
formed.  The  pulse  continuously  narrows  and  finally  becomes  a  picosecond 
pulse.  Without  considering  the  pulse  formation  process,  a  steady  state 
self-consistent  solution  in  the  cavity  may  be  found  directly.6 

Two  theoretical  models  have  been  advanced  for  pure  passive  mode-locked  solid 
state  lasers.  One  model  is  based  on  the  statistical  fluctuation. 9  In  this 
model  the  mode-locked  pulse  sequence  grows  out  of  the  initial  noise  fluctua¬ 
tion  and  the  strongest  pulse  is  selected  from  the  pulse  sequence.  This  model 
satisfactorily  explains  the  fluctuation  and  instability  of  the  laser  but  can¬ 
not  determine  the  shape  and  width  of  the  pulse.  The  other  model  is  based  on  a 
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quasi-steady  state  analytic  theory1 3-1 5  which  gives  the  shape  and  width  of  the 
pulse  but  fails  to  explain  the  fluctuating  nature  of  the  laser. 

In  order  to  explain  the  behavior  of  pulsed  active-passive  mode-locked  solid 
state  lasers,  the  two  totally  different  theories  above  (one  statistical  and 
the  other  deterministic)  must  be  combined.  This  is  a  difficult  task  and 
perhaps  also  the  reason  that  no  theoretical  model  is  yet  available  to  date. 

It  was  pointed  out  by  some  authors14’16  that  a  passive  mode-locked  laser  can 
be  treated  with  nonstatistical  quasi-steady  state  analysis  after  the  gain 
medium  and  the  saturable  absorber  enter  the  saturation  region.  Based  on  this 
idea,  Zhu  Zhenhe  von  der  Linde17  treated  the  passive  mode  locking  first  with 
the  fluctuating  model  and  then,  after  the  gain  saturation,  with  the  quasi¬ 
steady  state  analysis.  In  this  paper,  we  carried  this  method  further  and 
propose  a  theoretical  model  for  active-passive  mode-locked  solid  state  lasers. 

II.  Theory 

This  section  describes  the  function  of  the  active  modulator  and  propose  a 
theoretical  model.  Figure  1  shows  schematically  an  active-passive  mode-locked 
laser  consisting  of  the  cavity,  the  amplifying  medium,  the  saturable  absorber, 
and  the  active  modulator.  For  simplicity  we  shall  restrict  our  discussion  to 
attenuative  modulators.  For  an  ideal  modulator  the  amplitude  transmission 
coefficient6  is  given  by 

o(/)  —  exp[—  B  sin2  (2ff/m/)3 

and  the  transmissivity  of  the  modulator  is 

Tmoiii)  —  exp  [—  26  smJ(2jr/M/)] 

where  fm  is  the  modulation  frequency  and  2f  is  the  modulation  depth. 

As  shown  in  Figure  1(a),  the  output  mirror  with  a  reflectivity  R  and  the 
modulator  may  be  viewed  as  a  single  unit  having  a  periodic  reflectivity: 

> 

Rm od(0  -  Tni  (  /  -  ^  •  R  •  rnod  (  /  +  y) 

where  d  is  the  distance  between  the  modulator  and  the  output  mirror,  and  c  is 
the  speed  of  light.  Substituting  (2)  into  (3),  we  have 

2?mod(<)  *“  R  •  exp{  —  25[1  —  ATcos(4»/m/)]} 

(4) 

where  the  constant  cos^4jt/„,  *  * 


(1) 

(2) 

(3) 


If  we  view  the  Fabry-Perot  cavity  in  terms  of  an  equivalent  ring  cavity  (see 
Figure  1(b)),  then  the  transmissivity  of  the  equivalent  modulator  is 
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Figure  1.  Schematic  Diagram  of  a  Pulsed  Active-Passive  Mode-Locked  Solid 
State  Laser 


Tm,e(<)  —  exp{—  25  [1  —  ATcos(4jr/m/)]} 


(5) 


and  the  corresponding  equivalent  amplitude  transmission  coefficient  is 

—  exp{— 5[1  —  Kcos(4jrfm/)]}.  (6) 

Equation  (5)  shows  that  the  frequency  and  phase  of  the  equivalent  modulator 
are  the  same  as  those  of  the  actual  modulator  and  only  the  waveform  is 
slightly  different.  Figure  2  shows  the  transmission  function  of  the  equiva¬ 
lent  modulator  for  d  =  0  and  d  =  L/10  (L  is  the  cavity  length)  and  a  modula¬ 
tion  depth  of  26  =  0.25.  As  can  be  seen,  the  difference  in  the  modulation 
waveforms  is  very  slight  and  the  effect  of  d  on  mode-locking  is  expected  to 
be  small.  Most  of  the  calculations  are  therefore  conducted  for  d  =  0.  For 
d  =  0,  Equations  (5)  and  (6)  can  be  simplified  to  (7)  and  (8),  respectively; 

.  Tn.eCO  —  exp[—  45  sinJ  (2 xfmt)l  (7) 

«eff(0  —  «p[~  25  sin*  (2jt/m/)].  (8) 
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Figure  2.  Transmission  Function  Tm  e(t)  of  the  Equivalent  Modulator 

Modulation  depth  26  =  0.25,  curve  1  for  d  =  0  and  curve  2  for 
d  =  L/10 


We  again  make  the  assumption  that,  in  an  active-passive  mode- locked  laser,  the 
pulse  sequence  originates  from  the  initial  noise  fluctuation.  We  first  con¬ 
sider  the  problem  from  the  viewpoint  of  the  fluctuation  model  in  Reference  11. 
In  this  model,  only  the  variation  of  the  average  light  intensity  is  calculated 
for  the  initial  phase.  After  the  stimulated  emission  far  exceeds  the  spontane¬ 
ous  emission,  we  assume  an  initial  pulse  group  consisting  of  M  pulses.  The 
intensity  distribution  of  the  M  pulses  is  logarithmic10  and  the  magnitudes  of 
the  pulses  appear  in  a  random  order.  Now  we  add  an  active  modulator  in  the 
cavity.  Tm  e(t)  is  equivalent  to  a  "window"  function  and  each  round  trip 
passage  in  the  cavity  is  equivalent  to  one  multiplication  by  the  window  func¬ 
tion.  When  the  stimulated  emission  far  exceeds  the  spontaneous  emission,  the 
laser  field  has  made  thousands  of  passages  in  the  cavity  and,  because  of  the 
window  function,  the  order  of  the  pulse  magnitudes  is  no  longer  totally  random. 
Obviously  the  probability  for  a  large  pulse  to  appear  in  the  center  of  the 
window  is  greater  and  the  intensity  distribution  of  the  pulses  is  no  longer 
logarithmic. 

For  very  small  26,  the  effects  of  active  modulation  is  also  small.  For  a 
given  pulse  sequence  the  intensity  distribution  does  not  deviate  greatly  from 
a  logarithmic  function  and  the  order  of  the  pulse  magnitudes  remains  essen¬ 
tially  random.  Using  widely  different  pulse  distributions,  we  computed  the 
intensity  distribution  for  two  limiting  cases:  (1)  logarithmic  function 
of  the  noise  distribution,  and  (2)  logarithmic  function  multiplied  by  the  Nth 
power  of  the  window  function  (7).  Here  N  is  the  number  of  passages  of  the 
laser  field  in  the  cavity.  We  found  that,  as  long  as  6  is  sufficiently  small, 
all  the  calculated  results  are  very  close  to  each  other. 

As  6  increases,  the  behavior  of  pulse  development  gradually  changes  from 
random  to  deterministic  and  the  calculated  results  for  different  pulse 
distribution  begin  to  differ  considerably.  This  shows  that  the  positional 
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distribution  of  the  pulses  of  different  magnitude  can  no  longer  be  considered 
random  when  the  stimulated  emission  is  much  stronger  than  spontaneous  emission. 
It  is  easy  to  imagine  that,  for  large  6,  large  pulses  will  occur  at  the  center 
of  the  window  and  small  pulses  will  appear  at  the  edges  of  the  window.  We 
therefore  choose  a  distribution  when  the  pulse  magnitude  is  greatest  at  the 
center  of  the  window  and  the  pulses  gradually  become  smaller  away  from  the 
center.  The  intensity  distribution  function  of  such  a  pulse  distribution  is 
difficult  to  calculate.  Tentatively  we  take  it  to  be  a  logarithmic  function 
multiplied  by  the  Nth  power  of  the  window  function,  (7)  or  (5).  Naturally  the 
half-width  of  such  a  distribution  will  be  somewhat  smaller  than  the  actual 
half -width.  When  6  is  sufficiently  large  and  the  half -width  of  the  Nth  power 
of  the  window  function  is  much  smaller  than  the  half-width  of  the  logarithmic 
function  is  much  smaller  than  the  half -width  of  the  logarithmic  function  of 
the  noise  distribution,  it  may  be  shown  that  the  calculated  results  obtained 
by  modifying  the  fluctuation  model  and  taking  the  position  and  intensity  dis¬ 
tribution  above  are  the  same  as  that  obtained  by  using  the  pure  active  mode¬ 
locking  model.  (After  such  modification  the  fluctuating  model  has  in  fact 
lost  its  statistical  meaning.)  In  the  pure  active  mode-locking  model  the 
laser  field  begins  with  a  uniform  single  mode  field  in  the  cavity  and  makes 
many  passages  in  the  cavity.  In  each  passage  the  laser  field  must  go  through 
the  active  modulator,  the  amplification  medium,  and  the  saturable  absorber. 

Based  on  the  discussion  above,  we  propose  the  following  model.  In  the  early 
stage  of  mode-locking,  we  use  the  "modified  fluctuation  model"  and  calculate 
only  the  variation  of  the  average  light  intensity.  When  the  stimulated  emis¬ 
sion  far  exceeds  the  spontaneous  emission,  we  choose  a  collection  of  M  pulses 
and  take  the  intensity  distribution  of  the  pulses  to  be  the  logarithmic  func-- 
tion  of  the  noise  distribution  multiplied  by  the  Nth  power  of  the  window  func¬ 
tion  given  in  (7)  or  (5).  N  is  the  number  of  passages  already  made  by  the 
laser  field  in  the  cavity.  The  sizes  of  the  pulses  are  the  largest  at  the 
center  of  the  window  and  gradually  decrease  to  the  two  sides.  Calculations 
are  then  made  using  the  equation  given  in  Reference  11  with  an  additional 
active  modulation  term.  After  gain  saturation  is  reached,  calculations  are 

w  17 

made  using  the  quasi-steady  state  analytic  theory. 

Although  it  is  difficult  to  assess  how  well  the  model  describes  the  behavior 
of  the  pulsed  active-passive  mode-locked  solid  state  laser,  the  discussion 
above  has  shown  that  the  model  correctly  describe  the  two  limiting  cases: 

(1)  When  6  approaches  zero,  the  model  approaches  the  theoretical  model11  for 
a  pure  passive  mode-locked  laser;  (2)  when  6  becomes  extremely  large  and  the 
effect  of  the  saturable  absorber  becomes  correspondingly  small,  the  model 
approaches  the  theoretical  model8  for  a  pure  active  mode-locked  laser.  We 
therefore  conducted  computations  for  the  general  case  using  this  model  and 
compared  the  results  with  experiments. 

III.  Analytic  Model 

We  examine  the  behavior  of  a  laser  pulse  that  propagates  one  cycle  in  the 
ring  cavity  shown  in  Figure  1(b). 
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1.  Pulse  Parameters 


We  assume  the  electric  field  of  the  laser  pulse  to  be  of  the  following  form: 


«?(/)  - 


(9) 


where  E(t)  is  the  complex  amplitude,  and  fp  is  the  center  frequency.  For 
simplicity,  the  amplitude  is  normalized  so  that  the  pulse  intensity  (power 
intensity)  is  given  by 

/(0  -  I  E(0 1*. 


E  (t)  is  the  pulse  entering  the  amplification  medium,  the  peak  value  of  the 
pulse  occurs  at  5  =  0  and  Ei<f)  is  the  Fourier  transform  of  Ex(t). 

2.  Amplification  Medium 

It  can  be  shown  that  the  amplitude  gain6  of  a  uniformly  broadening  laser  is 

U1> 

where  ga  is  the  saturated  amplitude  gain  at  the  center  of  the  line  (f o )  and 
Afa  is  the  atomic  linewidth. 

If  E2(t)  is  the  pulse  emerging  from  the  gain  medium,  its  Fourier  transform  is 

Ei(/)-fjC/)-Si(/).  (12) 


3.  Saturable  Absorber 

Using  the  four— level  system  shown  in  Figure  3  in  describing  the  absorption 
characteristics  (t23  a  0,  *  0)  of  the  saturable  absorber,  we  deduced  an 

expression  for  the  absorption  coefficient1  * 

where  k0  is  the  unbleached  absorption  coefficient,  Is  is  the  saturation 
intensity,  x(i.e.,  t34)  is  the  relaxation  time  of  the  saturable  absorber. 

For  the  special  case  of  fast  relaxation  (t-*0)  ,  (13)  may  be  simplified  to10’18 


1  +  (/(<)//,)* 


(14) 
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Figure  3.  Four-level  Energy  Diagram  of  the  Saturable  Absorber 
The  amplitude  transmission  coefficient  through  the  saturable  absorber  is 

T «-  exp  [—  ~  *(*)] .  ( 15  ) 

Clearly  the  pulse  after  passing  through  the  absorber  is 

EM  “  T«»(0  •  EM.  (16) 

The  pulse  intensity  I(t)  appearing  in  Equations  (13)  and  (14)  is 

/(0-*|Bi(OI*  (17) 

where  b  is  the  beam  cross  section  in  the  gain  medium  and  the  saturable 
absorber. 

4.  Active  Modulator 

The  amplitude  transmission  coefficient  aeff(t)  of  the  laser  pulse  passing 
through  the  modulator  is  given  by  Equations  (8)  or  (6).  Since  the  laser  pulse 
suffers  a  certain  time  delay  (let  the  time  delay  of  the  pulse  envelope  be 
B)  in  passing  through  the  gain  medium  and  the  saturable  absorber,17  the  peak 
of  the  pulsers  is  no  longer  at  t  =  0  but  at  t  =  B  when  the  pulse  goes  through 
the  modulator.  The  pulse  coming  out  of  the  modulator  is  therefore 

E«(0  B)  *  CO.  (18) 

An  analytic  expression  for  B  cannot  be  deduced  easily  but  its  value  may  be 
obtained  by  simulation  calculation. 

The  ideal  modulation  frequency,  that  is,  the  modulation  frequency  in  the 
absence  of  detune,  is 
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(19) 


where  L  is  the  optical  length  of  the  cavity.  This  frequency  is  different  from 
the  detune-free  modulation  frequency6  of  an  active  mode-locked  laser 


f. 


+ 


(20) 


In  considering  the  effect  of  the  detune  error  in  the  modulator,  the  detune  is 
defined  as  the  frequency  shift  from  the  ideal  modulation  frequency: 

A (21) 

This  definition  is  consistent  with  the  one  given  by  Harris  and  McDuff.19  When 
fm  >  f mo ,  the  shift  is  negative,  and  when  fm  <  fmo»  the  shift  is  positive. 

5.  Completing  the  Ring 


Finally,  we  take  into  account  the  time  delay  Tr  —  2L/c  incurred  in  completing 
the  ring  and  introduce  an  effective  loss  coefficient  y ,  and  obtain  the  pulse 
after  one  cycle: 

E'M  —  J  r)  \E<(/  “  T').  (22) 


Here  y  includes  the  output  loss  of  the  output  mirror  and  all  the  loss  in  the 
cavity. 

6.  Change  in  Gain  Coefficient 

The  equation  satisfied  by  the  gain  coefficient  g  (g  =  2ga)  may  be  derived^ 
from  the  rate  equation^6  of  the  four— level  laser.  In  the  saturation  region 
the  equation  is  simplified  to 


AL-mm 

At 


2  a.l 
hv 


(23) 


where  a,  is  the  cross-section  for  stimulated  emission  and  hv  is  the  photon 

a 

energy. 

Let  the  light  intensity  and  the  corresponding  gain  in  the  cavity  on  the  kth 
passage  by  1^.  and  g^,  respectively,  then  it  is  easy  to  show  that 

2(7*41  Tjrr 

gk+x-gk-h  w*  (24) 

where  W k  is  the  pulse  energy  in  the  cavity,  given  by  the  following  equation 
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IV.  Modified  Fluctuation  Model 


The  calculation  steps  and  the  meaning  of  the  symbols  in  the  various  equations 
are  given  in  detail  in  Reference  11  and  will  not  be  repeated  here.  We  now 
describe  the  modification  of  the  model. 

Based  on  Section  II,  the  relative  intensity  of  the  pulses  in  the  initial  group 
of  pulses  is  as  follows 


(»-3,  4,  •••,  M ) 

(26) 

.  .  >■  ■  - —Li 

Inn 

C(M,  P), 

(27) 

Here  C(MP)  is  a  constant10’12  determined  by  the  number  of  pulses  M  and  the 
probability  P.  The  definition  of  N  is  already  given  in  Section  II. 

Equation  (14)  in  Reference  11  is  then  modified  to  be 

'<«•'  - '« + '«  h  - r  -  TTte  -  M  ('  •  -iv1)] 

(i  —  1,  2,  M).  (28) 


Equations  (26)  and  (28)  are  valid  only  when  the  modulator  has  no  frequency 
shift  error.  In  the  presence  of  a  frequency  detune  error,  a  pulse  may  pass 
through  the  modulator  at  the  center  of  the  window  in  a  certain  passage,  but, 
as  the  pulse  circulates  in  the  cavity  repeatedly,  the  time  at  which  the  pulse 
passes  through  the  modulator  gradually  deviates  from  the  center  of  the  window. 
We  therefore  modify  Equations  (26)  and  (28)  as  follows 


and 


(29) 


h+i.i  “  h.i  +  ht.i  —  f  —  — ~7T - 45  sin3  I*  • 

>  1  -T  Oik  A  L 


+  bh.i 

0-1,  2,  M ). 


- 1 _ Lm _ 

2  M 


]} 


(30) 


when  Afm  =  0,  (29)  reduces  to  (26)  and  (30)  reduces  to  (28). 


By  combining  Sections  III  and  IV,  we  have  a  complete  theoretical  model. 
Analytic  solutions  are  extremely  hard  to  obtain  but  computer  simulations 
may  be  used  instead.  In  the  companion  article,  we  present  the  computer 
simulation  results.  The  theory  and  the  experiment  are  in  good  qualitative 
agreement,  suggesting  that  our  model  is  a  reasonable  one.21 
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[Text]  Abstract:  Computer  simulations  are  carried  out  using 
the  theoretical  model  presented  in  the  previous  paper  for 
pulsed  active-passive  mode-locked  Nd-glass  and  NdsYAG  lasers. 

The  dependences  of  some  characteristics  of  active-passive 
mode-locked  lasers,  such  as  mode-locking  threshold,  output 
stability,  pulse  width,  and  pulse  asymmetry,  on  the  modula¬ 
tion  depth  and  detuning  of  the  active  modulator  and  the 
relaxation  time  and  concentration  of  the  saturable  dye  are 
obtained  in  the  simulations.  The  theoretical  results  are  in 
good  qualitative  agreement  with  the  experiments.  The  simu¬ 
lations  explained  some  of  the  experimental  observations  and 
provided  further  understanding  of  the  physical  processes 
involved  in  active-passive  mode-locking. 

I.  Computer  Simulation 

We  have  conducted  computer  simulation  for  pulsed  active-passive  mode-locked 
solid  state  lasers  using  the  theoretical  model  presented  in  the  previous 
paper.1 

1.  Simulation  Calculation  Procedure 

First  we  calculate  the  development  of  pulses  and  the  increase  of  the  gain  for 
the  pulse  sequence  using  the  modified  fluctuation  model.1  When  the  gain  g 
reaches  its  maximum  value,  the  intensity  and  gain  of  the  largest  pulse  are 
reached.  Using  the  model  described  in  Section  III  of  the  previous  paper,  the 
largest  pulse  is  allowed  to  repeatedly  pass  through  the  laser  cavity  in  the 
simulation  and  its  evolution  is  followed.  The  maximum  intensity  and  gain 
recorded  earlier  are  used  as  the  initial  values  in  the  simulation,  and  other 
pulses  are  treated  as  a  background.  Using  Reference  2,  the  small  signal  solu¬ 
tion  is  first  obtained.  In  other  words,  we  intentionally  kept  the  pulse 
intensity  at  a  low  value.  In  this  case,  a  steady  state  solution  can  always 
be  found  after  the  pulse  makes  a  large  number  of  circulations  in  the  cavity 
no  matter  what  its  starting  pulse  shape  is.  The  small  signal  solution  is 
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used  as  the  initial  pulse  (Ei(t))  in  the  simulation  and  the  justification  is 
given  in  Reference  2. 

The  computer  simulation  requires  repeated  Fourier  transforms  and  inverse 
Fourier  transforms.  Here  we  use  the  pulse  envelop  E(t)  for  the  time  interval 
(-T  )  /2  ^  t  <  t(ij/2  and  substitute  the  function  E(t)  with  N  discrete  sampling 
values  En  =  I“%(tn),  where  tn  =  nAtj  n  =  -N/2  +  1.  N/2;  and 

At  =  Tu/N.  In  the  meantime  the  frequency  spectrum  E(f  -  fo)  of  the  pulse  is 
also  replaced  by  the  discrete  spectrum  En  =  I"'1  E(fn-  fo).  where  fn  =  fo  +  nfw; 
n  =  -N/2  +  1,  ••*,  N/2;  and  fu  =  1TW.  By  doing  so,  we  are  able  to  use  FFT 
(fast  Fourier  transform)  for  rapid  transformation  of  the  sampled  En  between 
the  time  domain  and  the  frequency  domain. 

The  protocol  of  the  computer  simulation  is  as  follows.  The  initial  pulse  En 
is  stored  numerically  and,  using  the  equations  given  in  Section  III  of  Refer¬ 
ence  1,  the  signals  after  passing  through  each  component  are  calculated 
until  one  cycle  is  completed.  The  signal  after  one  cycle  is  again  stored 
numerically.  Using  Equation  (24)  in  Reference  1,  the  new  gain  is  computed. 
Using  the  new  gain  we  compute  the  new  amplitude  gain  coefficient  (Equation  (11) 
in  Reference  1).  The  process  is  then  repeated  for  hundreds  of  times. 

2.  Aliasing 

In  the  FFT  calculation  we  assumed  that  the  discrete  sampling  values  are 
periodic  in  both  the  time  domain  and  the  frequency  domain  and  therefore 
subject  to  the  aliasing  effect.3*4  In  order  to  prevent  aliasing,  we  multi¬ 
plied  the  time  domain  signal  F(tn)  and  the  frequency  domain  signal  F(fn)^  ^ 
a  window  function.  We  have  tried  several  commonly  used  window  functions  » 
and  found  that  although  they  can  suppress  aliasing,  they  also  distort  the 
pulse  shape. 

Actually,  the  aliasing  effect  may  be  neglected  as  long  as  the  time  interval 
T  is  much  longer  than  the  pulse  width  and  the  frequency  bandwidth  Nfw  is 
much  greater  than  the  bandwidth  of  the  pulse.  We  used  a  Tw  equal  to  10  times 
the  pulse  width  and  found  that  aliasing  is  negligible  in  most  cases.  Under 
certain  circumstances,  however,  aliasing  may  be  negligible  in  the  first  few 
cycles  but  gradually  increases  to  such  a  degree  that  the  simulation  process 
must  be  absorbed.  Further  investigations  are  needed  to  find  an  effective  way 
for  eliminating  alias. 

We  need  a  T^  equal  to  10  times  the  pulse  width  instead  of  using  a  window 
function  in  our  simulation.  As  a  result  we  cannot  calculate  the  parameters 
that  would  lead  to  an  increase  of  aliasing.  Fortunately,  useful  results  can 
already  be  obtained  within  the  parameter  range  we  used. 

Like  in  Reference  2,  we  also  choose  N  =  512  in  our  calculation  and  the  results 
are  accurate  enough. 
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3.  Parameter  Values 


Table  1  shows  the  parameter  values  used  in  our  calculation  for  Nd  glass  and 
Nd:YAG  lasers. 

Tab le  1 . 


Medium 

O.l  hv 

Ccm’/J) 

*■1. 

(cm-) 

(MW/cm1) 

b 

T, 

(“0 

Nd  glass 

0.16 

200 

50 

I 

10 

NdtYAG 

4.27 

6 

50 

2 

10 

As  explained  in  Reference  6,  the  beam  cross-section  ratio  b  =  2  is  better  than 
b  =  1  for  a  Nd:YAG  laser,  we  therefore  take  d  to  be  2.  Depending  on  the  exper¬ 
imental  situation,  the  variable  parameters  are  the  loss  coefficient  Y,  the 
absorption  coefficient  Ko  of  the  saturable  absorber,  the  relaxation  time  T, 
the  modulation  depth  26  and  the  detuning  Afm. 

II.  Calculation  Results  and  Discussion 

1.  Mode-locking  Threshold 

The  mode-locking  threshold  changes  considerably  when  an  active  modulator  is 
introduced  to  a  passive  mode-locked  solid  state  laser.  Figure  1  shows  the 
dependence  of  the  mode-locking  threshold  on  the  modulation  depth  of  Nd  glass 
and  Nd:YAG  lasers  at  different  dye  concentrations.  In  both  cases  the  threshold 
decreased  but  by  different  amounts. 


Figure  1.  Mode-locking  Threshold  as  a  Function  of  Modulation 
Depth  (numbers  by  the  curves  are  the  Ko  values) 

(a)  Nd  glass,  y  =  0.35,  no  detune,  normalized  to  Ko  =0.7  and  6=0 

(b)  Nd:YAG,  y  =  0.2,  no  detune,  normalized  to  Ko  =  1.3  and  6=0 
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Under  the  action  of  the  modulator*  the  laser  field  in  the  cavity  soon  changes 
from  the  initial  random  noise  to  become  orderly  and  the  largest  pulse  is  quick¬ 
ly  chosen,  and,  understandably,  leading  to  a  lower  threshold.  As  the  modula¬ 
tion  depth  increases,  however,  the  total  cavity  loss  also  increases  and  the 
mode- locking  threshold  increases  again  after  reaching  a  minimum.  This  pheno¬ 
menon  is  evident  in  Figures  1(a)  and  1(b). 

Figure  1  also  shows  the  effect  of  the  dye  concentration  on  the  threshold. 

Within  a  certain  range,  the  mode-locking  threshold  of  a  Nd  glass  active- 
passive  mode-locked  laser  increases  as  the  dye  concentration  increases.  To  a 
certain  extent,  the  dye  concentration  does  not  affect  the  mode-locking  thresh¬ 
old  of  a  Nd:YAG  active-passive  mode-locked  laser.  In  general,  a  Nd-YAG  laser 
requires  a  higher  dye  concentration  than  the  Nd  glass  laser. 

Figure  2  shows  the  change  of  the  threshold  as  a  function  of  detune  for  differ¬ 
ent  modulation  depths .  As  expected ,  the  threshold  increases  with  detune . 

Also,  the  greater  the  modulation  depth,  the  more  the  threshold  increases  with 
detune.  This  implies  that  the  accuracy  requirements  on  the  modulation  fre¬ 
quency  and  cavity  length  are  greater  at  deeper  modulation  depth. 


Detune  (kHz) 

(a) 


Detune  (kHz) 

U>) 


Figure  2.  Mode-locking  Threshold  Versus  Detune  (Afm)  .  All  thresh¬ 
olds  are  normalized  with  respect  to  the  threshold  6=0 


(a)  Nd  glass,  y  =  0. 35 ,  k0  =  0.7 ,  26  for  curves  1,  2,  3 
and  4  are,  respectively,  0,  0.02,  0.06  and  0.2 

(b)  Nd : YAG ,  y  =  0.2,  k0  =  1.3,  26  for  curves  1,  2,  3  and 
4  are,  respectively,  0,  0.05,  0.15  and  0.3 

Figure  3  shows  a  set  of  normalized  SHER-n*  curves  for  Nd :YAG.  (SHER  stands 
for  Second  Harmonic  Enhancement  Ratio,  and  t|*  is  the  mode-locking  superthresh¬ 
old.)  As  can  be  seen,  the  mode-locking  threshold  becomes  extremely  sharp  with 
active  modulation,  a  unique  feature  of  the  active-passive  mode-locked  laser. 

At  a  high  dye  concentration,  the  passive  mode-locking  threshold  is  already 
very  sharp,  the  addition  of  the  modulator  makes  it  even  sharper,  but  not  as 
sharp  as  that  shown  in  Figure  3.  The  threshold  of  an  active-passive 
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V  ■  ■  ■■ 

Figure  3.  Normalized  SHER  Versus  N*  Curves  of  a  Nd:YAG  Laser  Under 
Peak  Output.,  Values  of  26  are  shown  on  the  curves, 
y  =  0.2,  Ko  “  1.0,  and  no  detune 

mode-locked  Nd  glass  laser  also  sharpens ,  but  not  as  much  as  that  of  the 
Nd:YAG  laser.  The  corresponding  curves  of  the  Nd  glass  laser  are  not  shown 
here . 

2.  Stability 

Active-passive  mode-locking  is  known  to  improve  the  stability  of  lasers,  again 
confirmed  in  our  calculations.  Figure  3  shows  that  the  mode-locking  window 
widens  with  increasing  modulation  (also  true  for  Nd  glass  lasers) ,  indicating 
an  increase  in  stability.  In  order  to  quantify  the  stability  improvement,  we 
calculated  the  percent  changes  in  the  output  pulse  strength  when  the  pimping 
rate  is  changed  by  3  percent .  We  call  this  percent  change  the  output  devia¬ 
tion. 

Figure  4  shows  the  output  deviation  as  a  function  of  the  modulation  depth. 
Figure  5  shows  the  output  deviation  versus  detuning.  Figures  4  and  5  show 
that  the  output  stability  of  an  active-passive  mode-locked  laser  improves 
with  modulation  depth  and  deteriorates  with  detuning. 

In  producing  an  active -passive  mode-locked  laser,  one  always  strives  for  a 
lower  threshold,  a  higher  stability,  and  a  lower  sensitivity  to  modulation 
frequency  and  cavity  length.  Our  modulator  is  not  the  greater  the  better, 
instead,  there  exists  an  optimum  depth,  the  same  is  true  for  the  dye  concen¬ 
tration.  These  optimal  values  depend  on  the  cavity  loss,  the  beam  cross- 
section  ratio,  the  relaxation  time,  and  other  specific  requirements.  General¬ 
ly  speaking,  the  modulation  depth  for  a  Nd:YAG  laser  should  be  greater  than 
that  for  a  Nd  glass  laser  and  the  dye  concentration  should  also  be  greater. 

In  addition,  our  calculations  also  showed  that  Nd:YAG  lasers  can  achieve  a 
higher  output  stability. 
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(a)  (b) 


Figure  4.  Output  Deviation  Versus  Modulation  Depth 

(a)  Nd  glass,  Y  =  0.35,  K0  =  0.7,  no  detune 

(b)  Nd : YAG ,  Y  =  0.2,  K0  =  1.3,  no  detune 


Figure  5.  Output  Deviation  Versus  Detune 


(a)  Nd  glass,  Y  =  0.35,  Ko  =  0.7,  26  -  0.1 

(b)  Nd:YAG,  y  =  0.2,  K0  =  1.3*  26  =  0.1 

3.  Pulse  Shape  and  Width 

The  pulse  shape  of  an  active-passive  mode-locked  laser  is  Gaussian 
(exp (-ln2  •  t2/t2)).7  Treating  the  dye  as  an  inertialness  saturable  absorber , 
the  small  signalPsolution  of  a  passive  mode-locked  laser  in  a  hyperbolic 
secant:  (sech2  •  (0.881  •£-) ) . 2 » 8  Calculations  show  that  the  small  signal 
solution  of  an  active-pasiive  mode-locked  laser  (with  inertialess  dye)  is 
still  a  hyperbolic  secant.  The  situation  is  somewhat  different  for  the  Nd :YAG 
laser  under  certain  conditions.  The  waveform  of  an  active-passive  mode-locked 
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laser  is  therefore  determined  primarily  by  the  saturable  dye  and  not  so  much 
by  the  modulator. 


Since  the  relaxation  time  of  the  saturable  dye  is  not  infinitely  small,  the 
small  signal  solution  is  an  asymmetric  pulse  with  a  slower  rise  and  a  faster 
drop.9  Such  a  pulse  may  be  represented  by  the  following  empirical  formula  in 
Reference  2. 


•$cchJ  [°;  (i) + (£)’  +.°;  (i)  1 (/  <-0) 


where  |a3|  «|a2|«  a3  and  |a3|  «  |a^|  «a{.  The  values  of  ai,  a2,  a3, 

a[,  a'2>  and  aj  depend  on  the  laser  parameters. 

Our  main  interest  is  the  true  output  pulse  shape,  that  is,  the  pulse  shape  when 
the  laser  pulse  intensity  reaches  a  maximum  in  the  cavity.  For  an  inertialess 
dye,  the  calculated  pulse  shape  is  basically  Gaussian.  For  a  NdtYAG  laser, 
the  tail  is  lower  than  that  of  a  Gaussian  only  when  the  modulation  depth  is 
large.  But  in  actuality  all  days  have  a  finite  relaxation  time  and  the 
computed  waveform  is  asymmetric,  as  shown  in  Figure  6,  with  a  steep-leading 
edge  and  a  more  gradual  tail.  Each  a  waveform  may  be  represented  by  the 
empirical  equation  (1) ,  but  the  signs  of  a2  and  a2  are  opposite  to  those  of 
the  small  signal  solution.  For  the  discussion  on  the  pulse  width  and  assymetry 
to  follow,  we  introduce  a  quantity  Ras  known  as  the  asymmetry  factor.  Figure  6 
shows  the  definition  of  the  pulse  width  Tp  and  the  asymmetry  Ras. 


Figure  6.  A  Typical  Asymmetric  Pulse  and  the  Definition 
of  the  Pulse  Width  (Tp)  and  the  Asymmetry  ®as 

In  the  simulation  calculation  we  observed  some  interesting  changes  of  the 
pulse  shape  in  the  cavity,  see  Figure  7.  The  pulse  shape  starts  with  a 
gradual  leading  edge  and  a  steep  tail  and,  just  before  the  pulse  intensity 
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Number  of  transits  in  cavity 

Figure  7.  Changes  of  Pulse  Shape  in  the  Cavity  of  an 

Active-Passive  Mode-locked  Nd  Glass  Laser  With 
a  Finite  Relaxation  Time  Saturable  Absorber. 

On  the  460th  pass  through  the  cavity,  the  pulse 
intensity  reached  a  maximum.  The  pulse  intensi¬ 
ties  are  arbitrary. 


Figure  8.  Pulse  Width  (solid  line)  and  Asymmetry  (dashed  line)  Changes. 

(Pulse  intensity  reaches  a  maximum  on  the  Nth^y  passage.) 

(a)  Nd  glass  laser,  dye  relaxation  time  T  =  10  ps, 
modulation  depth  26  =  0.06,  Nmax  =  460 

(b)  Nd :YAG  laser,  T  =  10  ps ,  26  =  0.1,  Nmax  =  76 

reaches  its  maximum,  the  asymmetry  reverses  and  the  leading  edge  becomes 
steeper.  (The  changes  in  Figure  7  are  for  a  Nd  glass  laser  but  a  Nd : YAG  shows 
the  same  asymmetry  reversal.)  The  cause  of  this  phenomenon  should  be  explained 
here.  For  a  saturable  absorber  with  a  finite  relaxation  time,  Equations  (13) 
and  (15)  of  Reference  1  should  be  used  in  the  calculation  of  the  amplitude 
transmission  coefficient.  This  transmission  function  has  two  effects  on  the 
simulated  pulse,  one  is  to  slow  down  the  decay  of  the  pulse  tail  and  the  other 
is  to  delay  the  pulse  nonlinear ly.  The  latter  dominates  when  the  pulse  intens¬ 
ity  is  small  and  the  result  is  a  longer  leading  edge  and  a  steep  tail,  like  the 
first  pulse  in  Figure  7 .  As  the  light  intensity  increases ,  the  dye  is  gradual¬ 
ly  bleached  and  the  first  effect  becomes  more  important,  causing  the  asymmetry 
to  reverse.  At  this  time  the  nonlinearity  of  the  dye  causes  rapid  changes  in 
the  pulse  width  and  the  asymmetry,  see  Figure  8.  As  a  result,  the  pulse  shape 
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of  the  output  pulse  sequence  changes.  The  early  pulses  have  a  gradual  rise 
and  a  steep  drop,  the  middle  and  late  pulses  have  a  steep  rise  and  a  gradual 
drop.  The  pulse  widths  are  also  different,  the  early  pulses  are  wider  than 
the  late  pulses. 

In  addition,  the  pulse  delay  caused  by  the  saturable  dye  is  nonlinear. 

Figure  9  shows  the  nonlinear  increase  of  the  delay  as  the  number  of  passages 
increases.  Figure  9  is  based  on  a  Nd:YAG  laser,  but  the  situation  is  simi¬ 
lar  for  a  Nd  glass  laser.  Strictly  speaking,  the  separation  between  two  adja¬ 
cent  pulses  is  not  a  constant,  although  the  difference  is  minute.  For  example, 
for  a  separation  of  10  ns ,  the  variation  is  no  more  than  a  few  picoseconds , 
too  small  a  change  to  detect  experimentally. 


Figure  9.  Pulse  Delay  Versus  Number  of  Transits  for  an  Active- 
Passive  Mode-locked  Nd:YAG  Laser  With  a  Finite  Relax¬ 
ation  Time  Saturable  Absorber.  The  pulse  intensity 
reaches  a  maximum  on  the  76th  transit. 


Figure  10.  Pulse  Width  (solid  line)  and  Asymmetry  (dashed  line)  Versus  the 
Absorption  Coefficient  of  the  Saturable  Dye 
T  =  10  ps,  26  =  0.1,  no  detune. 

(a)  Nd  glass,  y  =  0.35  (b)  Nd:YAG,  y  =  0.2 
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Relaxation  time  (ps)  Relaxation  time  (ps) 

(»)  <b) 

Figure  11.  Pulse  Width  (solid  line)  and  Asymmetry  (dashed  line)  Versus 
Relaxation  Time  of  the  Saturable  Dye 
26  =  0.1,  no  detune 

(a)  Nd  glass,  y  =  0.35,  Ko  =  0.7 

(b)  Nd : YAG ,  y  =  0.2,  K0  =  1.3 


<b) 


Figure  12.  Pulse  Width  (solid  line)  and  Asymmetry  (dashed  line)  Versus 
Modulation  Depth 
x  =  10  ps ,  no  detune 

(a)  Nd  glass,  y  =  0.35,  K0  =0.7 

(b)  Nd : YAG ,  y  =  0.2,  K0  =  1.3 
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Figure  13.  Pulse  Width  as  a  Function  of  Detune 
T  =  10  ps 

(a)  Nd  glass,  y  =  0.35,  K0  =  0.7,  26  =  0.06 

(b)  Nd : YAG ,  y  =  0.2,  K0  =  1.3,  26  =  0.15 

Figures  10  through  13  show  the  changes  of  the  pulse  width  and  the  asymmetry 
versus  the  concentration  and  relaxation  time  of  the  saturable  dye  and  the 
modulation  depth  and  detune  of  the  modulator.  One  observes  that  the  relaxa¬ 
tion  time  of  the  dye  has  the  greatest  effect  on  xp  and  Ras.  The  simulation 
results  show  that  the  saturable  dye  plays  an  important  role  in  the  pulse 
formation  process  and  affects  the  pulse  shape  and  width. 

The  following  formula  is  often  used  in  the  calculation  of  the  pulse  width10’11: 


where  xpa  is  the  pulse  width  of  a  purely  active  mode-locked  laser,  and  Tpp  is 
the  pulse  width  of  a  purely  passive  mode-locked  laser.  Based  on  this  equation, 
the  pulse  width  xp  of  an  active-passive  mode-locked  laser  should  be  less  than 
xpp,  but  this  contradicts  with  experimental  results.10  In  our  calculation,  we 
found  the  increase  in  pulse  width  of  an  active-passive  mode-locked  laser,  see 
Figure  12.  The  effect  of  the  modulator  should  cause  the  width  to  narrow,  so 
why  does  the  width  of  the  output  pulse  broaden?  In  analyzing  our  simulation 
process  we  found  that  introducing  an  active  modulator  lowers  the  mode-locking 
threshold  and  the  initial  values  of  the  gain  and  the  light  intensity  used  in 
the  calculation  in  Section  3  become  smaller.  As  a  result,  the  small  signal 
solution  pulse  width  increases  and  so  does  the  output  pulse  width.  If  the 
initial  values  of  the  gain  and  the  pulse  parameters  were  intentionally  set  to 
the  values  for  a  pure  passive  mode  locking,  then  the  calculated  results  show 
a  narrowing  of  the  output  pulse.  Our  computed  pulse  widths  for  a  Nd  glass 
laser  and  a  Nd : YAG  laser  are  generally  in  good  agreement  with  the  experiment. 
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4.  Distance  Between  the  Modulator  and  the  Output  Mirror 

In  order  to  study  the  effect  of  the  distance  between  the  modulator  and  the 
output  mirror,  we  carried  out  calculations  while  varying  d  from  zero  to  1/10 
of  the  cavity  length.  Results  show  little  changes  of  the  threshold,  stability, 
and  pulse  width,  thus  confirming  the  prediction  in  Reference  1  and  conforming 
the  experimental  result  in  Reference  10. 

III.  Conclusion 

We  have  proposed  a  theoretical  model  for  a  pulsed  active-passive  mode-locked 
solid  state  laser  and  conducted  computer  simulation  for  Nd  glass  and  Nd:YAG 
lasers.  Calculated  results  confirmed  experimental  observed  phenomena  of  a 
lowered  mode-locking  threshold,  an  improved  stability,  and  an  increased  pulse 
width  (over  that  of  a  purely  passive  mode-locked  laser).  The  theoretical 
results  are  in  good  qualitative  agreement  with  the  experiment.  Since  simpli¬ 
fying  assumptions  were  made  in  the  model ,  it  is  unrealistic  to  expect  quantita¬ 
tive  agreement. 

We  have  obtained  the  effects  on  mode-locking  by  parameters  of  the  saturable 
dye  and  the  modulator.  These  results  will  be  important  in  the  improvements 
of  an  active-passive  mode-locked  laser.  Our  model  successfully  described  the 
physical  process  of  active-passive  mode-locking.  With  the  addition  of  an 
active  modulator,  the  pulse  development  changes  from  a  random  behavior  to  a 
deterministic  behavior  as  the  modulation  depth  increases.  This  is  the  essence 
of  the  improved  stability. 
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OPTICAL  ADJUSTMENT  FOR  EOT2000  INFRARED  RANGEFINDER  DETAILED 
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25  Jun  85  pp  33-39 

[Excerpt]  In  the  article  "A  Study  of  the  Four-Axes  Problem  of  the  EOT2000 
Infrared  Rangefinder"  published  in  Vol  5  of  CEHUI  TONGBAO,  we  made  an  initial 
study  of  the  four  axes  relationship  of  EOT.  In  this  paper  we  discuss  the 
mechanism  and  the  optical  adjustment  of  the  EOT. 

The  optical  adjustment  of  the  EOT  differs  from  that  of  an  ordinary  theodolite 
and  involves  the  adjustment  of  the  optical  paths  of  the  four  axes. 

The  EOT  is  a  stand-alone  apparatus  and  the  disassembly  and  adjustment  must  be 
performed  strictly  according  to  a  set  procedure  to  insure  the  optimum  operat¬ 
ing  condition  of  the  apparatus.  Here  we  describe  the  disassembly  procedure 
for  the  infrared  optics  and  their  operating  principles. 
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Figure  1.  The  Optical  Micrometer 

Key: 

1.  Balance  spring  for  eliminating  gravity  effects 

2.  Micrometer  travel  stop 

3.  Glass  plate 

4.  Field  lens  for  reading  window 

5.  Fitting  prism 

6 .  Micrometer  rack 

7,11.  Double  micrometer  optical  wedge  assembly 

8.  Horizontal  axis  tilt  calibration  hole 

9.  Micrometer  handwheel  linkage 

10.  Optical  path  divider  plate 

12.  Lens  offset  correction  knob  fbr  H/V  changeover  prism  (13) 

14.  Rotation  axis  of  prism  13 

15.  Vertical  image  rotator  prism 

16.  Objective  assembly  for  horizontal  reading  microscope 
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Figure  2.  Structure  and  Function  of  the  Jogging  Mechanism 


Key: 

1.  Horizontal  axis  tilt  correction  plate 

2.  Three  set  screws  of  the  correction  plate 

3.  Four  set  screws  of  the  vertical  jogging  arm 

4.  Horizontal  clamping  plates 

6.  Horizontal  axis  set  screw  with  notch 

7.  Fixed  plate  with  a  conical  surface  and  connected  to  the  horizontal  axis 

8.  20.  Universal  joint 

9.  Axle  frame 

10.  Vertical  brake  tightening  adjustment  screw,  situated  on  the  shoulder  of 
part  (11) 

11.  Axle  with  a  notch  and  linked  to  the  universal  joint 

21.  Adjustment  screw  for  setting  the  brake  tension  of  the  aiming  mechanism. 
Adjustment  is  done  through  the  opening  (37) 

22.  Coupling  linkage  for  horizontal  brake 


74 


Key : 

1.  Base  for  the  automatic  compensator 

2 .  Pendulum 

3,  17.  Vertical  plate  illumination  prism 

4.  Ridge  prism  for  the  1:1  imaging  system 

5,  14.  Objective  assembly  and  base  for  the  1:1  imaging  system 

6,  21.  Base  for  the  large  optical  assembly 

7 .  Weight 

8,  20.  Infrared  reflector  and  base 

9.  Pendulum  agate 

10,  11.  Pendulum  amplitude  limiter 

12,  13.  Cylinder,  piston  and  piston  adjustor  for  air  damping 
15.  Pendulum  coefficient  adjustment  plate 

18.  Horizontal  dial  illumination  prism 

19.  Objective  assembly  for  vertical  plate  microscope 
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Figure  4.  Graduation  Method  for  the  EOT  Dial 
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Figure  5.  Horizontal  Axle  Structure 


Special  tool  to  disassemble  the  horizontal  axle  (M21  x  0.5) 

Infrared  reflection  system 

+.  Clamping  plates  and  six  set  screws  for  the  vertical  disk 
3.  Vertical  dial  and  protective  glass 
Infrared  reflection  lenses  and  their  base 
Set  screw  for  Part  2. 

Eight  set  screws  to  fix  the  horizontal  axle  and  adjustment  limit  screws 
Reading  microscope  system  rotates  toward  the  prism  base 
Horizontal  axle 
Telescope  body 


Figure  7.  Clamping  Ring 


Figure  8.  Principle  of  the  Suspended  Elastic  Pendulum 
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Figure  9.  Photo  Showing  the  Horizontal  and  Vertical  Axles  of  the  EOT  2000 
Being  Adjusted 

Key: 

1.  Single  angle  reflector 

2.  Prism  cover 

3.  Telescope 

4.  Cleaning  hole  for  the  under  surface  of  the  horizontal  dial  plate 

5.  Set  screw  of  the  base  (6)  of  the  optical  detector  (16) 

6.  Small  aperture  grating 

8.  Adjustable  spot  light  and  base 

9.  6V  2.5A  step-down  transformer 

11.  Power  pilot  light 

12.  Brightness  adjustment  knob  and  brightness  pilot  light 

14.  Electrical  plug  for  optical  detector  (16) 

15.  Second  cover  of  the  EOT  left  arm 
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a  is  the  [angular]  spacing 
between  two  graduation  marks 

Figure  10.  By  rotating  the  aiming  mechanism  and  by  locating  the 

position  of  the  cross  hair  in  the  bright  spot  coordinate, 
the  center  of  rotation  of  the  cross  hair  is  determined. 


Figure  11.  Infrared  Reflector  and  Its  Adjustment  System 

Key: 

1.  Infrared  reflector  base 

2.  Infrared  reflector  mirror  (2  mm  glass  plate  coated  with  transmission/ 
reflection  film) 

3.  4.  Limiting  screws  for  the  two  adjustments 

5.  Set  screws  for  base  1  and  base  7 

6.  Spring  reed 

8,  10.  2  mm  diameter  steel  spheres 

9.  Set  screw  for  the  mirror  base  and  the  telescope 
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CONFIGURATION  AND  LONGITUDINAL  WIND  TUNNEL  TESTING  OF  FORWARD-SWEPT  WINGS 

Mianyang  KONGQIDONGLIXUE  XUEBAO  [ACTA  AERODYNAMICA  SINICA]  in  Chinese  No  2, 

Jun  85  pp  61-68 

[Article  by  Ye  Weiliang  [0673  3555  2733]  and  Xiao  Renxi  [5135  0086  3556]  of 
the  Nanjing  Aeronautical  Institute:  "Analysis  of  Aerodynamic  Configurations  of 
Forward-Swept  Wings  and  Low-Speed  Longitudinal  Testing"] 

[Text]  Abstract 

Based  on  low-speed  wind  tunnel  testing,  calculation  of  load  distribution  and  the 
transonic-supersonic  area  rule,  this  article  analyzes  the  aerodynamic  character¬ 
istics  of  aircraft  with  a  configuration  of  forward-swept  wings  [FSW  hereafter]. 
It  suggests  measures  for  improving  the  aerodynamic  characteristics  of  FSW  and 
uses  the  principle  of  symmetry  to  solve  problems  related  to  the  application  of 
the  area  eddy  method  and  the  area  rule  calculation  method  for  FSW.  The  experi¬ 
mental  research  and  calculations  show  that  FSW  have  rather  good  aerodynamic 
characteristics  and  that  a  winglet  at  the  wing  tip  can  improve  the  lift  and 
lift-drag  ratios  of  FSW.  FSW  more  closely  approximate  an  optimum  load  distri¬ 
bution  than  do  aft-swept  wings  (ASW  hereafter].  There  are  obvious  improvements 
in  airflow  separation  at  the  FSW  base  from  using  a  close-coupled  canard  and  a 
fairing  panel  at  the  wing  base.  FSW  more  closely  approximate  an  optimum  load 
distribution  than  do  ASW.  The  axial  cross-sectional  load  distribution  on  an 
FSW  with  a  canard  more  easily  approximates  an  optimum  load  distribution  and 
thereby  can  result  in  an  obvious  reduction  in  zero-lift  wave  drag. 

I .  Int  ro due  t ion 

Aircraft  with  an  FSW  configuration  now  are  in  the  flight  trial  and  research 
stage.  The  American  X-29  aircraft  is  a  typical  example.  Although  the  use  of 
the  favorable  external  shape  of  FSW  to  overcome  transonic  drag  divergence 
appeared  at  the  same  time  as  ASW,  the  problem  of  divergent  static  aerodynamic 
elasticity  remained  unsolved  or  in  a  state  of  stagnation  for  a  long  period. 
Developments  in  composite  and  active  control  technologies  have  provided  a 
foundation  for  full  utilization  of  the  superiority  of  an  FSW  aerodynamic  con¬ 
figuration.  The  rearward  placement  of  the  wing  base  provides  flexibility  for 
the  overall  configuration,  while  a  matching  canard  makes  the  entire  aircraft 
more  closely  suited  to  transonic-supersonic  area  rule  laws  and  reduces  zero- 
lift  wave  drag.  Moreover,  trim  lift  is  high,  drag  is  low  and  the  angle  of 
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attack  at  stalling  is  large.  This  means  that  there  is  a  large  range  for  angle 
of  attack.-*-  Landing  and  motion  characteristics  are  good.  Separation  is  made 
difficult  at  the  wing  base  of  an  FSW  and  aileron  operation  characteristics  are 
good.  Its  load  distribution  along  the  direction  of  spread  approximates  an 
optimum  load  distribution  under  conditions  of  equal  weight.  FSW  have  positive 
benefits  for  meeting  aerodynamic  and  structural  requirements. 

The  low-speed  wing  tunnel  testing  in  this  article's  study  of  FSW  configuration 
aerodynamic  characteristics  included  comparative  tests  of  FSW  and  ASW  wing 
units,  manometric  testing  of  FSW  and  longitudinal  low-speed  dynamometric  testing 
of  FSW  with  canards  and  fairing  panels.  We  also  used  the  Woodward  Method^  to 
claculate  the  load  distribution  of  FSW  and  the  Eminton  Transonic-Supersonic 
Area  Rule  Method^  to  calculate  the  wave  drag  of  FSW  units.  The  calculations 
indicate  that  these  two  computational  methods  are  suitable  for  FSW. 

II .  Longitudinal  Low-Speed  Aerodynamic  Characterisitcs 

The  experiments  were  done  in  a  1  meter  inlet  diameter  low-speed  wind  tunnel. 

The  wind  speed  was  26.47  m/sec  and  the  Reynolds  number  Re  =  2.9  x  10^.  There 
were  four  models  in  the  group.  Figure  1  illustrates  the  external  shape  of 
one-half  the  dynamometric  models.  The  external  shape  of  the  half  of  the  mano¬ 
metric  model  was  the  same  as  model  No  3. 

Figure  1. 

External  Shape  of  Models 


1.  Flow  trend  analysis 

Fluorescent  microfibers  and  a  camera  were  used  to  illustrate  flow  trends  at  the 
wing  surface.  FSW  have  slip  motion  characteristics  in  the  direction  of  the 
wing  base  in  comparison  with  ASW,  as  shown  in  Figure  2  [not  reproducible]. 

During  no  air-flow  separation  at  the  wing  surface  with  a  low  angle  of  incidence, 
S-shaped  flow  lines  formed  effects  at  the  wing  base  and  wing  tip  of  the  FSW,  and 
the  aerodynamic  load  distribution  was  concentrated  at  the  base.  When  a  =  9°, 
separation  eddies  at  the  front  edge  were  very  obvious.  These  separation  eddies 
expanded  gradually  moving  from  the  wing  tip  to  wing  base.  The  axis  of  the 
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eddies  expands  with  the  angle  of  incidence  and  they  slope  toward  the  rear  edge. 
Moreover,  lateral  flow  increases  the  thickness  of  the  boundary  layer  in  the 
area  of  the  wing  base.  When  a  =  17°,  the  air  flow  at  the  wing  base  is  com¬ 
pletely  separated,  but  does  not  separate  at  the  wing  tip.  There  was  complete 
separation  at  the  wing  tip  of  an  ASW  at  this  time. 

Installation  of  a  canard  caused  an  obvious  improvement  in  the  flow  field  at 
the  base  of  the  FSW.  This  was  especially  true  of  a  close-coupled  canard. 
Photographs  at  different  angles  of  incidence  indicate  that  with  the  effective 
interference  of  canard  body-eddy  removal,  there  is  almost  no  lateral  flow  on 
the  main  wing  when  a  =  5°.  Front-edge  eddies  at  the  main  wing  caused  partial 
separation  at  the  outer  wing  segment  when  a  =  9°.  Local  separation  occurs 
only  in  the  middle  section  of  the  main  wing  when  a  =  17°.  Without  the  canard, 
the  region  of  separation  is  smaller  and  there  is  no  separation  at  the  wing  base 
or  tip.  Installation  of  the  canard  at  different  angles  is  even  more  effective 
in  improving  the  flow  field  of  the  main  wing. 

A  fairing  panel  at  the  wing  base  of  an  FSW  also  effectively  improves  the  flow 
field.  A  fairing  panel  at  the  wing  base  causes  the  air  flow  to  flow  laterally 
toward  the  outside  and  thereby  stops  boundary  layer  accumulation  formed  by 
lateral  flow  of  the  outer  wing  air  flow  toward  the  inside.  In  combination 
with  effective  interference  of  the  canard,  the  region  of  local  separation  is 
reduced  even  more,  as  shown  in  Figure  3  [not  reproducible], 

2.  Dynamometric  results 

a)  Lift  characterisitcs.  An  FSW  with  a  moderate  or  small  angle  of  incidence 
has  a  linear  lift  curve.  The  separation  eddy  destruction  at  the  front  edge 
when  the  angle  of  incidence  is  large  results  in  a  declining  slope  in  the  lift 
line  and  the  curve  starts  to  bend.  After  addition  of  a  canard,  the  effective 
interference  of  the  body-eddy  removal  by  the  canard  controlled  air  flow  separa¬ 
tion  at  the  wing  base  and  caused  the  linearity  of  the  lift  curve  to  improve. 

A  slight  curve  in  the  lift  curve  appears  only  after  the  angle  reaches  17°,  as 
shown  in  Figure  4.  The  slope  of  the  lift  line  of  an  FSW  is  greater  than  that 
of  an  ASW.  The  ability  of  eddy-body  removal  by  a  close-coupled  canard  to  con¬ 
trol  air  flow  separation  on  the  surface  of  the  main  wing  is  greater  than  that 
of  a  far-coupled  canard.  When  the  lift  curve  of  a  far-coupled  canard  starts 
to  bend,  the  lift  curve  of  a  close-coupled  canard  still  maintains  an  excellent 
linear  relationship.  For  an  FSW  with  a  fairing  panel,  a  small  fairing  panel 
has  better  lift  than  a  large  panel.  The  small  panel  causes  local  improvements 
in  air  flow  at  the  wing  base,  with  the  wing  surface  in  the  middle  section  of 
the  main  wing  being  subjected  mainly  to  canard  body-eddy  removal.  This  results 
in  rather  good  lift  characteristics.  A  fairing  panel  that  is  too  large  causes 
stronger  slip  of  the  air  flow  from  the  wing  base  toward  the  outside,  which  is 
not  good.  Selection  of  the  size  of  a  fairing  panel,  therefore,  deserves  dis¬ 
cussion.  A  winglet  at  the  wing  tip  caused  an  obvious  increase  in  lift  of  the 
FSW. 

b)  Lift-drag  characteristics.  FSW  have  less  drag  than  ASW  when  the  angles 
of  incidence  are  identical.  Because  the  lift  of  ASW  is  greater  than  FSW, 
however,  the  lift-drag  ratio  of  ASW  is  greater  than  that  of  FSW  within  a  range 
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Figure  4. 
Lift  Curves 


of  moderate  angles  of  incidence.  When  there  is  a  large  angle  of  incidence, 
the  more  rapid  air  flow  separation  of  ASW  causes  the  lift-drag  ratio  to  be 
nearly  the  same.  Canards  have  obvious  effects  in  improvement  of  the  lift-drag 
ratio  of  FSW,  while  the  results  are  smaller  with  ASW,  as  is  shown  in  Figure  5. 
This  is  because  of  the  greater  frontwash  effects  of  FSW  compared  with  ASW. 
Moreover,  the  canard  body-eddy  removal  improves  the  results  of  FSW  wing  surface 
flow  fields  rather  well.  When  the  angle  of  incidence  reaches  17°,  there  still 
is  no  obvious  air  flow  separation  on  the  main  wing,  whereas  rather  large  air 
flow  separation  regions  begin  to  appear  at  the  wing  tip  on  ASW,  causing  a  rapid 
increase  in  drag.  Canard  placement  affects  load-drag  characteristics.  Under 
moderate  angles  of  incidence,  the  lift-drag  ratio  of  close-coupled  canards  is 
greater  than  far-coupled  canards.  For  FSW  with  fairing  panels,  the  lift-drag 
ratio  of  small  panels  is  better  than  that  of  large  ones.  When  a  is  within  a 
range  of  less  than  10°,  drag  is  reduced.  The  improvements  in  lift-drag  ratios 
are  shown  in  Table  1. 


Figure  5. 

Lift-Drag  Ratio  Curves 


Cy/Cxl 


a  FSW  +  canard 
vASW  +  canard 


•  FSW^  +  wing^tLarge  panel 
Small  panel 


Far-coupled 
^  canard 
Near-coupled 


Table  1. 


a  (degrees)  1 

2 

3 

5 

7 

9 

%  decrease  in  25.1 

22.10 

17.56 

7.6 

1.69 

0,73 

drag 

%increase  xn  Q7 

lift-drag  ratio9*8' 

32.70 

23.35 

8.46 

3.00 

2.77 
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The  maximum  lift-drag  ratio  of  the  small  panel  was  Kjnax  6.37.  A  winglet  at 
the  wing  tip  improved  the  lift-drag  ratio  of  FSW. 

c)  Moment  of  force  characteristics.  When  an  ASW  has  a  large  angle  of  incidence, 
air  flow  separation  at  the  wing  tip  causes  a  tendency  to  loft.  FSW  do  not  have 
this  shortcoming.  Because  of  the  effects  of  separation  eddies  and  air  flow 
separation  at  large  angles  of  incidence  at  the  front  edge,  however,  the  rela¬ 
tionships  in  the  curves  n^Cy  have  poor  linearity.  After  installation  of  a 
canard,  moment  of  force  characteristics  are  improved.  The  improvement  is  most 
obvious  on  FSW.  When  there  is  a  large  angle  of  incidence,  the  tendency  toward 
lofting  without  an  ASW  still  remains.  The  fairing  panels  had  no  obvious  effects 
on  moment  of  force  characteristics. 

III.  Measurement  of  Pressure  Distribution  and  Load  Computations 

The  difference  in  flow  trends  between  FSW  and  ASW  causes  the  absolute  negative 
pressure  value  of  the  base  part  cross-section  to  be  greater  than  the  tip  part 
cross-section.  The  lift  at  the  wing  base  is  rather  concentrated  toward  the 
front  edge.  This  was  quite  evident  in  measurements  of  surface  pressure  distri¬ 
bution  on  the  wing  profile  when  a  =  11°,  as  shown  in  Figure  6a.  The  back¬ 
pressure  distribution  is  related  to  the  down  wash  distribution.  Naturally,  the 
eddies  caused  downwash  in  the  region  of  the  rear  edge  to  be  greater  than  at  the 
front  edge.  As  the  angle  of  incidence  is  increased,  the  absorption  peak  at  the 
front  edge  increases  in  the  beginning.  When  the  angle  of  incidence  becomes 
large,  air  flow  separation  eliminates  the  absorption  peak  at  the  front  edge  and 
increases  the  load  on  the  rear  edge.  Figure  6b  shows  that  the  pressure  co¬ 
efficient  on  the  No  2  [model]  cross-section  changes  as  the  angle  of  incidence 
is  increased. 

Figure  6.  Pressure  Coefficient  Distribution 
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Under  conditions  of  no  air  flow  separation,  the  surface  eddy  method  can  be  used 
to  calculate  the  load  distribution  of  an  FSW.  Independent  calculations  were 
done  on  an  FSW  using  a  manometric  model.  The  results  show  that  the  back  load 
distribution  conforms  fairly  closely  to  the  test  results.  When  the  angle  of 
incidence  is  small,  say  a  =  5°,  the  load  coefficient  ACp  becomes  less  than  0.06 
except  at  the  front  edge.  When  using  the  surface  eddy  method  for  computations 
related  to  FSW,  although  there  is  not  change  in  the  basic  shape  of  the  angle, 
the  sequence  of  mutual  decrease  in  the  angle  point  solution  at  normal  turbulence 
speed  is  the  opposite  of  ASW,  however.  The  left  wing  of  FSW  is  divided  into  n 
blocks.  The  angle  point  solution  for  the  surface  eddies  of  block  j  for  the 
normal  turbulence  speed  W  of  the  i  control  points  (see  Figure  7)  is 


In  the  formulas, 


f  +  (2.2) 

4  n 

(2.3) 

4«t  r 


Xx  x  +  0*  y 

fi\/  (x  —  X xy)1  +  l\  z 1 


(2.4 


M  i  ~  s  h  ~ 1 


x 

Tr 


(2,5 
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At  subsonic  speeds,  K— 1,  8^=1—  1^  —  X  ,  +  3  •  A,  is  the  tangent  value 

of  the  angle  of  sweep  of  the  front  edge  of  the  small  fairing  panel,  and  d  = 

\/ x^  +  $2r2,  while  r^  =  y2  +  z^. 

Determination  of  normal  flow  should  correspond  to  the  boundary  condition: 

C(^)„  +  Or*)i,lyi~(a+dzldx)Vm  (2.6) 

>  - 1 

The  surface  eddy  intensity  Y>  is  derived  by  solving  this  linear  equation.  The 
load  coefficient  of  the  number  i  panel  is: 

ACp,-  —Cpy  Cp±.  ~2y,'  (2.7) 

The  No  3  model  FSW  and  the  corresponding  ASW  were  used  for  calculation  of  the 
lift  distribution  in  the  direction  of  spread  for  FSW,  Cyl-b/CybA  and  the  lift 
coefficient  distribution  in  the  direction  of  spread,  Cy^/Cy.  It  more  closely 
approximates  an  optimum  load  distribution  under  conditions  of  equal  weight  than 
do  ASW.  When  compared  with  an  elliptical  distribution,  this  type  of  distribu¬ 
tion  has  a  load  distribution  that  is  concentrated  at  the  wing  tip.  The  results 
of  the  computations  are  shown  in  Table  2. 

Table  2.  Comparison  of  Angle  of  Spread  Lift  Distribution 


2 z/i  KAngle.of  sj 

^read® 

0.2 

0.4 

0.6 

0.8 

1.0 

C  y  b/CybA 

FSW 

1.350 

1.275 

1.095 

0.872 

0.601 

0 

ASW 

1.108 

1.095 

1.055 

0.956 

0.775 

0 

C  V /Cy 

FSW 

1.008 

1.050 

1.060 

1.012 

0.890 

0 

ASW 

0.840 

0,921 

1.032 

1.108 

1.133 

0 

IV.  The  Use  of  the  Transonic  and  Supersonic  Area  Rule  To  Calculate  the  Zero- 
Lift  Wave  Drag  of  a  Forward-Swept  Wing  Unit 

Because  there  often  were  certain  controlling  cross-sections  in  the  overall 
configuration,  the  optimum  cross-sectional  distribution  for  minimum  zero-lift 
wave  drag  for  an  aircraft  was  derived.  The  Eminton  method  assumes  that  x  =  0, 
x  =  1  and  that  x  =  Ki(i  =  l,2....n).  When  the  cross-sectional  area  is  assumed 
to  be  S(0)  =  N,  S(l)  =3  and  S(Ki)  =  A±,  the  S(x)  distribution  for  minimum  zero- 
lift  wave  drag  and  the  corresponding  minimum  zero-lift  wave  drag  value  Cxo  can 
be  derived.5  In  order  to  facilitate  computation,  the  length  of  the  fuselage  is 
assumed  to  be  1.  When  the  reference  area  is  1,  the  zero-lift  wave  drag  co¬ 
efficient  can  be  written  as: 

Cx „  =  — M‘  [V(*,).S'(*,)ln|*,-*,|</*.d*.  (3.1) 

2?r  J  o  J  o 
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The  above  formula  can  be  converted  to  a  numerical  solution  as: 

.  ;  co  jo/10  /'■  *'•  *• 

S'(*)  =  2  a'sin  <0<6<x)  (3.2) 

The  relationship  between  8  and  x  is : 

*=JL  (i  cos  6) 

A 

°r  —  —  I  S' (x)sin  r  6  d  6  r  =  1,2— 

By  such  conversion,  the  corresponding  area  distribution  of  minimum  zero-lift 
wave  drag  is : 

s <*>  =N+  (B-N)U(x)  +  2  *‘P<X>  *<>  (3.3) 

1-1  * 

In  the  formula,  ; 


U (x)  =  1  [cos’*1  (1  -  2x )  -  2(1  -  2x )  y/x(l  -  x)  ] 

ft 


n,~  ,  \  j  t  *+*,  ~2xx,  +2y/xx,(l-x)  (l-x,)  , 

/*<*,*,)—  T u - *.)•  I n  rrj;-2-^;-2^~;a-i7r7T-"7 + 


+  2  (*  +  *,■  -  2x  x,)s/ xx,  (l  -  x)  (l  -  x,) 

7  2 

When  x  =  x,  then  P(x,x^)  =  4xt(l  -  x^) 

The  in  the  formula  can  be  determined  through  the  following  linear  equation 
when  the  given  area  distribution  condition  x  =  x±  so  that  S(xi)  =  A,-  (i  - 
1,2 . n) : 


x,)  =  (A,~  IV)  ~  (B-N)U(Xi) ,  ;=1,  2  —  n 


t  -  t 


(3.4) 


After  \±  is  derived,  then  zero-lift  wave  drag  is: 


C,„  =  ±  (fl-/V)‘  +  *  2  i  T, 


I  -  l  I  -  l 


(3.5) 


The  remaining  question  is  how  to  derive  the  Mach  tangent  area  of  FSW  for  each 
corresponding  value  of  M.  The  Mach  oblique  tangent  plane  of  FSW  can  be  con¬ 
verted  to  the  former  point  on  the  symmetrical  ASW  for  calculation.  The  front 
edge  of  the  FSW  is  the  rear  edge  of  the  ASW.  The  function  f(v)  for  listing  the 
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thickness  of  the  wing  profile  is  the  opposite  of  that  for  ASW  and  there  is  an 
opposite  sequential  order.  At  the  same  time,  the  axial  coordinates  of  the 
derived  Mach  tangent  plane  also  are  symbolic  opposites. 

The  normal  projection  cross-sectional  area  for  the  Mach  tangent  plane  of  a 
certain  thickness  is  comparison  with  ASW  is: 

f  y upper 

A,  =-  2  \  z  tly  (3.6) 

J "  lower 

The  z  in  the  formula  is  the  half- thickness  of  the  cut  wing  area,  as  shown  in 
Figure  8. 


Figure  8. 

Mach  Tangent  Area 


Given  that  v  =  — ,  that  v0  =  ^  and  that  n  =  then  j  is  the  maximum  thick- 

c  „  j  to 

ness  ratio.  When  f(v)  =  is  inserted  into  formula  3.6,  then 


As~  +S0C»  (AT  +  Voy 

f’uppe/(v) 

J-  (K±vydv 

lower 


X 


(3.7) 


In  the  formula, 


K  =  ±m(i+  ASSII) 

\  tg  A  y 

m  =  p.  tgA , 

V-/Q 
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The  0  is  the  corresponding  meridian  angle  of  the  Mach  oblique  tangent  plane. 

The  A  is  the  angle  of  sweep  of  the  front  edge. 

The  signs  in  the  formula  are  set  by  having  the  symbol  above  correspond  to  the 
right  wing  and  the  symbol  below  to  correspond  to  the  left  wing.  The  integral 
upper  and  lower  limits  are  determined  according  to  actual  conditions.  When 
intersecting  with  the  front  edge,  v  -  0  and  when  intersecting  with  the  back 
edge,  then  v  =  1.  When  intersecting  with  the  side  edge,  vt  =  xt/Ct. 

As  for  the  body  of  rotation  of  the  fuselage,  the  normal  projection  of  the  Mach 
oblique  tangent  area  can  be  computed  using  the  normal  tangent  cross-section. 
When  the  aircraft  wings  are  at  M  >  1,  the  cross-sectional  area  is  the  average 
value  of  the  normally-projected  Mach  oblique  tangent  plane  corresponding  to 
each  meridian  angle.  In  this  way,  we  can  calculate  the  zero— lift  wave  drag  of 
an  FSW  aircraft  fitted  with  a  canard  and  the  equivalent  body  of  rotation  cross- 
section.  Moreover,  it  can  be  compared  with  ASW.  The  fact  that  an  FSW  config¬ 
uration  more  easily  approximates  an  optimum  cross-sectional  area  means  that  its 
has  less-zero-lift  wave  drag  than  do  ASW  configuration  aircraft.  The  results 
of  the  calculations  using  model  No  3  are  that,  when  M  =  1.2,  then  the  FSW 
CXo  =  0.02009  while  the  ASW  Cxo  =  0.02425,  so  it  is  20  percent  less.  The 
equivalent  cross-sectional  distribution  is  shown  in  Figure  9. 


Figure  9. 

Equivalent  Cross-Section  Along 
the  Axial  Distribution  At  Mach 
1.2 


V.  Conclusions 


1.  Separation  occurs  easily  at  the  wing  base  of  FSW.  Flow  is  rather  good  at 
the  tip,  however.  FSW  have  less  drag  and  lift  than  ASW  at  identical  angles 
of  incidence.  The  lift-drag  ratio  of  FSW  is  less  than  that  of  ASW  at  moderate 
angles  of  incidence.  A  winglet  at  the  wing  tip  can  improve  the  lift  and  lift- 
drag  ratio  of  FSW.'7 


2.  Canards  have  obvious  benefits  in  improvement  of  air  flow  separation  at  the 
wing  base  of  FSW.  Close-coupled  canards  have  better  deflectional  results. 

When  a  canard  is  installed,  the  degree  of  improvement  in  the  lift-drag  ratio 
of  FSW  is  greater  than  for  ASW.  Canards  have  obvious  benefits  in  maintaining 
a  linear  relationship  with  moment  of  force. 

3.  A  fairing  panel  can  improve  the  flow  field  at  the  wing  base  of  an  FSW. 

The  fairing  panel  matched  with  the  canard  should  be  of  appropriate  size.  Small 
panels  had  a  higher  lift-drag  ratio  on  FSW  in  the  tests  than  did  large  panels. 
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4.  There  is  a  great  load  on  the  wing  base  of  FSW  and  lift  is  rather  concen¬ 
trated  at  the  front  edge.  Computations  indicate  that  the  load  in  the  direction 
of  spread  more  closely  approximates  an  optimum  distribution  than  for  ASW. 

5.  A  canard  configuration  on  FSW  facilitates  a  distribution  that  achieves 
the  transonic-supersonic  area  rule.  According  to  calculations  for  model  No  3, 
the  zero-lift  wave  drag  can  be  reduced  by  20  percent  compared  with  ASW. 
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JPRS-CST-85-036 

SCIENTISTS  AND  SCIENTIFIC  ORGANIZATIONS  23  October  1985 


CHINA'S  FIRST  SYMPOSIUM  ON  TECTONOPHYSICS  HELD 

Beijing  DIZHI  LUNPING  [GEOLOGICAL  REVIEW]  in  Chinese  Vol  31,  No  4,  Jul  85 
p  377 

[Text]  In  order  to  promote  the  interaction  between  the  subbranches  of  earth 
sciences,  a  sumposium  on  tectonophysics  and  earth  crust  deformation  was  held 
in  Beijing  on  10-13  March  1985.  The  symposium  was  sponsored  jointly  by  the 
Tectonogeology  Committee  of  the  China  Geology  Society  and  the  Seismogeology 
Committee  of  the  China  Seismology  Society.  Tectonophysics  is  an  interdisci¬ 
plinary  science  based  on  and  developed  from  traditional  geology  and  geophysics 
in  the  last  few  decades.  Its  main  scope  of  study  is  the  physical  nature  of 
structural  deformations  on  various  scales  of  magnitude.  It  is  related  to 
geodynamics  and  geomechanics  but  has  a  different  emphasis.  Both  provide 
scientific  basis  for  mineral  resources  development  and  for  minimizing  earth¬ 
quake  damages . 

In  the  meeting,  more  than  90  papers  were  presented  by  100  participants  repre¬ 
senting  24  units  from  the  State  Seismological  Bureau,  the  Chinese  Academy  of 
Sciences,  institutes  of  geology,  coal,  petroleum,  chemical  engineering,  and 
mining,  and  institutes  of  higher  learning.  Topics  discussed  in  the  symposium 
include  structural  deformation  on  various  magnitude  scales,  the  inception  and 
the  temporal  and  spatial  distribution  of  earthquakes,  and  the  mine  formation 
mechanism.  Presentations  are  based  on  field  survey,  earth  stress,  deformation, 
and  seismological  observations,  and  physical  and  mathematical  modelling.  This 
symposium  demonstrated  the  recent  developments  and  trend  of  tectonophysics  in 
China.  Notably  features  are:  (1)  Geophysical  and  geodetic  survey  of  tech¬ 
nology  has  greatly  improved  the  depth,  breadth,  and  resolution  of  field  obser¬ 
vations,  and  helped  to  quantify  the  study  of  structural  deformation.  Papers 
presented  in  the  symposium  include  studies  of  the  structural  movement  at 
Wenanbaxian,  the  horizontal  deformation  field  and  stress  field  of  the  Beijing 
region,  the  radial-arcuate  orthogonal  seismic  zone  in  northern  China,  small 
scale  convection  in  the  upper  mantle  and  the  stress  distribution  of  the  sub- 
duction  zone  of  plates.  These  studies  made  use  of  geological  deformation  and 
stress  data  and  seismology  and  satellite  gravity  data.  (2)  The  wide  appli¬ 
cation  of  mathematical  mechanics  principles  and  mathematical  modelling  shows 
that  tectonophysical  research  in  China  has  made  considerable  progress  toward 
quantification.  Numerical  and  analytical  methods  have  been  applied  to  the 
study  of  fracture  propagation,  -shaped  fracture,  scoop-shaped  fault,  tri¬ 
angular  structure  and  microstructures,  and  earthquake  threats  and  mine 
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formation  mechanisms.  Most  significant  development  is  the  application  of  the 
finite  element  method  and  elasto-plastic  modelling  and  viscoelastoplastic 
modelling.  Such  applications  have  elevated  the  mathematical  modelling 
tectono-deformation  into  a  new  level.  (3)  More  and  more  attention  is  given 
to  the  uniqueness  and  complexity  of  structural  deformations  with  respect  to 
the  medium  environment  and  spatial  and  temporal  distribution  and  to  t 

importance  of  high  temperature  and  high  pressure  rock  mechanics.  As  a  result 
importance  ox  ng  f  developed  triaxial  experimental  techniques  for 

°  cSCunder  high  ;reSsura^d  «  high  temperature.  Triaxial  testing  apparatus 

has  been  constructed  using  liquid,  gas  or  solid 

ThCsolid^medium3 triaxial6 testing  equipment  has  also 

study  of  the  strongly  constrictive  and  rigid  environment  of  the  interior  of 
tS  earth  Progresses  in  the  investigation  of  fault  material,  porosity  pres- 
sure  and  creep!  fracture  and  friction  in  rocks  and  developments  in  acoustic 
emission  monitoring  and  localization,  Reiser  effect,  laser 

j  pccentric  model  studies  hhve  also  attracted  attention.  (4)  A  uniqu 
^r^cStteitonophysical  studies  in  China  is  the  combination  of  the 
macroscopic  phenomena  and  the  microscopic  mechanisms.  A  number  of  papers 
dealt  with  the  deformation  damage  and  stress  estimate,  the  study  ofmedu 
nonuniformity,  and  the  causes  for  the  formation  of  low  speed  and  low  resi 
tance  layers  in  the  earth  crust.  These  studies  have  made  use  of  microstruc 
tural  analysis,  dehydration,  phase  change,  and  partial  fusion  at  high  tempera¬ 
ture  and  high  pressure,  and  thermodynamic  calculations. 

The  participants  of  this  symposium  include  the  old,  the  middle-aged  and  the 
young  generations.  Presentation  made  by  many  young  researchers  has  bro^bt 
I!n-Uo  the  meeting.  The  meeting  was  well  organized  and  proceeded  in 
a  fait  and  active  pace!  in  order  to  promote  the  development  of  tectonophysics 
In  SSiTlid  S  facilitate  academic  exchange  and  cooperation,  the  representa¬ 
tions  recommended  to  establish  a  tectonophysics  specialty  group  under  th® 
joint  leadership  of  the  Seismogeology  Committee  of  the  v  xeeom- 

and  the  Tectonogeology  Committee  of  the  China  Geology  o  y 
mended  a  second  Tectonophysics  Symposium  in  the  fourth  quarter  of  19  7. 
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QIU  Xinfang  [6726  0207  5364] 

WANG  Xu  [3076  2485] 

ORG:  None 

TITLE:  ''Mutagenicity  Test  of  ABS  and  AS  and  Their  .Inhibitory  Effect  On 
Mitomycin  C  Induced  SCE  Frequency" 

SOURCE:  Shanghai  FUDAN  XUEBAO  (ZIRAN  KEXUE  BAN)  [JOURNAL  OF  FUDAN  UNIVERSITY 
(NATURAL  SCIENCE)  in  Chinese  Vol  24,  No  2,  Jun  85  pp  121-126 

TEXT  OF  ENGLISH  ABSTRACT: 

Both  alkylbenzene  sulfonate  (ABS)  and  alkyl  sulfonate  (AS)  are  major  active 
components  of  current  synthetic  detergents.  The  potential  mutagenicity  of  ABS  and 
AS  were  investigated  with  Salmonella/microsome  test,  no  significant  or  consistant 
difference  was  detected  between  ABS  (5 — 30pg/ml)  or  AS  (0.01 — 0.5pg/ml)  and 
control.  Also  there  was  no  difference  whether  the  metabolic  activation  system  of  89- 
mix  was  used  or  not.  Treatment  of  ABS  (l~10ng/ml)  or  AS  (0.001 — 0.05pg/ml) 
alone  caused  no  significant  increase  of  SCE  frequency  in  Chinese  hamster  Wg3-h 
cells.  Treatment  of  Wg3-h  cells  with  various  doses  of  Mitomycin  C  (MMC)  from 
0.005  to  0.05pg/ml  caused  a  significant  dose-dependent  increase  of  SCE  frequencies. 

The  increase  was  significantly  inhibited  by  concurrent  addition  of  a  fixed  dose  of  ABS 
(5pg/ml)  or  AS  (O.Olpg/ml).  Furthermore,  Wg3-h  cells  exposed  to  MMC  at  a  fixed 
dose  of  0.025|ig/ml  plus  various  lower  ABS  doses  from  0.01~0.0005pg/ml  or  to  ABS 
at  a  fixed  dose  of  0.005  (xg/ml  plus  various  MMC  doses  from  0.05 — 0.005  pg/ml  all 
showed  a  dose-dependent  inhibitory  effect  on  MMC-induced  SCE.  These  results  indica¬ 
te  that  under  certain  conditions,  ABS  and/or  AS  suppress  the  SCE  frequencies  induced 
by  MMC.  We  suggest  that  ABS  and/or  AS  might  have  some  effect  on  the  protection 
of  human  against  environmental  harmful  factors.  (Paper  received  30  May  84) 
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TITLE:  "Preliminary  Studies  On  the  Preparation  of  a  Ribosomal  Vaccine  From 
Bacteria  S",02"^  and  Its  Antitumor  Activity" 

SOURCE:  Shanghai  FUDAN  XUEBAO  (ZIRAN  KEXUE  BAN)  [JOURNAL  OF  FUDAN  UNIVERSITY 
(NATURAL  SCIENCE)  in  Chinese  Vol  24,  No  2,  Jun  85  pp  133-139 

TEXT  OF  ENGLISH  ABSTRACT: 

This  paper  is  the  first  report  on  the  ribosomal  vaccine  from  Gram  negative  bacteria 
(S-Oa-l)  in  China.  Various  techniques  such  as  ultrasonic  disruption,  differential 
centrifugation  and  gel  filtration  were  used  for  its  isolation.  The  results  showed  that  the 
ribosomal  vaccine  could  stimulate  the  immune  response  of  spleen  in  mice,  and  suppress 
the  proliferation  of  transplanted  S-180  tumor  cells  in  the  peritoneal  cavity  of  mice. 

When  used  in  combination  with  the  incomplete  Freund’s  adjuvant,  the  inhibitory 
effect  of  the  vaccine  became  much  higher  than  that  of  the  vaccine  used  alone.  And  the 
toxicity  of  ribosome  was  lower.  The  results  could  provide  some  facts  for  the  clinical 
preparation  of  antitumor  immunopotentiators.  (Paper  received  on  14  November  84) 
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AUTHOR:  HU  Youliang  [5170  0645  5328] 

WANG  Huifang  [3769  1920  2455] 

ZHOU  Lirong  [0719  7787  5554]  o 

et  al .  .  , 

ORG:  Institute  of  Chemistry,  Chinese  Academy  of  Sciences 

TITLE:  "Determination  of  the  Chemical  Composition  and  the  Crystal  Modifica¬ 

tion  of  the  Ziegler-Natta  Catalyst  Prepared  from  the  Alkoxy  Titanium  Chloride" 

SOURCE:  Beijing  GAOFENZI  TONGXUN  [POLYMER  COMMUNICATIONS]  in  Chinese  No  1, 

Feb  85  pp  54-57  '  / 

TEXT  OF  ENGLISH  ABSTRACT:  The  elemental  analysis  data  and  corresponding  gas 
chromatography  indicated  that  the  empirical  compositions  of  the  intermediate 
compounds  obtained  by  reducing  Ti(OCitH9)Cl3  with  AlEt2Cl  are  TiCOC^Hg )C12  (I), 
and  the  catalysts  obtained  through  the  heat  treatment  of  TKOC^Hg  )C12  by 
means  of  TiCli,  have  the  formula  of  TiCl3  •  0.02  A1C13  •  0.05  (n-CitH9)20  • 

0.02  C^HgO  (II).  The  X-ray  diffraction  pattern  showed  that  (I)  is  an  amorphous 
material  and  (II)  is  S~TiCl3.  (II)  has  almost  the  same  crystal  modification 
(20  (deg)  =  15.13,  33.03  and  51.52)  and  average  size  of  crystallites  (~100  A) 
as  the  Solvay  catalyst  (III)  prepared  from  TiCl4.  (Paper  received  on 
28  October  1983. ) 
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GUAN  Jiayu  [7070  1367  3768]  !- 

QIAN  Renyuan  [6929  0086  0337] 

ZHANG  Pengfa  [1728  2590  4099] 
et  al . 

ORG:  GUAN  and  QIAN  of  the  Institute  of  Chemistry,  Chinese  Academy  of 

Sciences;  ZHANG,  et  al'.-,  of  the  Institute  of  Atomic  Energy,  Chinese  Academy 
of  Sciences 

TITLE:  "Effect  of  the  State  of  Macromolecular  Aggregation  on  the  Uranium 

Nuclear  Fission  Tracks  in  PET  Films" 

SOURCE:  Beijing  GAOFENZ I  TONGXUN  [POLYMER  COMMUNICATIONS]  in  Chinese  No  1, 

Feb  85  pp  62-65 

TEXT  OF  ENGLISH  ABSTRACT:  235U  fission  track  observations  have  been  performed 
for  PET  samples  stretched  and  heat  treated  under  different  conditions.  These 
damage  tracks  may  be  revealed  and  made  visible  under  an  ordinary  optical 
microscope  by  treatment  with  KOH  solution  that  rapidly  and  preferentially 
attacks  the  damaged  material.  As  the  stretching  and  annealing  of  PET  films 
causes  change  in  the  state  of  macromolecular  aggregation  in  the  film,  the 
diffus ion— etching  rate  of  the  KOH  solution  will  be  changed  accordingly. 
Anisotropy  in  the  diffusion-etching  rate  gives  us  anisotropy  in  the  diameter 
of  the  track  cross  section.  (Paper  received  on  11  November  1983.) 
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LIU  Xiaoming  [0491  2556  2494] 

ZHAO  Huashan  [6392  5478  1472;  deceased] 

ORG:  Beijing  Institute  of  Chemical  Technology 

TITLE:  "A  New  Electron  Microscopy  Staining  Method  for  Morphological  Studies 

of  Rubber-Modified  PVC — The  Sulfochlorides-Heavy  Metal  Staining" 

SOURCE:  Beijing  GAOFENZI  TONGXUN  [POLYMER  COMMUNICATIONS]  in  Chinese  No  1, 

Feb  85  pp  66-69 

TEXT  OF  ENGLISH  ABSTRACT:  A  sulfochlorides-heavy  metal  staining  method  is 
developed  for  morphological  studies  of  rubber-modified  PVC  (PVC/CPE,  PVC/EVA, 
PVC/MBS,  PVC/ACR,  etc.).  The  staining  procedure  is  as  follows:  Specimens 
react  with  HS03C1  below  -5°C,  then  are  embedded  in  epoxy  resin  and  mierotomed. 
The  ultrathin  sections  are  stained  with  U02Ac2  or  OsOi,  for  TEM.  This  method 
is  simpler  and  more  effective,  and  the  period  of  staining  is  much  shorter, 
than  those  reported  in  the  literature.  (Paper  received  on  17  November  1983.) 
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AUTHOR:  JIA  Lei  [6328  4320] 

CHEN  Sailu  [7115  6357  3873] 

WANG  Xiugang  [3769  4423  1511] 

ORG:  Institute  of  Chemistry,  Chinese  Academy  of  Sciences 

TITLE:  "Preparation  of  Activated  Carbon  Fiber  (ACF )  Containing  Nitrogen 

Atom  and  Its  Absorption  Behavior" 

SOURCE:  Beijing  GAOFENZ I  TONGXUN  [POLYMER  COMMUNICATIONS]  in  Chinese  No  1, 
Feb  85  pp  70-72 

TEXT  OF  ENGLISH  ABSTRACT:  Polyacrylonitrile  (PAN)  fiber  was  activated  at 
various  temperatures  for  varying  lengths  of  time  in  different  fluxes  of  C02. 
In  this  experiment,  the  best  conditions  for  activation  were  a  temperature  of 
850°C ,  time  period  of  2  hours,  and  C02  flux  of  4  1/min,  which  produced  an  ACF 
having  a  specific  surface  area  of  900-1000  m2/g  and  absorbing  CCI4  60-70  per¬ 
cent  (in  weight).  It  has  been  found  that  the  PAN  ACF  Possesses  the  same 
properties  as  the  viscose  ACF  made  in  Japan  through  the  determination  of 
absorption,  desorption  and  pore  diameter  distribution.  In  addition,  the 
pore  diameter  distribution  is  narrow  and  the  peak  of  distribution  is  about 
30  A.  It  has  a  larger  specific  surface,  enabling  easy  absorption  and  quick 
desorption  as  compared  with  the  conventional  activated  carbon.  It  is  being 
regarded  as  a  new  absorbent.  (Paper  received  on  20  December  1983.) 
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BAI  Mingzhang  (4101  2494  4545] 

SUN  Lijuan  [1327  7787  1227] 
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ORG:  Institute  of  Elemento-Organic  Chemistry,  Nankai  University 

TITLE:  "Study  of  Molybdenum  Complexes  Catalyzed  Copolymerization  of 

Propargyl  Alcohol  and  Propargyl  Chloride" 

SOURCE:  Beijing  GAOFENZ I  TONGXUN  [POLYMER  COMMUNICATIONS]  in  Chinese  No  1, 
Feb  85  pp  73-76 

TEXT  OF  ENGLISH  ABSTRACT:  The  polymerization  of  acetylenes  has  been  studied 
extensively  since  the  polymers  obtained  are  a  kind  of  prospective  material  r 
of  organic  conductor  and  semiconductor. 

It  was  found  that,  in  the  presence  of  molybdenum  complexes,  propargyl  alcohol 
and  propargyl  chloride  copolymerized  with  vigorous  reactions.  The  polymer  I 
obtained  was  insoluble  in  all  solvents  studied,  and  polymer  II  was  soluble 
in  methanol,  chloroform,  acetone,  etc.,  with  a  molecular  weight  of  700-1000. 

.  "t  "  .  '•  #  . 

The  authors  have  tried  using  nine  complexes  of  Ti,  Zr,  Fe,  Co  and  Ni  as  the 
catalyst,  but  no  reaction  has  taken  place.  (Paper  received  on  16  December 
1983.) 
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AUTHOR:  SUI  Chengyun  [7131  6134  0663] 

ORG:  Guangzhou  Research  Institute  of  Non-Ferrous  Metals 

TITLE:  "Optimal  Carrier  Gas  Flow  Rate  in  Inductively  Coupled  Plasma  Atomic 

Emission  Spectrometry" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13 
No  5,  20  May  85  pp  321-324 

TEXT  OF  ENGLISH  ABSTRACT:  It  is  shown  experimentally  that  the  optimal  carrier 
gas  flow  rate  of  ICP-AES  is  determined  by  ICP  power  and  excitation  potential 
of  the  analytical  spectral  lines  and  is  independent  of  the  type  of  nebulizing 
system  used  and  its  performance.  The  range  of  carrier  gas  flow  rate  attainable 
from  the  nebulizing  system  must  match  the  ICP  parameters.  (Paper  received  on 
13  February  1984.) 
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OUYANG  Zhongxue  [2962  7122  1350  1331] 
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YANG  Xi  [2799  3556] 
et  al. 

ORG:  XIE,  OUYANG  and  WEN  of  the  Department  of  Chemistry,  Jiangxi  Univer¬ 

sity;  YANG,  et  al.,  of  Dexing  Copper  Mine,  Jiangxi  Copper  Company 

TITLE:  "Studies  of  Supporter  Extraction  Spectrophotometry.  II.  Determination 

of  Trace  Gold  in  Ores" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13  No  5, 
20  May  85  pp  324-329 

TEXT  OF  ENGLISH  ABSTRACT:  Gold  in  the  sample  solution  is  extracted  with  a 
soft  foam  loading  with  0.20  ml  MIBK,  and  then  color  is  developed  in  a  weak 
acidic  medium  with  4,4'-bis  (dimethylamino)-thiobenzophenone  (TMK) .  The 
resulting  Au-TMK  complex  can  be  measured  directly  by  dual-wavelength  spectro¬ 
photometry.  The  method  has  been  applied  to  the  determination  of  Au  down  to 
0.02  g/t  in  ores  with  satisfactory  results.  (Paper  received  on 
29  February  1984.) 
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SONG  Junfeng  [1345  0193  1496] 

CAO  Yuanzhe  [2580  0337  0772]  ' 

HU  Yinhua  [5170  5593  5478] 

ORG:  Department  of  Chemistry,  Northwest  University 

TITLE:  "Polarographic  Study  of  Metal-4-(2-Pyridylazo)-Resorcinol  (PAR) 

Complex.  III.  Catalytic  Polarographic  Wave  in  Niobium  (V)-Tartaric  Acid- 
PAR-NaC103  System" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13  No  5 
20  May  85  pp  337-339 

TEXT  OF  ENGLISH  ABSTRACT:  In  a  supporting  electrolyte  of  1.6  x  10_5M 
PAR-0. 01M  tartaric  acid-2. 2M  NaCl03  (pH  1.9±0.3)  niobium  (V)  produces  a 
sensitive  catalytic  wave  with  peak  potential  at  -0.90  V  (vs  SCE).  The  wave 
height  is  proportional  to  the  concentration  of  Nv(V)  in  the  range  of  1  ppb 
to  80  ppb.  Ta  in  4000-fold  excess  relative  to  0.02  ppm  Nb  does  not  interfere 
The  method  has  been  used  for  geochemical  samples.  (Paper  received  on 
22  May  1984.) 
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LI  Yunhua  [2621  0061  5478] 

SHI  Huiming  [0670  1979  2494] 

ORG:  Department  of  Chemistry,  Nankai  University 

TITLE:  "Optimization  Design  Investigation  of  Spectrophotometric  Analysis 

by  Weighted  Centroid  Simplex  Method" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13 
No  5,  20  May  85  pp  344-349 

TEXT  OF  ENGLISH  ABSTRACT:  An  improved  simplex  optimization  procedure  based 
on  the  weighted  centroid  method  has  been  developed  which  shows  promise  as  a 
method  for  accurate  and  rapid  search  of  an  optimum  region. 

We  applied  this  procedure  to  a  nonlinear  equation  and  two  color  reaction 
systems.  The  proposed  procedure  has  been  shown  to  yield  the  correct 
experimental  conditions  for  two  systems.  (Paper  received  on  12  June  1984.) 
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TITLE:  "Conductometric  Determination  of  Crown  Ether  Content" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13  No  5, 
20  May  85  pp  349-353 

TEXT  OF  ENGLISH  ABSTRACT:  A  new  method  for  the  determination  of  crown  ether 
content  is  described.  Let  a  Rind  of  crown  ether  dissolve  in  a  certain  amount 
of  low  dielectric  constant  solvent  chloroform  and  add  a  solid  metal  salt 
(such  as  alkalimetal  picrate)  in  excess.  The  chosen  salt  is  hardly  soluble 
in  the  solvent,  and  capable  of  complexation  with  the  crown  ether  for 
determination.  After  the  reaction  equilibrium  is  attained,  the  electric 
conductance  of  the  system,  L,  is  measured.  Then  the  concentration  of  crown 
ether,  R,  can  be  drawn  by  the  equation:  R/L  =  QL  +  P ,  where  Q  and  P  are 
constants.  (Paper  received  on  14  June  1984.) 
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Sciences 

TITLE:  "Determination  of  Nd3+  in  the  Laser  Crystal  Y3Al5012:Nd3+  by 

Magnetic  Balance" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13  No  5, 
20  May  85  pp  354-356 

TEXT  OF  ENGLISH  ABSTRACT:  Magnetic  susceptibilities  of  laser  crystal 
YAG:Nd3+  were  determined  at  10°C  by  magnetic  balance.  The  increased  linearly 
with  the  increasing  concentration  of  Nd203.  Using  the  linear  equation 
Nd203(wt -percent)  =  aXg  +  b,  the  concentration  of  Nd203  in  YAG:Nd3+  can  be 
determined  rapidly  by  the  magnetic  susceptibility  (Xg)  measured.  Solid 
powder  samples  may  be  used  directly  without  decomposition,  separation  or 
other  time-consumption  processes.  (Paper  received  on  25  June  1984.) 
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TITLE:  "Studies  of  a  New  Color-forming  System.  III.  The  SystemHg(ll)-5- 

(3,5-Dichloro-2-Pyridylazo)-2,4-Diaminotoluene-Sodium  Dodecylbenzene 
Sulfonate" 

SOURCE:  Changchun  FENXI  HUAXUE  [ANALYTICAL  CHEMISTRY]  in  Chinese  Vol  13  No  5, 
20  May  85  pp  362-365 

TEXT  OF  ENGLISH  ABSTRACT:  The  optimum  conditions  for  the  color  reaction  of 
Hg(H)-3,5-Cl2-PADAT-SDBS  have  been  studied  systematically.  In  a  tartrate 
buffer  solution  of  pH  3.8  to  5.0,  Hg(ll)  forms  a  red  ternary  complex  with 
3,5-Cl2-PADAT  in  the  presence  of  SDBS .  The  maximum  absorption  is  measured 
at  520  nm  with  molar  absorptivity  being  5.11  x  104  1  mol-1  cm-1.  The  ratio 
of  Hg(ll)  to  3,5-Cl2-PADAT  to  SDBA  in  the  ternary  complex  has  been  found  to 
be  1:2:1.  A  linear  relationship  is  obtained  in  the  concentration  range  of 
2.5  to  70  yg/25  ml.  The  method  is  simple,  rapid  and  accurate,  and  can  be 
applied  to  the  spectrophotometric  determination  of  mercury  in  ores.  (Paper 
received  on  3  January  1984.) 
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TITLE:  "The  Great  Achievements  of  Nanjing  University's  Scientific  Research 

in  the  35  Years  Since  the  Founding  of  the  PRC" 

SOURCE:  Nanjing  NANJING  DAXUE  XUEBAO  (ZIRAN  KEXUE)  [JOURNAL  OF  NANJING 
UNIVERSITY  (NATURAL  SCIENCES  EDITION)]  in  Chinese  Vol  20  No  4,  Dec  84 
pp  605-613 

TEXT  OF  ENGLISH  ABSTRACT:  The  great  achievements  of  this  university  in 
scientific  research  during  the  35  years  since  the  founding  of  the  PRC  under 
the  leadership  of  the  Party  are  described  in  this  paper.  The  university  has 
gone  through  four  stages  of  scientific  research. 

I.  Initial  Period  (1949-1956): 

Most  of  the  teachers  selected  their  own  research  topics  based  on  the  academic 
speciality  of  the  old  teachers.  Although  there  were  not  many  scientific 
research  projects  and  topics,  they  reflected  a  significant  level  of  academic 
research  and  showed  particular  characteristics; 

II.  Period  of  Overall  Expansion  (1956-1966): 

Scientific  research  developed  overall  from  decentralized  and  self-selected 
research  to  unified,  organized  and  planned  research; 

III.  Period  of  Setbacks  (1966-1976): 

Because  of  10  years  of  internal  disorder,  scientific  research  in  this 
university  was  hit  by  an  unprecedented  setback.  Under  extremely  difficult 
circumstances,  may  teachers  persisted  in  their  work  and  did  not  allow  their 
scientific  research  to  completely  stop; 

IV.  Flourishing  Period  (1976-present): 

Under  a  corrected  line,  by  emancipating  the  mind  and  inspiring  enthusiasm,  the 
scientific  research  of  this  university  has  taken  on  a  new  and  flourishing 
look.  A  series  of  gratifying  achievements  have  been  accomplished  by  unifying 
the  methods  of  basic  research  and  applied  research,  and  by  developing  inter¬ 
disciplinary  research. 

When  we  look  to  the  future,  the  responsibilities  and  tasks  are  heavy  and  the 
road  is  long.  This  university  is  determined  to  follow  the  guidelines  of 
educational  reform,  "facing  modernization,  the  world  and  the  future,"  and 
stride  forward  toward  the  great  goals  of  education  and  scientific  research. 
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Structure  Oxalate" 

SOURCE:  Nanjing  NANJING  DAXUE  XUEBAO  (ZIRAN  KEXUE)  [JOURNAL  OF  NANJING 
UNIVERSITY  (NATURAL  SCIENCES  EDITION)]  in  Chinese  Vol  20  No  4,  Dec  84 
pp  786-790 

TEXT  OF  ENGLISH  ABSTRACT:  The  tetragonal  structure  oxalate  was  prepared  by 
coprecipitation  with  (NH4)2C204  solution  as  the  precipitant.  The  processes  of 
formation  and  thermal  decomposition  of  these  compounds  were  studied  by  X-ray 
diffraction,  electron  microscope  observation  and  magnetic  analyses.  The 
following  conclusions  were  obtained:  (1)  The  tetragonal  structure  oxalate 
is  transformed  from  an  a-type  oxalate  in  (NH4)2C204  solution.  (2)  The  thermal 
decomposition  of  these  compounds  is  similar  to  that  of  ol~  and  3-type  oxalates. 
(Paper  was  received  on  10  April  1983.) 
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Ions" 

SOURCE:  Nanjing  NANJING  DAXUE  XUEBAO  (ZIRAN  KEXUE)  [JOURNAL  OF  NANJING 
UNIVERSITY  (NATURAL  SCIENCES  EDITION)]  in  Chinese  Vol  20  No  4,  Dec  84 
pp  791-795 


TEXT  OF  ENGLISH  ABSTRACT:  As  a  new  electrochemical  concentrative  material 
the  carbon-mat  electrode  is  recommended  for  the  purification  of  wastewater 
containing  chromium  ions.  After  chemical  treatment,  the  primary  collecting 
efficiency  of  the  carbon-mat  electrode  reaches  almost  98  percent,  which  makes 
it  ideal  for  the  purification  of  wastewater.  The  experiment  can  be  carried 
out  easily  and  the  operation  is  simple,  successful,  rapid  and  has  good 
reproducibility.  All  experimental  data  were  determined  through  atomic 
absorption  spectrometry. 
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Waveguide  Lens"  ' 

SOURCE:  Tianjin  TIANJIN  DAXUE  XUEBAO  [JOURNAL  OF  TIANJIN  UNIVERSITY] 
in  Chinese  No  1,  Jan  85  pp  20-25 

TEXT  OF  ENGLISH  ABSTRACT:  The  aberration  of  the  spherical  geodesic  optical 
waveguide  lens  and  the  light  intensity  distribution  at  its  focal  plane  have 
been  calculated  by  ray  tracing  and  Fourier  Transformation  techniques.  The 
optimum  aperture  of  this  lens  and  the  minimum  size  of  the  focal  spot  can  also 
be  determined  conveniently  with  our  computer  program.  In  addition,  the 
aberrations  of  the  lens  is  corrected  with  the  twice  ion -exchange  method  inside 
the  sphere.  Both  sphere  and  plane  waveguides  in  our  lens  are  of  a  single 
mode,  and  the  measured  effective  index  ratio  of  their  TE0  mode  is  1.009, 
while  the  focal  spot  of  16  jjm  is  obtained  at  a  relative  aperture  of  0.2. 

The  experiment  results  are  in  good  agreement  with  the  theoretical  analysis. 
(Paper  received  on  21  May  1984.) 


110 


AUTHOR:  SUN  Rend e  [1327  0117  1795] 

SUN  Dong  [1327  2767] 

LU  Dengping  [6424  4098  1627] 
et  al. 

ORG:  All  of  Tianjin  University  :  .v";..-  ‘ 

TITLE:  "The  Application  of  Image  Processing  by  a  Microcomputer  in  Welding 

Electric  Arc  Research" 

SOURCE:  Tianjin  TIANJIN  DAXUE  XUEBAO  [JOURNAL  OF  TIANJIN  UNIVERSITY] 
in  Chinese  No  1,  Jan  85  pp  102-106  i.  ^ 

TEXT  OF  ENGLISH  ABSTRACT:  Image  processing  by  a  microcomputer  is  a  new 
development  in  computer  applications .  Welding  electric  arc  is  a  heat  source 
with  high  temperature  and  strong  light,  and  is  hard  to  investigate  by 
ordinary  methods.  This  paper  describes  how  the  temperature  field  of  welding 
electric  arc  was  measured  by  .  the  use  of  pictures  taken  by  an  ultrared  camera 
and  a  microcomputer.  The  interrelationship  between  the  welding  electric 
current  and  welding  temperature  field  was  obtained  when  D.C.  TIG  welding 
was  applied  to  stainless  steel.  Under  the  arc  with  an  Ar+N2  protecting  gas, 
a  clear  image  and  arc  structure  was  produced  by  the  false  color  method.- 
(Paper  received  on  5  April  1984.) 
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TITLE:  "Research  On  the  Properties  of  Blue  Tungsten  Oxide  Prepared  By 

Intraraolecularly  Reducing  Method" 

SOURCE:  Shanghai  FUDAN  XUEBA0  (ZIRAN  KEXUE  BAN)  [JOURNAL  OF  FUDAN  UNIVERSITY 

(NATURAL  SCIENCE)  in  Chinese  Vol  24,  No  2,  Jun  85  pp  141-146 

TEXT  OF  ENGLISH  ABSTRACT: 

A  new  method  for  preparing  blue  tungsten  oxide  was  d«scribed.  With  APT  as  raw 
material,  named  IRT  (Intramolecularly  Reducing  Tungstate)  was  synthesized.  Then, 

IRT  was  heated  in  Nt  atmp.  at  620  ±  30°C,  for  0.5  hour  to  give  blue  tungsten  oxide 
and  the  blus  tungsten  oxide  was  investigated  by  X-ray  diffraction,  I.  R.  spectroscopy, 

S.  E.  M.  and  chemical  analysis. 

The  crystal  of  the  blue  tungsten  oxide  was  simple,  WO,J0.  Purity>99.99%  and 
surface  area  22~30m*/g.  Other  properties,  K+-absorption  and  true  density  were 
reported  too. 

The  blue  tungsten  oxide  prepared  by  this  way  is  an  ideal  intermediate  material  for 
manufacturing  products  of  tungsten.  (Paper  received  on  3  March  84) 
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SOURCE:  Shanghai  FUDAN  XUEBAO  (ZIRAN  KEXUE  BAN)  [JOURNAL  OF  FUDAN  UNIVERSITY 
(NATURAL  SCIENCE)  in  Chinese  Vol  24,  No  2,  Jun  85  pp  147-152 

TEXT  OF  ENGLISH  ABSTRACT: 

1 1-Tongsto-niobo-phosphoric  acid  was  prepared  by  mixing  proper  amounts  of 
sodium  niobate,  tungstate  and  phosphate  iD  alkaline  solution,  followed  by  acidification 
with  sulfuric  acid  to  pH  0.95,  and  then  boiling  for  15  minutes.  Isolation  of  the  solid 

was  performed  through  formation  of  ethyl  ether  adduct.  Dilute  sulfuric  acid,  0.2 M, 
may  be  used  for  further  purification  through  recrystallization. 

U.  V.  absorption  spectra,  I.  R.  spectra  and  DTA  all  indicate  the  resemblance 
between  H^PNbW^O*).* H,0  and  H^PW.A,).*  H.O.  (Paper  received  on  3 
March  19.84) 
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TITLE:  "Computer  Study  of  the  Influence  of  System  Parameters  on  Interface 

Shape  During  Bridgman  Growth"  '  <>  ' 

SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  129-134 

TEXT  OF  ENGLISH  ABSTRACT:  Finite  element  computer  analysis  was  used  to 
investigate  the  influence  of  interface  position,  temperature  gradient,  growth 
rate  and  latent  heat  on  the  interface  shape  during  crystal  growth  by  the 
Bridgman  method.  The  interface  shape  was  found  to  depend  critically  on  the 
distance  from  the  interface  to  the  point  where  the  temperature  gradient  changes 
steeply.  Heat  transfer  brought  about  by  moving  the  crucible,  as  expected,  has 
a  tendency  to  make  the  interface  more  concave.  (Paper  received  on  15  August 
1984.) 
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SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  159-163 

TEXT  OF  ENGLISH  ABSTRACT:  The  lithium  chloro-alumino-borate  glasses  were 
immersed  in  buffer  solutions,  and  the  change  in  pH  value  of  the  solutions  with 
time  was  measured.  The  reaction  process  of  glasses  in  the  solutions  and  the 
change  of  ion  exchange  of  glasses  with  their  compositions  were  analyzed.  It 
was  found  that  the  ion  exchange  ability  of  glasses  has  a  good  relationship  with 
the  ionic  conductivity.  It  follows  that  the  ionic  conductivity  of  glasses 
mainly  depends  on  the  mutual  forces  of  Li+  ions  with  anions  or  anion  groups, 
and  this  supports  the  weak  electrolyte  theory  on  the  ionic  conduction  of  glass. 
(Paper  received  on  12  December  1983.) 
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TITLE:  "Investigation  of  Network  Structure  for  Niobium-Borate  Glasses" 

SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  168-173 

TEXT  OF  ENGLISH  ABSTRACT:  The  network  structures  of  glasses  of  2Nb205’ 
4B203-4K20  (NBK)  and  1 .5Nb205 -4. 5B203 • 3 .5K20 • 0 .5Ti02  (NBKT)  systems  have  been 
investigated  by  IR  spectra  analyses  and  the  radial  distribution  functions 
obtained  from  X-ray  diffraction  measurements.  It  is  found  that  the  coordina¬ 
tion  numbers  of  Nb  and  B  atoms  are  5.8— 5.6  and  3  or  4  respectively. 

Experimental  results  indicate  that  the  average  bond  length  of  the  Nb-0  bond 
is  1.88-la92  A,  and  the  distances  between  Nb-B,  0-0  and  Nb-Nb  are  3.11,  3.17 
and  3.85(A)  respectively.  Based  on  infrared  spectroscopic  analyses,  the  peaks 
for  the  bond  vibration  of  glasses  are:  1440  cm-1  and  1220  cm  1  (0-B=), 

1040  cm-1  or  920  cm-1  (0-B=) ,  1330  cm-1  (Nb-O-B)  and  870  cm  1  ,  730  cm  (O-Nb). 
In  other  words,  the  characteristic  absorption  of  B04  and  Nb06  coexist  with 
B03  in  the  network  structure.  In  this  case,  the  fraction  N4  of  B04  depends 
on  the  contents  of  K20  and  Nb20s.  This  is  similar  to  the  results  obtained 
in  the  system  K20-B203  by  P.J.  Bray  using  the  NMR  method. 

A  model  of  the  NBK  glass  network  structure  has  been  proposed  based  on  the 
analyses  and  experimental  data.  The  network  structure  consists  of  6  and  4 
circular  structural  units  formed  by  NbOg,  B04  and  B03.  (Paper  received  on 
18  May  1984.) 
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SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  178-182 

TEXT  OF  ENGLISH  ABSTRACT:  Fluid  dynamics  for  the  cj)2.2m  astronomical 
telescope  lens  have  been  studied.  Several  problems  are  presented  in  this 
paper,  such  as  the  critical  condition  for  the  vortex  of  the  glass  liquid 
occurring  in  a  single  crucible,  the  calculation  of  the  lowering  of  the 
liquid  surface  and  casting  period,  the  effect  of  the  end  of  a  cylindrical 
pipe  on  the  calculation  and  the  selection  of  the  simulated  criterion  number. 
(Paper  received  on  6  October  1983.) 
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SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  183-189 

TEXT  OF  ENGLISH  ABSTRACT:  An  assumption  that  absorption  of  60Co  narrow  beam 
y-rays  in  glass  or  other  materials  merely  depends  on  the  mass  per  unit  area 
perpendicular  to  the  y-ray  is  made  in  this  paper.  According  to  this 
assumption,  a  relationship  between  the  linear  absorbing  coefficient,  the  lead 
equivalent  and  density  is  deduced.  For  some  glasses  the  given  theoretical 
values  of  the  lead  equivalent  are  compared  with  the  experimental  values .  It 
is  proven  from  the  comparison  that  the  theoretical  values  are  in  agreement 
with  the  experimental  ones.  Therefore,  the  measurement  of  the  lead  equivalent 
and  linear  absorbing  coefficient  of  glasses  can  be  simply  converted  to  the 
measurement  of  their  densities.  In  addition,  similar  problems  with  other 
materials  are  discussed  in  this  paper.  (Paper  received  on  24  May  1984.) 
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SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  190-197 

TEXT  OF  ENGLISH  ABSTRACT:  C.W.  Brindley,  et  al.  (1959)  considered  the 
principal  crystalline  phase  in  the  exothermic  reaction  of  clay  mineral 
kaolinite  at  about  980°C  as  Al-Si  spinel  (2A1  20  3* 3Si02  or  3A1203  -2Si02 ) 
and  pointed  out  that  the  Al-Si  spinel  is  similar  to  that  of  y~A1203,  making 
it  possible  to  mistake  it  for  y-Al 203  in  the  case  of  inaccurate  measurements. 
F.C.  Loughnan,  et  al.  (1981)  found  that  the  kaolinite  flint  clay  could  be 
transformed  into  mullite  and  chris tobalite  during  the  self-burning  of  coal 
in  New  South  Wales,  Australia,  while  in  the  low  temperature  zone  metakaolinite 
occurred.  They  also  considered  the  above  fact  to  have  been  discovered  by 
themselves.  However,  the  formation  of  the  Al-Si  spinel  was  not  ascertained. 

Ma  Jinping,  Ba  Qinglian,  Jing  Xigui  and  the  author  (1982)  observed,  near  a 
highway  in  Hequ  County,  Shanxi  Province,  that  carbonaceous  clay  had  been 
whitened  by  scorching  during  the  self-burning  of  coal  within  the  carboniferous 
coal  series.  Later,  through  the  author's  study,  it  was  determined  that  the 
temperature  of  the  carbonaceous  kaolinite  had  reached  1000°C,  almost  that 
needed  to  form  Al-Si  spinel  near  the  exothermic  peak  of  kaolinite.  It  is  in 
agreement  with  data  obtained  by  F.W.  Brinkley  and  with  that  of  the  exothermic 
peak  of  kaolinite.  The  appearance  of  the  samples  resembles  kaolinite  clay 
heated  to  1000°C.  They  are  white  and  soft  with  micropores,  slightly  rough, 
and  stain  fingers  easily.  They  are  also  dark  in  luster  and  not  viscous. 

The  X-ray  data  from  a  powder  photo  (Jeney  photo,  10  hours)  are  as  follows: 
1.9673  &  (strong),  1.3919  A  (very  strong)  and  1.140  A  (weak)  in  width.  All  of 
the  curves  of  different  thermal  analyses  and  those  of  dehydration  are  straight 
lines.  The  result  of  chemical  analysis  is  similar  to  that  of  the  dehydrated 
kaolinite.  The  content  of  impurities  is  4.69  percent,  A/S  =  0.902.  Granular 
texture  and  micropores  are  shown  under  the  scanning  electronic  microscope. 

The  results  of  infrared  spectrum  analysis  are  the  same  as  those  from  the 
reaction  when  Suzhou's  kaolinite  was  heated  up  to  1000-1100°C  (Chen  Jingyan 
1963),  and  the  1300-960  cm-1  wide  peak  shows  the  curve  of  the  elastic  vibration 
of  (Si-O)-Si;  804  cm-1  shows  the  vibration  of  Si-O-Al,  and  469  cm-1  shows  the 
vibration  of  Si-0. 

The  burned  clay  is  always  associated  with  a  small  amount  of  mullite.  This 
resembles  the  reaction  of  kaolinite  heated  to  1000°C,  but  is  different  from 
halloysite  heated  to  1000°C  (no  mullite).  Sometimes  a  small  amount  of  halloy- 
site  (7  A)  is  formed  by  rain  action.  (Paper  received  on  26  December  1983; 
revision  on  31  July  1984.) 
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in  Chinese  Vol  13  No  2,  Jun  85  pp  217-223 

TEXT  OF  ENGLISH  ABSTRACT:  Sintering  behavior  of  natural  high  purity  dolomite 
and  magnesitic  dolomite  was  studied  with  calcinated  specimens  (D  and  MD)  and 
calcinated-hydrated  specimens  (DH  and  MDH).  The  mechanism  for  the  initial 
sintering  stage  of  MD,  DH  and  MDH  specimens  was  postulated  and  the  activation 
energy  calculated . 

Hydrate  specimens  are  found  to  be  much  more  susceptible  to  sintering  than 
oxides,  with  a  temperature  difference  of  about  400°C  for  the  same  level  of 
densification.  The  susceptibility  of  hydrates  to  sintering  is  expected  to  be 
related  to  the  presence  of  water  vapor,  which  gives  rise  to  some  kinds  of 
rapid  adsorption-desorption  of  OH-  between  two  solid  phases. 

Dense  grog  with  bulk  density  3.30-3.38  g/cm3  and  well-developed  microstructure 
was  obtained  by  sintering  DH  and  MDH  specimens  at  1600°C  for  2-15  hours. 

(Paper  received  on  9  June  1984.) 
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in  Chinese  Vol  13  No  2,  Jun  85  pp  224-233 

TEXT  OF  ENGLISH  ABSTRACT:  Most  chrysotile  asbestos  in  Shimian  County, 

Sichuan  Province,  is  in  the  form  of  slip-tilted  fibers,  but  occasionally 
cross  fibers  are  found.  One  hundred  twenty-one  samples  of  asbestos  fibers, 
including  commercial  fibers,  have  been  collected  and  systematic  miner alogical 
studies  on  the  samples  have  been  made  using  wet  chemical  analyses,  XRD,  DTA, 

TG,  IR,  EM  and  ED.  Physical  properties  of  asbestos,  such  as  tensile  strength, 
elastic  modulus,  volume  resistivity,  specific  magnetic  susceptibility,  specific 
gravity,  heat  resistance,  heat  conductivity,  specific  surface,  chemical 
resistance,  etc.,  have  been  measured,  providing  a  lot  of  data,  plots  and  tables 
for  industrial  application.  Asbestos  fibers  in  this  area  can  be  classified 
into  four  groups  and  seven  types,  according  to  the  impurities  in  fibers. 

In  addition  to  chrysotile  crystallinity,  the  physicochemical  condition  of 
mineral  formation  and  the  influence  of  impurities  on  the  physical  properties 
of  fibers  are  discussed  in  the  paper,  and  the  influences  of  mineral  processing 
and  heat  treatment  on  the  physical  properties  are  also  discussed.  Many 
experiments  show  that,  contrary  to  the  rule  commonly  accepted,  the  tensile 
strength  of  chrysotile  asbestos  in  the  Shimian  mine  is  10-70  percent  more 
than  that  of  the  original  fibers  (150-250  kg/mm2)  after  heat  treatment  at 
380-450°C,  so  that  the  value  of  the  asbestos  in  use  is  increased.  (Paper 
received  on  19  April  1984.) 
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in  Chinese  Vol  13  No  2,  Jun  85  pp  247_249 

TEXT  OF  ENGLISH  ABSTRACT:  Ferroelastic  properties  of  the  NdP5014  crystal  have 
been  determined  by  the  optical  method.  Experiment  results  of  the  micro-optical 
method  show  that  the  ferroelastic  transition  temperature  on  the  (100)  plane, 
is  96±1°C.  Extinction  angles  (6)  between  adjacent  domains  have  been  determined 
as  a  function  of  pressure.  When. a  pressure  of  6  kg/cm3. is  applied  on  the  (100) 
plane,  the  phase  transition  temperature  of  the  crystal  is  117±1°C.  Anomalous 
temperatures  for  birefringence  of  the  NdP5014  crystal  are  clearly  shown. at 
142°C  and  158°C,  using  the  results  of  laser  scanned  curve.  (Paper  received 
on  12  April  1984.) 


122 


AUTHOR: 


YE  Ruilun  [0673  3843  0243] 

WANG  Xinrong  [3769  2450  2837] 

ORG:  Wuhan  Institute  of  Building  Materials 

TITLE:  "Synthesis  of  a  New  Compound — CaeMgCSiO^) 4CI2" 

SOURCE:  Beijing  GUISUANYAN  XUEBAO  [JOURNAL  OF  THE  CHINESE  SILICATE  SOCIETY] 
in  Chinese  Vol  13  No  2,  Jun  85  pp  250-251 

TEXT  OF  ENGLISH  ABSTRACT:  Recently  a  new  compound  was  synthesized  in  the  CaO- 
Mg0-Si02-CaCl2  system.  The  compound  belongs  to  the  cubic  system,  with  regular 
octahedral  shape ,  density  D0}jS  =  3.08  g/em3  and  refraction  index  n  =  1.6910. 
The  crystal  is  stable  in  the  range  of  810-1260°C  and  in  incongruent  at  about 
1265°C.  From  chemical  analysis  its  molecular  formula  is  Ca8Mg(Si04)4Cl2 . 

From  powder  X-ray  diffraction,  the  unit-cell  parameter  is  a  =  15.06  A,  and  the 
crystal  lattice  is  a  face-centered  cubic  lattice.  By  using  the  four-circle 
diffractometer,  it  is  found  that  the  unit-cell  parameter  is  a  =15.0653  K, 
space  group  is  Fd3m,  and  eight  molecules  are  contained  in  a  cell  (Z=8) .  The 
crystal  structure  is  determined  by  a  direct  method  (SHELXTL).  (Paper  received 
on  29  January  1985.) 
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